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INTRODUCTION 


The  T-burner  has  been  employed  [or  sometime  to  define  linear 
response  characteristics  of  solid  propellants  to  fluctuations  in  environ- 
mental conditions.  Although  qualitative  assessment  of  relative  stability 
is  accomplished  without  undue  difficulty,  quantitative  determinations  are 
complicated  by  the  need  for  quantitative  definition  of  the  bosses.  Until 
1972  when  CulickO  ) first  deduced  the  presence  of  acoustic /mean  flow 
interactions  (A/MFI)  from  the  linearized,  one-dimensional  equations  of 
change,  losses  were  considered  determinable  from  properly  executed 
experiments.  However,  Culick's  analysis  introduced  the  flow  turning 
loss  and  vent  gain.  The  former  is  usually  negligible  in  burners  con- 
figured for  pressure  coupled  measurements.  Unfortunately,  the  latter 
is  not.  Because  the  flow  turning  gain  varies  with  mass  flow  rate,  it 
renders  direct  T-burner  loss  measurement  techniques  impotent. 
Therefore,  indirect  "measures"  of  this  term  were  effected  by  including 
it  within  the  statistical  data  correlation  procedure.  Unfortunately,  this 
indirect  procedure  is  not  altogether  practical.  The  lack  of  precision  in 
the  data  bits  combine  with  the  necessity  to  define  three  unknowns  to 
require  a very  large  (and  expensive)  data  base  in  order  to  define  the 
vent  gain  with  confidence. 

It  is  well  known  that  qualitative  features  of  nonsteady,  two- 
dimensional  flows  can  be  determined  by  employing  the  so-called  hydrau- 
lic analogy.  Consequently,  the  purpose  of  this  portion  of  the  program 
was  to  explore  the  vent  region  flow  field  with  the  hydraulic  analogy. 

Since  1969,  steady-state  composite  propellant  combustion  model- 
ing has  been  dominated  by  the  Beckstead,  Derr,  Price  (BDP)  model(2.3). 
Subsequent  development  of  the  basic  model  has  extended  it  to  bidisperse 
AP  propellants  with  aluminum^ ),  polydisperse  AP  propellants  with 
aluminum^5),  and  nitramine  propellants^  “ 8).  The  extensions  to  bidi- 
sperse AP  and  nitramine  propellants  were  made  by  assuming  all  particles 
burn  at  the  same  rate;  the  extension  to  a polydispersion  was  made  by 
"coalescing"  all  particles  above  a critical  diameter  into  a single,  mean 
particle  (those  below  the  critical  diameter  were  assumed  to  react  in  the 
surface  melt);  the  extension  to  nitramines  was  made  by  introducing  a 
surface  melt  criteria  and  different  particle  surface  constraints  with  and 
without  a surface  melt.  A major  difficulty  with  these  analyses  lies  with 
the  fact  that  although  mean  states  are  involved  the  mean  states  are 
postulated  rather  than  deduced.  In  short,  "are  the  mean  states  employed 
correct?"  In  1973  Glick^»  *®)  presented  a "complete"  statistical 
analysis  of  additive  free  composite  propellant  combustion  that  eliminated 
the  need  to  select  a mean  state  because  all  possible  states  were  explicitly 
included  in  the  analysis.  This  approach  appeared  to  be  a viable  alterna- 
tive to  the  mean  state  schemes  and  included  true  mixed,  polydisper6e 
oxidizer  capabilities.  The  purpose  of  the  combustion  modeling  aspect 
of  this  program  was  to  develop  the  "complete"  statistical  approach. 


TECHNICAL  DISCUSSION 


( OMRUSTION  MODELING 


Int  roduction 


A revolution,  wrought  by  combat  pilot's  demands  for  lov  missile 
exhaust  signature,  has  recently  occurred  in  propellant  formulations  for 
many  tactical  applications  of  solid  propellant  rockets.  This  demand  for 
low  signature  has  created  demands  for  new  binders,  new  oxidi/ers, 
new  additives,  and  new  formulations.  Moreover,  as  significant  amounts 
of  condensed  phases  are  verboten  because  they  enhance  visible  signature, 
damping  of  acoustic  disturbances  in  the  motor's  chamber  i.  decreased 
relative  to  that  for  a similar  formulation  with  condensed  phases  in  the 
products.  Consequently,  combustion  instability  has  returned  as  a real 
and  frequent  problem.  Unfortunately,  there  is  neither  a sufficient  base 
of  experience  nor  adequate  understanding  of  combustion  phenomena  (this 
extends  to  the  common  AP/hydrocarbon  propellants  too)  to  solve  low 
signature  propellant  formulation  problems  expeditiously  if  instability  is 
a factor.  As  a result,  the  propellant  formulator  laces  situations  in 
w'hich  rate,  exponent,  temperature  Sensitivity,  physicals,  processing, 
and  energetics  are  constrained*  and  stability  is  demanded.  In  the 
ballistics  arena  the  formulators  "weapons"  are  strand  burner,  T-burner, 
particle  collection  bomb,  and  ballistic  test  motor.  Evidence  that  these 
"weapons"  can  prevail  exists  in  every  operational  motor.  Unfortunately, 
as  time  progresses,  the  competition  becomes  ever  fie-cer  as  new  ingred- 
ients and  tighter  constraints  enter  the  fray.  Consequently,  our  past 
successes  may  be  more  preliminary  than  main  events.  In  summary, 
real  formulation  problems  exist  with  low  signature  propellants  if  stability 
requirements  are  imposed. 

The  overall  objective  of  this  work  is  to  construct  -,n  analytical 
model  describing  steady-state  combustion  of  non-metalli  ed,  polydisperse 
composite  propellant.  This  model  can  then  be  employed  within  a computer 
aided  propellant  formulation  strategy  to  assist  in  the  development  of 
propellants  capable  of  meeting  the  challenging  constraints  imposed  by 
today's  weapons  "systems.  Basically,  the  parameters  in  the  model  would 
be  defined  with  a small  baseline  of  data  for  the  oxidizer,  binder,  additive 
system  under  consideration.  With  a viable  model  this  "baselined  model" 
would  he  capable  of  accurately  predicting  ballistic  parameters  within  the 
baseline  propellant(s)  family.  Consequently,  either  the  propellant  formu- 
lator could  interact  directly  with  the  baselined  model  to  satisfy  design 
objectives  or  non-linear  optimization  techniques  could  be  employed  with 
the  baselined  model  to  define  an  optimal  (to  some  criteria)  formulation 
for  the  design  constraints. 

It  is  very  important  to  note  that  a steady-state  combustion  model 
can  also  define  nonsteady  (linear  and  nonlinear)  combustion  characteristics 


,:  As  might  be  expected,  the  lower  bound  of  these  constraints  often 
represents  the  state-of-the-art! 


in  the  low  and  mid-frequency  regimes  through  the  Z-N  procedure. 
Therefore,  the  aforementioned  computer  aided  propellant  formulation 
strategy  can  include  instability  constraints  as  well  as  steady-state  ones. 
Indeed,  extension  of  the  model  developed  herein  into  the  nonsteady 
regime  via  the  Z-N  procedure  is  currently  underway  at  Purdue  Univer- 
sity under  AFRPl.  funding. 

It  would  be  extremely  satisfying  to  state  that  the  modeling  aspect 
of  this  program  has  satisfied  all  of  its  objectives.  However,  it  hasn't  . 
Indeed,  we  have  reversed  our  field  in  the  combustion  modeling  area 
and  are  now  proceeding  in  a direction  considerably  different  from  the 
original  one. 


Fundamentals  of  Statistical  Combustion  Modeling 


The  "clay"  of  the  propellant  formulator  is  heterogeneous  solid 
propellant.  These  propellants  are  admixtures  of  solid  particles 
(oxidizers  and  additives)  and  binder  (polymer  and  additives)  and  the 
total  solids  content  is  usually  as  large  as  possible  because  of  energetics 
constraints.  Therefore,  the  solid  propellant  is  literally  a packed  bed 
of  polydisperse  particles  "filled"  with  binder.  The  admixture  is  most 
often  formed  by  blending  components  in  high  shear  mixers.  Therefore, 
it  is  expected  that  the  packing  is  random.  This  expectation  is  supported 
by  the  isotropy  shown  by  macroscopic  properties  (ballistic,  physical). 
However,  packing  fraction  measurements  indicate  a higher  percentage 
of  theoretical  (an  ordered  packing)  than  is  apparently  achievable 
through  complete  randomization.  At  first,  these  bits  of  evidence 
appear  to  conflict;  however,  harmony  is  restored  by  assuming  that  the 
solids  structure  consists  of  ordered,  multiparticle  aggregates  that  are 
themselves  randomly  arranged.  Therefore,  there  is  order  on  the 
scale  of  particle  size,  but  long  range  disorder.  It  is  interesting  to 
note  that  nature  has  chosen  precisely  this  arrangement  for  liquid 
structure  (molecules  represent  short  range  order  but  liquid  is  disorder- 
ed). Therefore,  this  structure  is  common  in  nature. 


o«oe*ec 


The  burning  surface  traverses  the  solid.  Therefore,  the  burning 
surface  must  reflect  characteristics  of  the  solid.  Although  the  chemical 
discretization  wrought  by  the  heterogeneous  structure  can  be  softened 
by  the  presence  of  surface  melts,  there  appears  little  possibility  that 
the  random  discreteness  wrought  by  the  solids  internal  structure  can 
be  completely  washed  out.  Therefore,  a random,  discrete  structure 
is  expected  on  the  burning  surface.  This  means  that  all  states  of  all 

particles  will  be  present  on  the  burning  surface.  By  all  states  it  is  r-z^TIT'  ■'  !'•  ' 

letTO^^^Tj 
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meant  that  there  is  a sequence  of  states  between  the  state  of  a particle 


t 


Ct 


as  it  emerges  at  the  burning  surface  and  the  state  of  that  particle  as 

it  "leaves"  the  burning  surface  (either  consumed  or  ejected)  and  that  

,,  ....  . , setaoeuce 

every  one  of  these  sequential  states  will  be  represented.  SeTwroo  SwetoCTJCE 

*The  Z-N  procedure^*  enables  one  to  utilize  a steady-6tate  model  Couvjm»TioO 

to  predict  nonsteady  behavior  as  long  as  the  reactive  zones  behave  in  a 
quasi-steady  fashion. 

**No  one  has  solved  the  theoretical  problem  of  the  maximum  packing  fraction 
obtain  able  with  a random,  polydispersion  of  spheres.  Therefore,  we  don’t 
know  if  theory/experiment  are  harmonious  or  not. 
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* 

The  hr6t  thing  to  notv  is  that  the  state  of  a particular  particle 
on  the  burning  surface  is  always  changing.  Therefore,  heterogeneous 
propellant  combustion  is  never  steady-state  in  the  classical  sense  that 
nothing  varies  in  time.  However,  rocket  motor  experience  clearly 
demonstrates  that  steady-state  condition  can  be  closely  approximated. 
Since  the  whole  is  the  sum  of  its  parts,  this  demonstrates  the  following: 
(11  the  motor  follows  the  mean  behavior  of  all  particles  inhabiting  the 
turning  surface  and  (b)  the  mean  phenomena  can  be  stationary  in  time. 
In  other  words,  ballistic  performance  of  a rocket  motor  and  ballistic 
parameters  of  heterogeneous  propellant  depend  upon  ensemble  means 
of  the  deflagrating  surface's  microstates  ! Therefore,  heterogeneous 
p r op ellant  combustion  phenomena  is  a stochastic  phenomena.  C ons e - 
cpu ently,  successful  attacks  on  heterogeneous  propellant  phenomena 
must  come  from  a combination  of  stochastic  and  more  conventional 
methods  of  analysis.  Although  interest  is  centered  herein  on  combus- 
tion phenomena,  this  philosophy  must  also  hold  for  both  processing  and 
physical  properties. 


. TCS. 


r, 


i-b 


Ballistic  parameters  represent  ensemble  means  of  microstates, 

In  practical  propellants  the  particulates  are  polydispe rse . Therefore, 
the  burning  surface's  microstates  will  be  spread  among  a distribution 
of  distinguishable  particles.  Ensemble  averaging  can  be  accomplished 
at  roughly  three  levels.  The  crudest  (and  hence  the  least  realistic)  is 
to  simply  average  all  microstates  into  a single  mean  state.  This  is  the 
one-dimensional  approach.  Since  this  means  that  the  multiplicity  of 
distinguishable  microstates  are  "crammed"  into  a unique  mean,  this 
model  can  predict  only  gross  phenomenological  characteristics  because 
all  information  pertaining  to  the  diversity  of  the  microstructure  has 
been  lost  in  the  averaging  process.  Therefore,  these  models  are 
inherently  incapable  of  predicting  detailed  particle  size  effects.  All 
one -dimensional  combustion  models  are  of  this  type^^  ~ The  next 

level  of  "averaging"  involves  averaging  of  the  microstates  associated 
with  distinguishable  particles.  In  this  "petite  ensemble"  method  diver- 
sity due  to  distinguishable  particles  appears  explicitly  while  diversity 
due  to  a distinguishable  particle's  microstates  appears  implicitly. 

Therefore,  "petite  ensemble"  models  can  predict  particle  size  effects.  , 

The  statistical  combustion  model  developed  under  this  program  and  that  ^ ' 

of  Miller,  et  al.  (20»  2 1 > fall  into  this  category.  The  statistical  combus - 
tion  models  of  Cohen,  et  al.  W and  Sammons(5)  lie  somewhere  betw'een 
this  and  the  former  category  because  they  employ  an  equal  rate  hypothesis. 

That  is,  these  models  "recognize"  distinguishable  particles  but  presume 
that  the  distinguishable  particles  all  have  the  same  regression  rate.  The 
ultimate  reality,  because  it  is  what  nature  does  comes  from  recognizing 
all  microstates  in  a grand  ensemble.  Glick^*  has  partially  formula- 
ted a combustion  model  of  this  type. 


lTa,re\  molt, 
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It  is  important  to  note  that  computational  effort  rapidly  escalates 
with  complexity.  At  the  one-dimensional  level  only  one  flame  structure  "■* — 

must  be  computed;  at  the  "petite  ensemble"  level  one  must  compute  as 
many  flame  structures  as  there  are  distinguishable  particles;  at  the  grand 


*Devote  the  "state  of  a particular  particle"  by  its  "microstate". 
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ensemble  level  as  many  flame  structures  as  there  are  microstates  must 
be  computed.  Moreover,  at  the  grand  ensemble  level  each  microstate  is 
always  nonsteady  and  there  is  no  reason  (other  than  mathematical  necessity) 
to  assume  the  condensed  phase  is  homogeneous.  Thus,  at  the  grand 
ensemble  level  one  is  faced  with  a virtually  overwhelming  computational 
burden.  On  the  other  hand,  "petite  ensemble"  averaging  provides  huer- 
istic  reasons  for  treating  a distinguishable  particle's  flame  structure  as 
ensemble-steady  and  the  condensed  phase  as  homogeneous.  Consequently, 
it  appears  that  the  "petite  ensemble"  approach  offers  the  most  quality 
information  per  unit  cost.  For  this  reason  it  will  prevail  for  a consider- 
able period  of  time.  This  is  the  theoretical  approach  followed  in  this 
work;  this  is  the  reason  that  effort  on  the  grand  ensemble  approach  was 
discontinued. 


Because  of  the  random,  discrete  structure  of  heterogeneous 
propellants,  the  burning  surface  possesses  a random,  discrete  structure. 
Attention  is  focused  herein  on  this  structure  and,  in  particular,  at  the 
basic  "element"  of  this  structure:  an  oxidizer  particle  and  its  "binder". 
These  fuel  surfacer/oxidizer  particle  pairs  are  the  "fundamental  parti- 
cles" of  statistical  combustion  modeling . Because  of  the  random  struc- 
ture of  the  solid,  the  population  of  fuel  surface /oxidizer  particle  pair 
microstate6  for  any  distinguishable  pair  will  be  random. 


/-  fu«l 

Fuf<-  SutfUcs  (ou 


The  purpose  of  statistical  analysis  is  to  relate  mean  behavior  to 
individual  microstate  behavior.  Since  the  whole  is  the  sum  of  its  parts, 
this  is  accomplished  by  summing  all  microstates  on  the  burning  surface. 
Consider  the  control  volume  sketched  in  Figure  1.  For  this  control 
volume  mass  and  energy  conservation  (neglect  kinetic  and  potential 
energy)  give^^ 


Jl 


^ ^ • As  - ^dv/tk  » Ae /At  ^ pALi*  At  + ^ A! 


- 
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where  the  fact  that  l <j.  I - O on  S^,  has  been  employed.  Assume 
that  S*  moves  such  that  the  volume  in  the  control  volume  is 
stationary  "}  , Then  to  an  excellent  approximation 

n yu  *«*.(.*)  . The  stored  energy  E is 
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Since 


"t 


anti  h = h 


the  moving  coordinate  system. 


' dS  3 % Q S 


Eq.  (1)  becomes,  in 


<h) 


c 


With  this  relation,  Eq . (3),  and  the  fact  that  V*  V(t),  Eq . (2)  becomes 

(.5) 


^edjjisVJt  + ^(-6*  - 4:,„W 

4 4 • i 
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Equations  (4)  and  (5)  are  "integral"  equations  that  constrain  the 
phenomena. 


Assume  that  there  are  Q microstates  (j  = 1,Q)  on  Sa 
Fqs.  (4)  and  (5)  become 


^ r 

rp L - 2 i dS  /S* 

AS 


. Then 


ct; 


and 


?,  I - A,*  j ™ \ =< 
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(.7) 


On  the  individual  microstate  level  q , hc  s,  m"  are  functions  of 
time  and  p , e are  functions  of  ^ and  t because  of  temperature 
gradients  and  the  discrete  structure  of  the  condensed  phase. 
Consequently,  if  all  distinguishable  microstates  are  to  be  "counted" 
(grand  ensemble  method)  at  least  Q simultaneous  ordinary  differen- 
tial equations  must  be  solved  to  define  r.  Since  Q is  a very  large 
number,  this  approach  is  currently  computationally  impossible. 

There  are  many  more  microstates  than  distinguishable  particles. 
Therefore,  to  reduce  the  computational  burden  define  a suitable  mean 
microstate  for  each  distinguishable  particle  and  then  "count"  distinguish- 
able particles.  This  is  the  petite  ensemble  approach.  Let  QD  be  the 

...  ....  U 


number  of  particles  and  N 


the  number  of  microstates  for  the 


distinguishable  particle.  Then  Q * and 

subdivision  of  microstates  Eqs.  (6)  and  (7)  become 


Fp  • I ( Z‘ 


:th 
With  this 
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In  both  equations  the  inner  sum  represents  the  sum  over  all  micro- 
states associated  with  the  j1^1  distinguishable  particle.  Since  these 
microstates  are  similar  to  those  of  an  appropriate  monodisperse 
propellant,  the  inner  sums  represent  ensemble  averages*  for  a 
sequence  of  monodisperse  pseudo-propellants.**  In  this  fashion  one 
is  lead  quite  naturally  to  the  concept  that  the  combustion  phenomena 
of  propellants  with  mixes,  polydisperse  oxidizers  can  be  related 
mathematically  to  burning  rates  from  monodisperse,  psuedo-propellants 
if  the  fuel  surface /oxidizer  particle  pairs  act  independently.  In  othe  r 
words,  a complex  poly -phenomena  can  be  reduced  to  a sequence  of 
simpler  mono-phenomena;  the  well  known  divide  and  conquer  strategy. 

Other  advantages  accrue  with  the  petite  ensemble  strategy. 
Because  ensemble  averages  do  not  fluctuate  randomly  in  time  (the 
summation  eliminates  the  "noise"  due  to  the  randomness)  (a)  only 
temporal  fluctuations  coherent  with  environmental  fluctuations 
need  be  accounted  for  and  (b)  all  monodisperse  pseudo-propellant 
condensed  phase  thermophysical  properties  and  thermal  wave  thick- 
nesses are  those  for  the  bulk  propellant  ^ sizes  are  small  compared 
to  the  thermal  waves  thickness.  The  fact  that  thermal  wave  thicknesses 
are  equivalent  for  all  pseudo-monodisperse  propellants  implies  that 
lateral,  conductive  energy  transport  in  the  condensed  phase  will  be 
small. 


Application  of  the  mean  value  theorem  for  integrals  to  Eqs.  (8) 
and  (9)  yields 


r ft. 
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♦Term  these  ensemble  averages  petite  ensemble  averages  to  show  that 
they  apply  to  a particular  distinguishable  particle. 

♦♦This  is  really  a mathematical  construction;  the  pseudo  prefix  denotes 
that  these  monodisperse  propellants  can  not  exist  in  nature. 

♦ ♦♦Recall  that  the  pseudo-propellant  is  a mathematical  construct.  As  the 
condensed  phases  structure  is  random,  a petite  ensemble  psuedo- 
propellant  mean  for  the  condensed  phase  must  include  all  possible 
arrangements  for  those  particular  distinguishable  particles.  The 
mean  of  this  sum  is  the  mean  of  all  possible  arrangements.  In  other 
words,  the  bulk  propellant. 
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where  the  bar  over  denotes  suitable  mean  values.  The  integral  energy 
equation  is  obviously  satisfied  when  each  summand  is  identically  zero. 
Since  the  summand  augmented  by  lateral  cf  ndensed  phase  energy 
transfer  is  the  mean  energy  equation  for  a fuel  surface /oxidizer 
particle  pair  and  lateral  energy  transport  is  small  in  the  condensed 
phase,  any  valid  solution  for  the  mean  microstate  of  a monodisperse 
p seudo -propellant  can  be  employed  with  Eq . (10)  to  compute  the  mean 
rate  of  a polydisperse  propellant. 
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If  the  planar  statistics  of  the  burning  surface  are  stationary,differentia- 
tion  of  Eq.  (12)  yields 

"b?  f'bC.)  * pt  / <>c  )]  ^ 03) 

where  ( ) denotes  any  environmental  variable  (p,  T„,  X flow,  etc.). 

Thus,  with  some  manipulation  and  the  appropriate  definitions 
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These  relations  show  that  aU  important  ballistic  parameters 
can  be  computed  Iron  pseudo. monodisperse  values.  They  also  show 
that  if  the  pseudo-monodisperse  propellant  mean  properties  ( ).  are 

bounded  the  polydisperse  mean  property  ( ) lies  within  these  bounds 


The  entral  task  in  the  petite  ensemble  approach  is  to  determine 
appropriate  monodisperse  pseudo-propellant  means.  Since  the  mean 
microstate  must  satisfy  the  constraints 
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These,  together  with  the  phenomenological  fuel  surface  /oxidizer  particle 
pair  combustion  model,  represent  the  major  instruments  for  defining 
appropriate  means.  For  this  purpose  these  integrands  can  be  inter- 
preted as  a value  for  a particular  microstate  and  dS  as  the  surface  area 
occupied  by  particles  with  that  microstate.  Assume  a reasonable 
functional  form  for  m"  • (t).  Then,  since  temporal  and  spatial  micro- 
states are  related  by  the  ergodic  surmise  (stationary  phenomena),  the 
mass  flux  from  any  microstate  can  be  determined.  The  number  of 
specific  microstates  is  defined  by  assuming  a burning  surface  topo- 
graph.. The  oxidizer  particle  size  distribution  defines  the  number 
density  of  j1*1  oxidizer  particles  (dNj/dV)  as  only  those  j1"  particles 
lying  within  ±Dj/2  from  the  burning  surface  can  intersect  it  and  all  ^ 
"depths"  of  intersection  are  equally  probable,  the  *' 

number  of  particles  in  any  specific  microstate  can  be  determined.  By 
combining  these  two  steps  mean  values  can  be  defined  in  terms  of 
parameters  in  the  forms  assumed  for  mass  flux  and  surface  topography. 

There  now  appears  to  be  several  ways  to  proceed.  They  depend  (» 
upon  whether  Eq.  (16)  is  employed  as  a relation  for  the  mean  state  or 
as  a tool  to  assist  in  defining  the  mean  state.  In  the  first  approach  the 
mean  state  is  defined  by  requiring  that  the  mean  microstate  be  an 
accessible  microstate  and  then  Eq . (16)  employed  with  the  physio- 
chemical  combustion  model  to  define  , etc.  Another  approach 

would  be  to  employ  certain  aspects  of  the  p^ysiochemical  combustion 
model  (i.  e.  , surface  pyrolysis  law(s)]  to  define  hc  g in  terms  of 
mp.  j (t)  and  then  Eq . (16)  to  define  q^  . By  requiring  the  mean 
microstate  to  be  an  accessible  microstate  m"  ,,  etc.  could  then  be 
defined.  ‘ 
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The  above  is  vague;  it  is  difficult  to  be  general  here  without 
living  vague.  An  important  point  is  that  surface  topography,  mjj  j(t), 
and  the  mean  state  cannot  all  be  arbitrarily  specified  if  a consistent 
model  is  desired. 

The  discussion  above  has  concentrated  on  the  situation  where 
fuel  su rface /oxidize r particle  pairs  are  independent  and  the  condensed 
phase  preheat  zone  is  much  thicker  than  the  largest  oxidizer  particle. 
Both  assumptions  are  untrue  to  some  degree.  Does  this  invalidate  the 
approach?  No,  it  doesn't,  but  it  does  make  it  more  complex  because 
inter-particle  interactions  must  be  considered  to  relax  these  assump- 
tions. Fortunately,  on  a deflagrating  surface  there  is  a finite 
"memory"  because  the  combustion  process  literally  "burns  up"  the 
"distant"  past.  Thus,  interactions  are  limited  to  nearby  particles. 
Interactions  will  occur  through  reactive  zone  interactions  that  modify 
q",  lateral  energy  transport  in  the  condensed  phase,  and  fluctuations 
in  F (yl  from  the  mean.  All  aspects  are  related  because  they  depend 
upon  the  probability  of  finding  another  particle  with  specific  microstate 
within  a certain  distance  from  the  particle  being  considered  and  the 
variation  of  the  interaction  with  distance.  Since  phenomological  laws 
are  known  for  interactions  and  the  number  density  of  specific  particles 
can  be  computed  because  the  media  is  random,  these  interactions 
should  be  calculable.  * This  represents  the  current  frontier  of  the 
theory. 

Review  of  Current  Models 

To  define  the  current  technical  position  the  status  of  steady- 
state  combus tion  modeling  must  be  surveyed.  Current  combustion 
models  come  in  two  forms;  one-dimensional  and  statistical.  One- 
dimensional models  are  inherently  incapable  of  quantitative  description 
of  ballistic  phenomena  because  heterogeneity  is  not  included  in  an 
operational  fashion.  Since  one -dimens ional  combustion  models  have 
been  reviewed  elsewheref*^  * only  statistical  models  will  be 

reviewed  herein. 


Hermance  Model  - Steady-state  statistical  combustion 
modeling  was  initiated  by  Hermance.  His  creative  step  was  to  combine 
a detailed  physiochemical  model  for  spherical  oxidizer  particle  combus 
tion  with  statistical  concepts  defining  the  mean  microstate.  Figure  2 
illustrates  the  physiochemical  model  for  a fuel  surface  /oxidize  r parti- 
cle pair.  The  unique  feature  is  a crevice  at  the  oxidizer  /binder  inter- 
face wherein  heterogeneous  reaction  between  gaseous  oxidant  and  solid 
fuel  occurs . 

Hermance  employed  the  continuity  equation 
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where  the  mass  flaxes  on  the  RHS  were  means  for  burning  fuel  surface/ 
oxidizer  particle  pairs  and  the  areas  were  for  all  pairs  on  the  burning 
surface.  It  was  subsequently  assumed  that  the  oxidizer  regression 
rate  was  a step  function  , 
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where  t = o when  the  particle  first  penetrates  the  burning  surface  and 
tjgn  *s  t*le  igoifion  delay.  Since  the  heterogeneous  surface  reaction 
requires  oxidizer  decomposition,  mjr  = o for  - - * 


i t 


ign- 


Therefore,  for  o s t < t.  there  is  only  fuel  vapor  flow  from  the 
surface  of  a fuel  surface  ^oxidize  r particle  pair.  Consequently,  with 
an  ignition  delay  tjgn  * it  clear  that  the  burning  surface  is 
populated  with  some  fuel  surface /oxidizer  particle  pairs  possessing 
mox  = msr  = °*  Applying  Eq.  (10)  to  the  burning  surface  and  noting 
that  moX  = msr  = ° for  some  fuel  surface /oxidize r particle  pairs 
yields 


or 
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Contrasting  Eqs.  (1?)  and  (19)  shows  that  the  former  employs  mean 
fluxes  for  ignited  particles  with  surface  areas  for  all  particles.  This 
is  not  consistent.  Since  Sp< ox>  ign  < SPj  ox  and  Sp>  s r#  jgn  < Sp>sr, 
the  rate  predicted  by  Hermance's  equation  will  be  too  large  if  all  other 
things  are  equal.  As  the  inconsistancy  varies  with  tign  and  tjgn  varies 
with  pressure,  pressure  exponent  will  also  be  effected.  However, 
when  tign  ~ ° ,this  inconsistency  vanishes.  Therefore,  Hermance's 
theory  is  limited  to  small  ignition  delays  by  this  omission. 

Hermance's  analysis  assumes  the  fuel  surface  is  a plane. 
Therefore,  the  number  of  particles  intersecting  the  surface  Sp  per 
unit  of  Sp  is 

dto/dS^  » b dh) /dV  do) 


For  the  spherical  particles  assumed  in  this  analysis  the  volume  fraction 
of  oxidizer  is  ^ , and  the  planform  of  each  intersection  of 

a particle  with  Sp  is  circular.  Since  the  fraction  of  Sp  occupied  by 
oxidizer  is  ^ (dSi0,  /d%  = J)and  the  mean  diameter  of  intersections 
between  Sp  and  the  oxidizer  particles  is  defined  by  0,  / dS>^  ofTrBVq.') 
(dto/d'ipj  i algeb raic  manipulation  gives 
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This  dimension,  the  mean  for  ail  particle  intersections,  was  employed 
to  define  the  characteristic  dimension  for  the  mean  deflag rating  6tate. 
However,  since  all  particles  are  not  burning,  this  assumption  is 


inconsistent  with  physical  reality  unless  t; 


sion  defines  characteristic  latitudes  on  the  spherical  particle. 
Assuming  that  the  fuel  regression  rate  is  steady. 
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Note  that  this  dimen- 
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Therefore,  if  L is  the  vertical  distance  from  the  north  pole  of  an 
oxidizer  particle  tc  a characteristic  latitude,  the  time  required  for  the 
binder  to  regress  to  a characteristic  latitude  is 


With  the  step  function  assumed  for  r^  (t)  it  is  clear  that,  if  tjgn  < tc,N 
intersections  with  ignited  particles  can  occur  in  both  northern  and 
southern  characteristic  latitudes.  However,  if  ^ c 

only  the  southern  characteristic  latitude  is  acceptable.  Finally,  if 

> *c,s  , neither  characteristic  latitude  is  acceptable!  The  latter 
situation  means  the  model  breaks  down.  The  first  situation  implies  the 
characteristic  mean  state  is  degenerate  (northern  and  southern  charac- 
teristic latitudes  are  both  accessible).  Hermance  assumes  intersections 
occur  only  at  the  southern  characteristic  latitude.  This  requires  > 

N . However,  aforementioned  inconsistencies  require  . 

Hermance's  model  is,  therefore,  inconsistent  under  all  conditions. 
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Hermance  assumes  the  mean  fissure  width  £ is  given  by 

€ = b-b 

Therefore,  the  mean  heterogeneous  reaction  area  per  particle  is 


and  the  total  heterogeneous  reaction  area  ratio  is 


S*,.  / Sp  3 TTO  (b-B)  du/cKp  UAl 

The  fuel  surface  is  planar.  Since  +■  ^ * ip  and  Syi0«.  /Sp  * 

~ l"  J..  can) 

♦ This  can  be  "cured"  by  accounting  for  the  degeneracy  of  the  mean  state 
when  »lV,  < *eM 
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Thus,  all  area  ratio's  in  Eq . (17)  are  defined. 

By  assuming  the  temperature  of  oxidizer  and  binder  are  equal, 
oxidizer  decomposition  is  an  equilibrium  process  and  the  decomposi- 
tion product  is  an  ideal  gas,  and  kinetic  expressions  for  the  heterogene- 
ous reaction,  the  hete  rogeneous  reaction  mean  mats  fbix  was  related 
to  the  surface  temperature.  By  assuming  a pyrolysis  relation  for  the 
binder  its  mean  mass  flux  was  also  related  to  surface  temperature. 

For  steady-state  combustion  binder  and  oxidizer  mass  flows  must  be 
in  proportion  to  the  ingredients.  Therefore, 
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C onsequently,  the  only  unknown  is  the  surface  temperature. 

The  surface  temperature  was  defined  by  solving  a one-dimensional 
energy  equation  with  heat  release  at  the  burning  surface  and  at  a gas 
phase  lame.  The  latter's  standoff  distance  was  computed  by  assuming 
a second  order  reaction.  Consequently,  a complete  set  of  equations 
relating  mean  rate  to  pressure  and  initial  propellant  temperature 
through  roughly  seventeen  "parameters"  was  obtained. 

Hermance  also  considered  propellants  with  polydisperse 
spherical  oxidizer.  To  extend  the  unimodal  analysis  to  this  situation 
it  was  implicitly  assumed  that  fuel  surface/oxidizer  particle  pairs  with 
different  size  oxidizer  particles  possess  the  same  deflagration  rate.* 
Therefore,  all  have  the  same  surface  temperature  so  that  the  entire 
burning  surface  can  be  treated  with  the  energy  equation  employed  with 
the  monodisperse  propellant. 

Hermance  computed  rate  versus  pressure  for  two  monodisperse 
polysulfide / AP  propellants  and  compared  the  result  with  data.  figure 
7 shows  that  the  model  is  capable  of  virtually  duplicating  rale /pressure 
characteristics  for  a specific  propellant  but  that  extrapolation  (without 
parameter  changes)  to  the  same  formulation  with  different  particle  size 
is  definitely  not  quantitative.  Since  the  model  is  not  self-consistent, 
this  behavior  is  not  surprising. 

Another  major  difficulty  w'ith  Hermance's  model  is  the  dominant 
nature  of  the  crevice.  In  the  words  of  Reference  2 , "This  is  unfortu- 

nate because  in  actuality  the  oxidizer  does  not  regress  as  (Hermance) 
assumed;  it  does  not  maintain  an  overall  spherical  shape  and  there  is 
little  evidence  of  crevice  formation.".  Reference  2 notes  further  that 
SEM  pictures  tend  to  deny  crevice  existence  and  show  that  the  oxidizer 
particles  are  recessed  at  high  pressures  and  protrude  at  low  pressure. 
Therefore,  Ref.  2 concludes  that  "it  seems  very  unlikely  that  the 
combustion  mechanism  could  be  dominated  by  a nonexistent  cusp  (crevice) 
md  the  associated  reaction".  This  physical  evidence  led  Beckstead, 

Derr,  and  Priced*-*)  to  discard  Hermance's  unit  physiochemical  model 
as  physically  unrealistic  and  develop  a more  realistic  one. 

* This  will  be  termed  the  equal  rate  hypothesis. 
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Figure  3 Comparison  Theory  and  Experiment 
Hermance's  Model  (Fef.  231 
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Because  the  model  is  neither  physically  realistic  nor  self- 
consistent  there  seems  to  he  little  reason  to  pursue  this  model  further. 

Beckstead,  Derr,  Price  Model^"*  ''  - Beckstead,  Derr,  and 
Price  (BDP)  basically  embedded  a realistic  physiochemii  al  model  for 
fuel  surface  /oxidizer  particle  pair  combustion  in  the  statistical  proce-  F.usc 
dure  developed  by  Hermance  to  upgrade  that  model.  Figure  4 illustrate 
the  physical  model  for  a fuel  surface/oxidizer  particle  pair.  The  unique 
leatures  are  a tri-flame  structure  (oxidizer  decomposition  flame,  final 
diffusion  flame,  and  primary  flame)  and  geometric  structure  for  the 
deflagrating  oxidizer  particle.  The  flame  structure  is  a natural  modi- 
fication of  the  extended  GC  F flame  structure  developed  by  Summerfield 
and  co-workers^-*'.  The  latter  is  buttressed  by  considerable  experi- 
mental data.  The  permissible  oxidizer  surface  geometry  also  lies 
within  bounds  observed  experimentally.  Consequently,  there  are  very 
good  reasons  to  believe  that  the  phy Biochemical  model  is  realistic.  Put.  ^ WotEL 

Beckstead,  Derr,  and  Price  write  the  continuity  equation  as 

■-  ^ 1 ***  '•VV)  * 1%)  U1) 

where  is  the  total  surface  area  of  the  burning  surface.  It  is  impor- 
tant to  note  that  the  means  on  the  RHS  are  for  deflag  rating  fuel  surface/ 
oxidizer  particle  pairs  while  the  areas  are  for  all  pairs.  It  was  subse- 
quently assumed  that  oxidizer  regression  rate  was  a step  function 

= ro’<  < * > V') 

where  t = o when  the  particle  first  penetrates  the  burning  surface  and 
lign  's  t^le  ignition  delay.  Clearly,  since  an  ignition  delay  >0 

is  assumed,  the  burning  surface  will  contain  both  ignited  and  non-ignited 
fuel  surface/oxidizer  particle  pairs.  Therefore,  application  of  Eq.  (3) 
yields 

rp^  --  R*"  --  /%)  ♦[  C3l) 

Contrasting  Eqs.  (29)  and  (31)  two  differences  are  noted.  First,  if  a 
mean  mass  flux  for  ignited  oxidizer  is  to  be  employed,  the  associated 
area  must  also  be  limited  to  ignited  oxidizer.  Therefore,  the  BDP 
model  is  inconsistent  in  exactly  the  same  fashion  as  He rmance's  model. 

Consequently,  the  BDP  model  is  limited  to  situations  where  . 

However,  note  that  in  the  BDP  model  that  mean  rate  is  referred  to  the 
total  burning  surface  S0  rather  than  its  planar  projection  Sp.  Since 
strand  and  motor  burning  rates  are  based  on  a planar  surface  (because 
SQ  isn't  known),  the  burning  rates  computed  bv  the  BDP  model  cannot 
be  compared  directly  with  experimental  data.  The  latter  inhibits 
utilization  of  the  model. 

* Both  of  these  defects  are  correctable. 
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(23 ) 

Following  Hermance  the  dimension  D was  employed  to  define 
t and  hence  the  mean  microstate  of  a deflagrating  fuel  6urface/oxidizer 
particle  pair.*  However,  unlike  Hermance  the  duality  of  the  mean  micro- 
state  when  tjgn  z t was  recognized.  It  is  important  to  note  that  this 
recognition  is  clear 'in  neither  Ref.  (2)  nor  Ref.  (3);  however,  it  is  clear 
in  the  computer  cone.  Unfortunately,  D is  a characteristic  dimension 
f°r  all  particles  on  the  burning  surface  and  not  for  the  deflag  rating 
particles  to  which  the  dimension  is  applied.  This  inconsistency  also 
limits  the  model  to  situations  where  tjsn  ~ o. 
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To  define  the  oxidizer  surface  area  it  is  assumed  that  the 
deflagrating  surface  of  an  oxidizer  particle  is  a spherical  "cap", 
the  altitude  of  the  cap  at  the  mean  microstate  is 
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and  it  is  further  assumed  that  rf  = r,  h and  subsequently  Sox/S0  can  be 
expressed  in  terms  of  rJjX  and  tjgn. 


The  flame  structure  was  quantified  by  assuming  that  reactions 
were  conce  itrated  at  flame  standoff  distances  based  on  kinetic  and 
diffusional  considerations.  Heat  feedback  to  the  burning  surface  was 
based  on  an  area  ratio  weighted  summation  of  the  heat  feedbacks  from 
each  flame.  This  implicitly  assumes  each  flame  in  the  triad  acts 
independently.  Unfortunately,  this  is  not  the  case  for  it  is  clear  from 
Figure  4 that  the  final  diffusion  flame  cannot  communicate  directly  with 
the  burning  surface;  its  effect  is  derived  from  modification  of  the 
oxidizer  decomposition  flame.  That  is,  heat  released  in  the  final 
diffusion  flame  increases  the  temperature  at  the  oxidizer  decomposition 
flame  which,  since  it  is  kinetically  controlled,  increases  its  reaction 
rate.  This  causes  the  oxidizer  decomposition  flame  to  move  closer 
to  the  burning  surface.  It  is  this  "thermal  compression"  of  the  oxidizer 
decomposition  flame  by  the  final  diffusion  flame  that  causes  the  heat 
feedback  via  the  final  diffusion  + oxidizer  decomposition  flame  path  to 
exceed  that  from  an  independent  oxidizer  decomposition  flame.  It  must 
be  concluded  that  the  analytical  description  of  the  flame  structure  is 
inconsistent  with  the  physiochemical  model.  As  with  the  other  inconsis- 
tencies, these  can  also  be  "cured". 
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Beckstead,  Derr,  and  Price  have  employed  the  model  to  compute 
ambient  rate  pressure  characteristics  for  "unimodal"  polysulfide /AP 
propellants  (the  same  propeUantB  employed  by  Hermance).  Figure  5 
presents  a comparison  of  theory  and  experiment.  Since  the  20  /u. 
propellant  was  selected  as  a baseline,  the  theory/experiment  compari- 
son is  essentially  exact.  However,  the  theory  is  qualitative  at  best 
when  extrapolated  to  the  2 00 /u.  propellant.  Condon(25)  has  tested  the 
models  capability  for  predicting  temperature  sensitivity  of  the  JANNAF 
standard  propellant.  Figure  6 illustrates  the  comparison  between 

♦Another  dimension,  the  mean  width  of  the  fuel  surface  surrounding  the 
oxidizer  particle  is  involved.  This  width  was  computed  by  assuming  that 
the  propellant  has  an  ordered  body  centered  cubic  structure.  As  its 
structure  is  random,  this  dimension  is  also  slightly  in  error. 
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Figure  6h.  Comparison  Theory  and  Experiment  - BDP  Model 
JANNAF  Standard  Propellant  (Ref.  25) 


theory  and  experiment  for  both  cr  p and  rate/pressure.  Notice  that  while 
rate/pressure  behavior  is  in  reasonable  agreement  with  experiment  the 
theoretical  u p differs  substantially  from  experiment.  These  divergences 
between  theory  and  experiment  are  not  unexpected  in  view  of  the  afore- 
mentioned inconsistencies  in  the  model. 

Contrasting  the  theory/experiment  comparisons  illustrated  by 
Figs.  3 and  5 demonstrates  that  both  the  Hermance  and  BDP  models  fit 
the  same  data  base  equally  well  (remember  Hermance  adjusted  para- 
meters to  the  200  micron  data).  This  is  a rather  surprising  result 
because  Hermance's  unit  physiochemical  model  is  physically  unrealistic 
while  BDP's  isn't.  Why  did  this  occur?  The  only  sensible  explanation 
is  that  both  theories  have  sufficient  "adjustable"  parameters  available 
to  provide  an  acceptable  fit  to  virtually  any  limited  data  base.  Since 
very  few  of  these  parameters  are  defined  with  precision,  numerous 
parameters  are  available  for  curve  fitting.  This  result  demonstrates 
in  a very  concrete  fashion  that  physical  realism  is  not  guaranteed  by 
a model's  ability  to  predict/cor  relate  experimental  data.  In 
mathematicalese  the  ability  to  fit  experimental  data  is  a necessary  but 
not  sufficient  condition  for  a model's  physical  validity. 

Miller,  Hartman,  and  Myers  Model^J'  - Miller,  Hartman,  and 
Myers  (MHM)  were  also  influenced  by  Hermance's  innovative  blending 
of  combustion  model  and  statistics.  However,  rather  than  following  the 
path  of  a detailed  combustion  model  plus  statistics  for  monodisperse 
propellant  they  employed  an  existing  qualitative  combustion  model  and 
extended  it  to  polydisperse  situations.  This  approach  was  both  innova- 
tive and  practical.  It  was  innovative  because  it  embeds  particle  size 
explicitly  in  the  resulting  model  in  an  operational  fashion.  It  was 
practical  because  all  real  propellants  are  polydisperse.  In  addition, 
their  approach  is  computationally  undemanding. 

Miller,  Hartman,  and  Myers  began  their  analysis  for  steady- 
state  combustion  with  the  continuity  expression 

rwi  - r pe  bp  * «- 


In  addition,  application  of  continuity  to  fuel  and  oxidizer  flows  yields 

“ ft  r;  \ : h r;  %,f  ^ 

* Pm  “ Pm7„  cssl 


where  the  subscript  p denotes  a planar  projection  of  the  area.  For 
steady-state  deflagration  the  planar,  areal  mean  regression  rates  are 
identical 

ft,  • ft  * F 
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C onseq  uently. 


Recognizing  that  the  burning  surface  is  composed  of  a polvdispersion 
of  Particles,  that  the  whole  is  the  sum  of  its  parts,  and  that  in  Reynold's 
rules  of  averaging  the  mean  of  a sum  is  the  sum  of  the  means(  ),  they 
lorm  the  mean  as  a sum  over  all  particles  on  the  burning  suriace.  Thus, 
'■  ith  the  j index  representing  a particular  size  of  particle 
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Comparison  of  Eqs.  (38)  and  (10)  shows  that  MHM  invented  the  "petite 
ensemble"  method.  Since 


r 

‘hi 


t w) 


Kq.  (38)  can  be  rewritten  in  terms  of  individual  planar  regression  rates 

as 


2. 

i 


(mo) 


Assuming  that  the  burning  surface  is  planar  and  the  particles 
—s Pfleric al  the  number  of  j particles  per  unit  surface  is 


V%  - 6>)i  Air 6^) 


(Ml) 


Howeve 


r,  volume  fraction  and  weight  fraction  i*r  are  related  by 


ftp  ) 


Since  the  mean  diameter  of  the  intersected  particles  is 


and 


(M?) 


(M«) 


'See  discussion  prereding  Eq.  (21). 

*Note  that  Eq.  (10)  in  Ref.  (20)  is  in  error  by  the  factor  2/3.  As  every- 
thing is  later  lumped  into  a constant,  this  makes  no  real  difference. 


24 


MHM  showed  that 
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Eq.  (40)  becomes 
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It  appears  that  Eq.  (47)  applies  only  to  spherical  particles 
because  that  assumption  was  employed  in  the  derivation.  However, 
that  is  not  the  case  because  it  can  be  shown  that 


\w,,  /%  = W»'l 

* 

irregardless  of  the  shape  of  the  particle.  Thus,  the  major  restrictive 
assumption  at  this  point  is  a planar  burning  surface. 

To  complete  the  analysis  MHM  chose  the  Summerfield  GDF 
model(26)  m relate  rj  to  Dj  and  environmental  variables.  However, 
since  in  that  model 


r 
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(41) 


it  is  clear  that 
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Miller,  Hartman,  and  Myers  noted  that  this  limit  was  physically 
unattractive  and,  therefore,  modified  the  GDF  equation  to 


* Trans  late  a plane  Sp  a distance  A X in  the  propellant.  The  volume 


swept  out  is  S„  A 


?. 


definition 
given  by  S 


’!»%  AX 
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The  volume  of  j oxidizer  in  this  volume  is  by 
However,  this  volume  of  oxidizer  is  also 


AX 


Thus,  S 


V,OK 
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v-ht-re  r # func  (D  ,).  It  is  important  to  note  that  r is  not  dependent 
upon  its°neighbors . ^Therefore,  the  MHM  theory  neglects  interactions. 


To  achieve  their  final  burning  rate  expression  they  assumed 


that 
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S i ug  titution  of  Eqs . (51),  (52),  intoEq.  (47),  therefore,  yields 

r - a v * ^ » ky  v'v  1 / ^ 

With  this  rate  expression  they  noted  that 
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Thus,  they  concluded  there  were  definite  bounds  on  the  amount  of  rate 
control  achievable  with  particle  size  manipulation.  In  addition,  they 
concluded  that  both  AP  size  and  distribution  are  important.  Therefore, 
rate  cannot  be  correlated  with  a single  weight  mean  diameter. 


Equation  (53)  was  employed  to  correlate  strand  rate/pressure 
data  from  sixteen  different  formulations  (particle  size  distribution  was 
the  "variable").  Table  1 presents  a tabulation  of  measured  and 
calculated  values.  In  most  cases  agreement  is  very  good. 

This  work  and  Hermance's  must  be  rated  as  the  pioneering 
efforts  in  statistical  combustion  modeling.  The  defects  of  this  model 
are  those  noted(planar  burning  surface,  no  interactions),  those  of  the 
GDF  model,  and  a total  neglect  of  mixture  ratio  variations  among  the 
particles.  The  latter  may  be  significant  because  small  particles 
should  be  fuel  rich  relative  to  the  norm.  This  occurs  because  such  a 
small  amount  of  fuel  is  required  that  competition  with  neighbors  is 
slight.  On  the  other  hand,  a large  particle  is  probably  fuel  lean  rela- 
tive to  the  norm  because  it  requires  so  much  fuel  that  it  must  compete 
with  its  neighbors  to  such  an  extent  that  its  appetite  is  never  satiated. 
Since  composite  propellants  are  always  fuel  rich  and  as  a rough  rule 
rate  degrades  as  one  moves  away  from  the  stoichiometric  condition, 
it  can  be  expected  that  the  larger  particles  will  have  higher  rates  and 
the  smaller  particles  lower  rates  than  predicted  by  the  MHM  theory. 


verv  8Mvy  Poev. 


"oiicTT^n  has  noted  that  rQ  should  correspond  to  r^p  (P»  TB)  and  that  if 
it  does  Eq.  (51)  correlates  exponent  break  data  for  AP  composite 
propellants . 
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TABLE  1 


COMPARISON  THEORY  AND  EXPERIMENT  - MILLER'S  MODEL  - 
CTPB/AP/F  e^O  ^ POLYDISPERSE  PROPELLANT  (REF.  20) 


Prossun* 

(psia)  300  300  1300  2000 


Rurnln>'  R.Ur  (in/:.c<) 


Pi  ope  1 1 . uu 

exp  . 

calc  . 

exp  . 

calc  . 

exp  . 

. 

> X p . 

calc 

TP-  1 

0.3- 

0.32 

0.40 

0.39 

0.56 

0.57 

0.83 

0.63 

TP-2 

0.32 

0.34 

0.38 

0.42 

0.75 

0.64 

1 .50 

0.71 

TP-3 

0.39 

0.35 

0.39 

0.44 

0.69 

0.69 

0.88 

0.76 

TP-4 

0.39 

0.35 

0.43 

0.45 

0.76 

0.69 

0.87 

0.77 

TP-5 

0.36 

0.37 

0.45 

0.47 

0.73 

0.74 

0.84 

0.84 

TP-6 

0.41 

0.39 

0.45 

0.50 

0.84 

0.81 

0.99 

0.91 

TP-7 

0.41 

0.43 

-0.54 

0.57 

1 .06 

1 .05 

1 .25 

1 .23 

TP-6 

0.46 

0.43 

0.62 

0.61 

1.23 

1 .22 

1 .46 

1 .48 

TP-9 

0.45 

0.44 

0.62 

0.63 

1.27 

1 .28 

1.50 

1.55 

TP-10 

0.54 

0.50 

0.72 

0.70 

1.33 

1.41 

1.55 

1.69 

TP-11 

0.52 

0.49 

0.70 

0.69 

1 .34 

1.40 

1 .57 

1 .68 

TP-12 

0.40 

0.45 

0.60 

0.61 

1.23 

1.18 

1.49 

1.39 

TP-13 

0.56 

0.51 

0.76 

0.73 

1 .48 

1.51 

1.75 

1 .84 

TP- 14 

0.52 

0.51 

0.73 

0.71 

1 .44 

1.46 

1.65 

1 .77 

TP- 15 

0.58 

0.51 

0.76 

0.70 

1 .40 

1 .43 

1 .67 

1.72 

TP-16 

0.54 

0.51 

0.72 

0.71 

1.38 

1 .42 

1 .62 

1.72 

i 
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On  this  basis  it  would  be  expected  that  propellants  with  wide  distribu- 
tions would  be  correlated  worst.  Propellants  TP-7  to  -9  and  TP-11  to 
-16  have  the  widest  distributions  and  tend  to  bear  out  the  above  suspic- 
ions . 


(4  ) 

Cohen,  Derr,  and  Price  Model  - The  BDP  model  is  limited 
to  additive  free,  monodisperse  propellants.  However,  all  "real" 
propellants  are  polydisperse.  Therefore,  the  BDP  model  has  little 
relevance  to  practical  propellants  per  se.  Cohen,  Derr,  and  Price 
partially  "remedied"  this  situation  by  extending  the  model  to  bi-disperse 
propellants  with  an  inert  additive. 

In  the  extension  all  basic  assumptions  of  the  BDP  model  are 
retained.  The  major  new  assumption  is  that  all  sizes  of  oxidizer 
particles  have  the  same  burning  rate.*  Therefore, 


rrr 


wht  re  k,  j are  in  the  same  set  but  k # j.  With  this  assumption  and  the 
i'.DP  model's  ignition  delay  criteria  the  spherical  cap  height  hj  can  be 
i omputed  for  each  oxidizer  particle  size  [see  Eq.  (32)]  . With  the  cap 
height  known  and 

- V i/i 1 Ci , 
i 4 

the  geometry  of  the  deflagrating  oxidizer  surface  is  known  and, 
therefore,  (S0x'^p,ox)j  computed  for  each  particle  size.  Since 
Sp  ox  . / Sp  = "j  , the  total  oxidizer  surface  area  is 

S../S,  - \ 

for  the  BDP  mode’  the  ratio  Sox/S^  is  desired.  As 

V -A  "S../V- 

is  known,  SQ/Sp  and,  subsequently,  SoX/Se  can  be  computed.  Therefore 
only  mJJx  needs  to  be  computed  to  define  r with  Eq.  (29). 

To  compute  moX  the  energy  balance  employed  by  BDP  with 
modifications  to  account  for  the  presence  of  inert  additives  is 
employed.  To  compute  heat  feedback  from  the  gas  phase  reaction 
zone  the  BDP  methodology  is  employed.  The  only  modification  being 
that  the  characteristic  fuel  dimension  is  modified.  It  is  computed  on 
the  assumption  that  all  particles  have  the  same  annular  thickness  of 
fuel  surrounding  them. 


'--Termed  the  equal  rate  hypothesis  herein. 


Computations  derived  from  the  model  were  presented.  Figure 
7 presents  a comparison  of  theory  and  experiment  for  the  ambient  rate/ 
pressure  characteristics  of  a specific  propellant.  Note  the  excellent 
ag  reement. 

This  extension  builds  upon  the  BDP  model  without  correcting 
any  of  its  inconsistencies.  Consequently,  the  "sins  of  the  father  are 
visited  on  the  son".  In  addition,  it  is  assumed  that  all  particles  possess 
the  same  rate  irregardless  of  size.  This  assumption  seems  highly 
unlikely  in  view  of  the  well  known  fact  that  particle  size  is  a major 
lactor  in  the  rate  control  of  AP  composite  propellants.  The  fact  that 
agreement  with  experiment  for  a specific  formulation  was  achieved 
should  not  be  given  too  much  weight.  Recall  that  Hermance's  invalid 
model  and  MHM's  model,  which  assumed  m^1  t mj^  , both  correlate 
data. 


Note  also  that  although  the  Cohen,  Derr,  Price  model  has  incon- 
sistencies and  is  limited  to  two  particle  sizes  that  these  inconsistencies 
can  be  removed  and  the  model  extended  to  a true  polydispersion  by 
following  the  path  indicated  above. 


Sammons  Model 


(5) 


The  Sammons  model  extends  the  BDP 


model  to  propellants  with  polydisperse  oxidizer  and  inert  additives. 

The  basic  path  followed  is  that  blazed  by  Cohen,  Derr,  and  Price. 
However,  several  unique  features  were  added.  First,  the  oxidizer  is 
segregated  into  two  classes:  subcritical  and  supercritical.  The  former 
is  assumed  to  undergo  condensed  phase  reaction  "at  the  burning  surface" 
thereby  modifying  the  surface  heat  release  term  heretofore  assumed 
constant  by  both  BDP  and  CDP.  The  supercritical  oxidizer  is  assumed 
to  burn  as  described  by  BDP.  Second,  the  ignition  delay  is  handled 
differently  than  by  Cohen,  Derr,  and  Price.  Third,  improvements  in 
the  computation  of  thermophysical  properties  and  the  short  flame  Burke - 
Shuman  solution  were  introduced. 


The  surface  geometry  of  the  supercritical  oxidizer  is  handled 
differently  than  by  Cohen,  Derr,  and  Price.*  Instead  of  computing  hj, 
Dj,  and  (Sox/Sogjp)j  for  each  particle  size  and  then  summing  an  average 
ignition  delay  is  computed  for  the  supercritical  oxidizer  as 


£ wr  * 
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where  tjgn  j = K0  D-  /p  is  the  expression  employed  by  BDP.  Then 
the  BDP  expression  for  h/D  is  employed,  i.  e. 


t/b  --  0.5  ( l t i/VTK'-  r.«/?  )+  r„,  >• 


(UOj 


*The  description  of  what  was  actually  done  is  very  unclear. 


29 


With  this  h/D  the  value  of  S0x^b  is  computed  with  the  BDP  expression. 


(1.0 


where  the  + and  ■ superscripts  denote  the  i terms  in  Eq . (60). 

Rate  is  basically  obtained  from  the  energy  balance  employed 
by  BDP  but  modified  to  account  for  inerts  and  6ubcritical  AP.  The 
inerts  are  treated  according  to  Cohen,  Derr,  and  Price.  However, 
details  of  the  subcritical  oxidizer  treatment  are  unclear.  The  criteria 
specifying  critical  oxidizer  particle  size  and  the  specific  surface 
energy  release  associated  with  the  subcritical  oxidizer  are  not 
discussed. 


Rate  pressure  characteristics  were  computed  for  several 
polydisperse  propellants  and  compared  with  data.  Results  show  that 
qualitative  agreement  exists. 

This  model  suffers  from  the  inconsistencies  inherent  in  the 
BDP  model.  Moreover,  it  incorporates  the  equal  rate  hypothesis. 
Additionally,  the  method  employed  to  obtain  the  average  Sox/S ^ has 
no  sound  physical  basis.  Finally,  the  description  of  this  model  is 
very  poor.  To  find  out  what  is  going  on  one  must  literally  decode  the 
computer  program  (an  arduous  task).  In  this  writer's  opinion  this  is 
not  worthwhile. 

Click's  "Grand  Ensemble"  Model^  ’ - Click  noted  that  the 

burning  surface  of  heterogeneous  propellants  was  an  ensemble  of 
different  microstates  and  criticized  the  BDP  theory  for  (a)  collapsing 
the  microstates  into  a single  mean  microstate  and  (b)  selecting  the 
dimension  characterizing  the  mean  microstate  from  purely  geometric 
arguments.  4 To  circumvent  this  averaging  process  a statistical 
procedure  was  devised  to  satisfy  continuity  and  thereby  compute  the 
mean  rate  in  terms  of  microstate  rates.  In  addition,  Glick  noted  that, 
if  the  surface  statistics  were  invariant  in  time,  the  mean  nonsteady 
properties  could  be  expressed  in  terms  of  nonsteady  properties  for 
single  flames. 

Two  di stribution  functions  Fox  and  Ff  were  introduced  such  that 
the  fraction  of  fuel  surface /oxidizer  pairs  with  + and 

4 £„»  5 £0,  + d€ol  is 

cbJ  - Wu  P,  ci£oi  cU*  (U) 


Assuming  that  the  mass  flux  of  oxidizer  was  given  by  the  functional 


*As  shown,  this  criticism  was  well  founded. 
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mn  ' /«;,  U.T.,  €„,  €4,  - ) H.*i 

the  mass  flow  of  oxidizer  from  pairs  with  is 

d w o*  (du/d%)  F,  d€ot  cU;  CM) 


rherefore,  the  total  oxidizer  flow  is  given  by  integrating  over  all  possi- 
ble and  . For  a monodisperse  situation  04 

For  there  is  no  reason  for  an  upper  bound;  however,  the  fact  that 

the  particles  cannot  interpenetrate  means  that  the  smallest  bit  of  fuel 
available  is  always  non-zero.  Therefore,  s 

Consequently, 
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If  mixture  ratio  is  preserved  (true  in  steady-state) 


It 

/Y*  ' 
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If  changes  in  time  are  quasi-steady,  mixture  ratio  is  invariant 
in  time.  Therefore,  if  the  surface  statistics  are  also  invariant  in  time, 
Fox  and  Ff  are  time  independent.  Consequently. 


Since  the  mean  pressure  coupled  response  function  is  defined  as 

<■  i _ « \ / r i i \ cu  6) 


aCd«*>  / ««  )/  (d^/^j 

manipulation  of  Eqs.  (67)  and  (68)  gives 
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It  was  noted  that  a similar  expression  held  for  the  velocity  coupled 
response  function  Ry  and  that  one -dimensional  functional  forms  for 
Rp  (and  Ry)  could  be  employed. 

To  evaluate  the  integrals,  Fox  ( fc„.  ),  Fj.  ( ),  and  m" 

need  to  be  evaluated.  To  evaluate  Fox  and  Ff  it  was  assumed  that  all 
accessible  and  were  equally  probable  and  constraints  on 

number  and  surface  area  applied.  It  was  demonstrated  that  under  these 
conditions  the  distribution  functions  possessed  the  functional  form  of 
the  Boltzman  factor.  Therefore,  FQX  and  Ff  were  computed  for  propel- 
lants with  monodispe  rse , spherical  oxidizer.  The  distribution  functions 
were  extended  to  propellants  with  polydisperse,  spherical  oxidizer. 

The  FQX  extension  was  rigorous;  the  Ff  extension  was  approximate. 

To  demonstrate  the  method  the  GCF  model^  ^ was  employed 
for  m"  and  the  rate/pressure  characteristics  of  the  resulting  model 
for  monodisperse  propellants  was  tested  against  data^"* The  com- 
parison showed  that  the  "statistical"  CD F and  CDF  models  gave 
almost  identical  results. 

This  work  represents  an  interesting  line  of  attack.  However, 
no  attention  was  given  to  energy  conservation  and  the  fact  that  the 
microstates  were  themselves  nonsteady.  These  considerations  invali- 
date utilization  of  the  GDF  model  in  this  statistical  framework.  Indeed, 
no  existing  model  can  be  employed  in  this  framework.  Finally,  the 
distribution  function  F0x  is  incorrect  because  all  €e,  are  not  equally 
probable;  all  depths  of  intersection  of  the  particle  with  Sp  are 

As  noted  in  Fundamentals  of  Statistical  Combustion  Modeling 
the  worst  defect  of  this  approach  is  the  computational  problem  associa- 
ted with  tracking  individual  microstates  through  time.  Because  of 
this  fact  alone,  this  approach  possesses  little  current  interest. 

f z -ft) 

Cohen's  Nitramine  Model'  - Virtually  all  of  the  statistical 

combustion  models  apply  to  AP  propellants.  However,  Cohen  has 
adapted  the  BDP  model  to  propellants  with  nitramine  oxidizers.  This 
development  has  (and  still  is)  proceeded  in  several  stages.  In  the 
first  stage,  the  BDP  model  was  modified  to  treat  monodisperse  nitra- 
mine propellants.  In  the  second  stage,  the  Cohen,  Derr,  Price  metho- 
dology was  employed  to  extend  the  monodisperse  model  to  a bidisperse 
situation  with  mixed  oxidizers.  In  the  third  stage,  active  binders  were 
explored. 

The  initial  nitramine  propellant  model  developed  by  C ohen^ 
was  limited  to  monodisperse  nitramine  propellants  and  basically  con- 
sisted of  a nitramine  unit  physiochemical  combustion  model  embedded 
in  the  BDP  statistical  framework.  The  unit  physiochemical  model  was 


*Suppose  the  particle  is  a rod  with  pointed  ends.  Obviously,  the  pro- 
bability of  =.0  is  now  much  less  than  that  for  G„  - tTt>  M 
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based  on  ballistic  cinephotographic  and  SEM  (extinguished  by  rapid 
depressurization)  investigations  of  nominal  195  micron  and  5 micron 
additive  free  monodisperse  HTPB/HMX  propellants  at  the  75%  total 
solids  level  over  a wide  pressure  range.  Consequently,  the  model 
incorporates  the  following  features;  no  energy  release  in  the  6ub-melt 
layer  (supported  by  SEM  results),  existence  of  oxidizer  and  binder 
melts  (supported  by  cinephotographic  evidence),  provision  for  termi- 
nation of  oxidizer  melt  (supported  by  cinephotographic  evidence),  and 
provision  for  deep  "flame"  penetration  into  the  propellant  via  "deep 
penetration"  HMX  crystal  deflagration  (supported  by  cinephotographic 
evidence).  Figure  8 illustrates  the  model  for  three  situations: 

Nitramine  melting,  nitramine  not  melting,  and  nitramine  not  melting 
with  deep  penetration.  When  the  nitramine  melts  it  is  assumed  that 
the  nitramine's  burning  surface  is  planar.  However,  when  there  is  a 
binder  melt,  that  melt  may  encroach  upon  the  nitramine  surface  as 
shown.  This  possibility  i6  embedded  by  a two  parameter,  purely 
empirical  expression.  When  the  nitramine  does  not  melt,  the  nitra- 
mine surface  geometry  is  determined  via  the  same  geometric  rela- 
tions (an  ignition  delay  is  utilized)  employed  by  BDP  with  the  exception 
that  the  recessed  state  is  unbounded.  That  is,  recession  is  allowed  to 
extend  beyond  the  boundaries  of  the  deflagrating  HMX  particle  to 
account  for  deep  penetrations.  In  this  extension  boundedness  is  intro- 
duced through  an  empirical  manufacturing  parameter. 

A question  of  prime  importance  to  the  model  is  "under  what 
conditions  does  HMX  melting  cease?"  This  is  answered  on  a thermal 
basis.  Basically,  an  HMX  particle  is  dropped  into  binder  whose 
temperature  time  history  simulates  the  thermal  wave  and  removed 
after  its  transit  time  (thickness  wave/burn  rate).  Obviously,  if  transit 
time  is  very  short  melting  will  not  occur.  Consequently,  this  proce- 
dure defines  a critical  burning  rate;  above  it  no  melting;  below  it 
melting.  By  identifying  exponent  break  point  with  HMX  melt  termination, 
Cohen  shows  that  this  procedure  correlates  break  point/particle  size 
data  in  the  40  - 300  micron  6ize  range;  data  in  the  < 40  micron  range 
is  not  correlated  (see  Figure  9 ). 

With  particle  geometry  defined,  the  energy  equation  basically 
closes  the  problem  by  defining  the  mean  surface  temperature.  The 
treatment  here  is  also  virtually  identical  to  that  employed  by  BDP  with 
the  exception  that  a:  endothermic  term  is  included  to  account  for  nitra- 
mine melting  (when  melting  occurs). 

Execution  of  the  model  shows  very  good  agreement  with  the 
rate/pressure  characteristics  of  the  experimental  propellants. 

Figure  10  illustrates  typical  results.  Table  2 presents  a 
list  of  the  standard  parameters  employed.  As  noted  previously,  these 
parameter  values  are  not  rigorously  defined.*  The  exponent  break 

* The  intention  here  is  not  to  "put  down"  this  analysis  because  all  detailed 
models  must  utilize  roughly  the  same  amounts  of  "parameters".  However, 
to  disregard  the  ability  of  parameter  adjustment  to  improve  data  correla- 
tion would  be  to  lose  contact  with  reality. 
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figure  8.  Cohen's  Nitramine  Surface  Geometry  (Ref.  6) 
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TABLE  2 

I'AUAMLTKRS  EMPLOYED  IN  COHEN  N'lTKA.VlINE  MODEL  (REF.  (,) 


a 

Weight  percent  oxidizer 

75.0 

D 

o 

Mean  diameter  of  oxidizer 

195(i,  1:(j 

P 

Propellant  density 

l.r'  . c c 

p( 

Binder  density 

0.03  g/cc 

P 

Oxidizer  density 

1 .9  1 g/cc 

ox 

Qf 

Binder  heat  of  decomp  -sitmn 

56 9 cal/g 

T 

K 

] : ‘pellant  flame  temperature 

167SfK 

M 

Gas  molecular  veight'1' 

17.0  g m 1 

V 

Diffusion  parameter1'1 

10*  cm*  - atm/ sec  - °K 

'!» 

£>?>.. ich:  <metric  ratio'1' 

53 

C 

D.ffurii  n parameter11  ’ 

17 

Gign 

Oxidizer  igni tion  del  iy  constant 

190  sec  atm°‘72/cmOJ 

Pign 

Oxidizer  ignition  delay  pressure -dependence 

0.72 

D » g n 

Oxidizer  ignition  delay  diameter-dependence 

0.8 

Oxidizer  heat  of  decomposition 

-225  cal/g 

qlm 

Oxidizer  heat  of  fusion 

132  cal,  g 

Ef 

Activation  energy  of  binder  decomposition 

16.9  Kcal/mol 

Activation  energy  of  oxidizer  decomposition 

50  Kcal/mol 

Af 

Pre factor  for  binder  decomposition 

2 99  g/c.-r'-scc 

A 

OX 

Prefactor  for  < xidizer  decomposition 

5 x I09g/c  n/  - sec 

6 

Order  of  diffusion  flame  reaction 

2 

6 

Order  of  oxidizer  flame  reaction 

2 

kPF 

Prefactor  for  diffusion  flame  reaction 

30  g/cu.1- sec -atm2 

k 

Prefactor  for  oxidizer  flame  reaction 

0. 2 4b  . -sec  - atm2 

OX 

T 

ox 

Oxidizer  flame  temperature 

3275°K 

X 

Gas  conductivity 

O.OOCJ  cal/cm  - sec -°K 

\ 

o 

Oxidizer  conductivity 

0.0i»049  cal/cm-sec-°8 

9 

Binder  conductivity 

0.00044  cal/cm- sec-°8 

c 

p 

Heat  capacity 

0.3  cal/g -°K 

K 

0 

Oxidizer  diffusivity 

0.001 1 cm2/  sec 

Kf 

Binder  diffusivity 

C.001  1 cm2/sec 

T 

m 

Oxidizer  melting  point 

2 78"C 

K 

Coefficient  for  melt  surface  structure- 
diameter  dependence 

1 

a 

Melt  surface  structure  diameter 
dependence 

0.25 

N 

Manufacturing  parameter 

3,  1 

n 

Diffusion  flame  pressure-dependence 

0.5 

T 

o 

Initial  propellant  temperature 

2 5°C 

P 

Pressure  range 

10  - 2200  atm 
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associated  with  the  nominal  195  micron  propellant  is  associated  with 
the  termination  of  HMX  melting  and  arises  from  the  sudden  elimination 
of  the  HMX  heat  of  fusion  from  the  energy  equation  (increases  rate) 
and  the  removal  of  the  melting  geometric  constraints.  The  latter  lets 
SQX  suddenly  increase  thereby  suddenly  increasing  rate. 

In  Ref.  7,  Cohen  extended  the  monodisperse  model  to  handle  a 
bidisperse  mixture  of  two  nitramine  ingredients.  That  is,  two  particle 
sizes  of  nitramine  A and  two  particle  sizes  if  nitramine  B were  consi- 
dered. This  extension  represents  a melding  of  Cohen,  Derr,  and 
Price's  extension  of  the  BDP  model  from  the  monodisperse  case  to  the 
bidisperse  case",  Cohen's  unit  nitramine  phy6 iochemical  model^  5), 
and  a "relaxation"  of  the  "equal  burning  rate  hypothesis".  The  latter 
was  necessitated  by  the  mixed  oxidizer  situation.  In  essence  Cohen, 

Derr,  and  Price's  bidisperse  procedure  (and  its  concomittant  equal 
rate  hypothesis)  was  applied  separately  to  each  nitramine  ingredient 
with  Cohen's  nitramine  unit  physiochemical  model  replacing  the  BDP 
unit  physiochemical  model  and  a separate  energy  equation  for  each 
nitramine  specie.  Mean  rate  was  computed  as  a weighted  sum  of 
specie  rates.  Thus,  this  development  represents  a hybrid  of  "equal 
rate  hypothesis"  and  unequal  rate  statistics.  Figure  11  illustrates 
typical  results. 

In  Ref.  8 Cohen  extended  the  preceding  model  to  include  active 
binders  and  bidisperse  two  specie  mixtures  of  nitramine /nitramine  and 
nitramine  / AP  oxidizers.  The  inactive  binder  situation  for  bidisperse 
nitramine /A P oxidizers  was  a logical  extension  of  the  nitramine/ 
nitramine  situation  and  was  achieved  by  simply  embedding  BDP's  unit 
AP  physiochemical  model  in  the  nitramine /nitramine  framework.  In 
other  words,  instead  of  computing  the  rates  of  bidisperse  nitramine  A 
and  bidisperse  nitramine  B and  then  employing  the  weighted  sum  as  the 
mean  burning  rate,  bidisperse  nitramine  B was  replaced  by  bidisperse 
AP Three  schemes  were  advanced  for  including  an  active  binder. 

"First,  it  was  assumed  that  the  regression  rate  of  the  binder  at  any 
pressure  was  its  intrinsic  value  (no  oxidizer  present).  Second,  the 
constraint  of  oxidizer -binder  continuity  was  relaxed.  Third,  in  the 
case  of  nitramine  oxidizers  it  was  assumed  that  there  was  no  diffusion 
flame  because  the  nitramine  and  binder  are  stoichiometrically  balanced.  " 

No  equations,  etc.  were  produced  to  detail  active  binder  inclusion  in  the 
model.  The  results  of  a sequence  of  calculations  was  presented  for 
nitramine /nitramine  with  active  and  with  inactive  binders  and  mixed 
nitramine/AP  with  active  and  with  inactive  binders.  Cohen  noted  that 
these  results  follow  experimental  trends. 

Cohen's  work  to  date  represents  a powerful,  pioneering  thrust 
toward  realistic  modeling  of  nitramine  propellants.  The  only  elements 
omitted  are  true  polydisperse  capabilities*  and  the  ability  to  treat 

*The  question  here  is  that  of  representing  a bimodal  mixture  of  polydisperse 
lots  by  two  mean  particle  sizes.  The  results  of  Miller,  Hartman,  and 
Myers(^)  suggest  that  one  cannot.  However,  this  was  for  AP  composites. 
In  any  event,  Cohen's  model  cannot  handle  true  polydisperse  blends. 
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catalysts  explicitly.  However,  the  model,  since  it  employs  BDP  metho- 
dology, possesses  the  noted  inconsistencies  associated  with  that  model. 
Moreover,  with  deep  penetrations  the  fuel  surface  is  no  longer  planar. 
Consequently,  the  statistics  of  the  surface  are  altered.  These  inter- 
actions were  not  treated  by  Cohen. 


Miller.  Donohue,  and  Peterson  Model  " In  this  work  the  model 
ol  Miller,  Hartman,  and  Myers  is  extended  to  include  interaction  effects. 
This  extension  was  stimulated  by  the  finding  that  the  MHM  theory  did  not 
correlate  particle  size  effects  well  in  uncatalyzed  HTPB/AP  propellants. 
In  particular,  it  was  found  that  emphasis  was  weighted  to  the  coarse  end 
of  the  distribution  rather  than  the  fine  end  of  the  distribution.  As  noted 
in  the  review  of  the  MHM  model,  consideration  of  mixture  ratio  effects 
will  tend  to  skew  the  weighting  in  the  direction  required. 


Interactions  among  particles  were  assumed  to  occur  from  flame 
interactions.  Specifically,  it  was  assumed  that  the  diffusion  flame  about 
an  oxidizer  particle  was  larger  than  the  particle  so  that 
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where  i . It  was  further  assumed  that  all  oxidizer  particles  lying  in 
the  jth  particles  annular  area 


Lk 


M 


- t.  ) /4 


Cftlj 


C.OOK1T  \ 


s{c 

would  be  "excluded".  If  N"  is  the  number  of  j particles  per  unit  planar 
surface,  the  total  "exclusion  area"  of  particles  is 


* 17  C^-0  bj  /4  (81) 

an  i*-^1  particle  to  be  excluded  by  a j^1  particle  it  must  fit  in  the 
1th  particle's  annular  "exclusion  area".  Thus, 


= \ is*) 

txcv_ob€t» 

With  Dj  increasing  with  j this  defines  a critical  j = jc,i  such  that  only 
jth  particles  with  j jc,i  wiH  exclude  an  i^1  particle.  Since  the 
number  of  ith  particles  per  unit  area  is  N"  , the  number  of  ith  particles 
excluded  by  _a_  j ^ jc  j particle  is 

Ab",«  = iru"  wj’cV-O  (*s) 

*This  means  that  the  contribution  of  the  excluded  particles  to  the  total  flow 
of  products  is  deleted.  This  is  patent  nonsense  because  mass  is  conserved. 
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The  total  number  of  i1'’1  particles  "excluded"  is  the  sum  over  all  j1*1 
particles  with  j *•  jc  j.  Thus, 

' (_’,T  <L  ki  b C6U^ 

A i i 

Consequently,  the  number  of  it*’1  particles  that  "contribute"  is 

--  kjV  --  b-trU-O/*]  Z M-  h l ceil 

It  is  clear  that  this  postulated  strategy  "skews"  effects  to  the  larger 
particles  because  large  particles  "exclude"  small  particles  and  not 
vice  ve  rsa. 

As  noted  previously  in  the  review  of  Ref.  (9)  lo"  oC  W-  b* 

Thus,  employing  the  rate  expression  of  Ref.  (9) 

F=  I F ur  I i-  X.UM  Z bit  / f «r  fW) 

L h.i  - *■ 

where  K is  a constant. 

Miller,  Donohue,  and  Peterson  note  that  this  is  only  a partial 
expression  because  particles  hiding  particles  can  be  hidden.  This  leads 
to  the  expression 

p T f “1  Fi  ' <,l  ? «i  ^ I''*  (M) 

where  is  an  empirical  constant.  They  have  successfully  emploved 

this  expression  to  correlate  ambient  rate/pressure /particle  size 
distribution  data  for  HTPB/AP  propellants. 

This  approach  represents  a pioneering  approach  at  pair/pair 
interactions  and  demonstrates  how  the  statistics  of  the  burning  surface 
can  lie  employed  to  deduce  interactions  within  the  framework  of  an 
assumed  interaction  model.  For  these  concepts  this  work  is  valuable. 
However,  "excluding"  particles  is  nonsense  because  they  are  still  there 
and  their  mass  must  go  into  gaseous  products.  In  order  words,  the 
postulated  interaction  model  is  not  physically  plausible.  As  mentioned 
previously,  there  are  sound  physical  reasons  to  assume  that  small 
particles  are  more  fuel  rich  than  large  particles.  Since  the  propellant 
as  a whole  is  underoxidized  and  rate  usually  falls  off  as  one  moves  away 
Irom  the  stoichiometric  condition,  rate  will  be  skewed  away  from  that 
predicted  by  the  GDF  model  by  mixture  ratio  effects.  As  Ebenezer, 
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Cole,  and  McAlevy  have  formulated  a mixture  ratio  sensitive  GDF 
model,  replacement  of  the  GDF  model  in  MHM's  theory  by  this  model 
would  make  more  sense  than  physically  implausible  interaction  mechan- 


DEVELOPMENT  OF  POL YDISPERSE  MODEL 


General  Comments  - The  work  to  be  reported  in  this  section  was 
conducted  over  a two-year  span.  The  goal  of  the  first  years  work  was 
primarily  to  embed  the  BDP  model  in  a "correct  statistical  formulation" 
thereby  effectively  extending  the  model  to  additive  free  propellants  with 
spherical  oxidizer.  The  goal  of  the  second  years  effort  was  basically 
to  extend  this  model  to  aspherical  particles  and  include  binder  melts. 

The  first  goal  was  accomplished  by  abandoning  the  original  grand 
ensemble  statistics  and  developing  the  petite  ensemble  method.  However, 
in  the  operation  of  this  code  "problems"  appeared.*  In  the  process  of 
exploring  these  "problems"  further  insites  into  the  petite  ensemble 
method  were  glimpsed  and  new  methodology  for  selecting  the  mean 
state  was  developed.  This  impacted  the  aspherical  particle  work  tough- 
ening that  problem.  During  investigations  of  spherical  problems  (late 
in  the  program  when  petite  ensemble  methodology  was  thought  to  be 
surprise  free)  a further  "glimmer"  into  the  potential  of  the  method 
was  obtained.  This  "glimmer"  showed  that  retention/explusion  of 
partially  consumed  oxidizer  particles  and  "deep  penetrations"  (see 
review  of  Cohen's  Nitramine  Model)  should  flow  naturally  from  a 
proper  combination  of  petite  ensemble  statistics  and  BDP  model. 

These  "enlightenments"  have  been  stimulating  but  troublesome;  they 
interfere  with  orderly  progress  and  reporting.  As  for  the  reporting, 
the  method  is  developed  with  generality  as  far  as  possible.  After  this 
point,  main  emphasis  is  placed  on  spherical  particles . However,  all 
partially  completed  work  is  included  in  appropriate  appendices  for 
future  reference. 


In  retrospect  there  should  have  been  no  "surprises".  However, 
it  just  wasn't  that  way.  The  spirit  of  the  first  "enlightenment"  is  f 

included  in  the  review  of  the  BDP  model.  The  second  "enlightenment" 
consisted  of  realizing  that  under  most  conditions  an  oxidizer  particle 
is  either  consumed  before  the  fuel  plane  reaches  its  south  pole  or  not.  For  11  ' ' 
latter  situation,  which  occurs  with  large  particles  at  low  oressures,  the  Re: 
partially  consumed  oxidizer  particles  are  either  "freed"  or  retained.  fV' 

The  retention  question  is  probably  answered  by  response  to  the  question 
"is  there  a surface  melt"?  If  there  is,  the  partially  consumed  parti- 
cles probably  stick  to  the  surface.  If  not,  they  are  probably  freed.  g 

If  they  stick,  there  are  more  oxidizer  particles  on  the  burning  surface 
than  computed  by  planar  statistics  because  planar  statistics  assumes 
particles  vanish  when  the  fuel  plane  reaches  a particles  south  pole. 

The  stuck  particles  are  different  than  the  original  particles;  they  are 
smaller.  Thus,  a monodispersc  propellant  can  become  effectively 
bi-disperse  under  certain  conditions.  On  the  other  hand,  if  partially 
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*See  Ref.  (28)  for  some  details. 
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consumed  solid  particles  are  freed,  oxidizer  flows  from  the  burning 
surface  in  both  gaseous  and  solid  (freed  particles)  forms.  As  the 
continuity  expression  employed  counts  only  the  gaseous  form,  that 
expression  needs  modification. 
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If  the  oxidizer  particle  is  consumed  before  the  fuel  plane 
reaches  its  south  pole,  a fuel  rich  depression  is  left.  Therefore, 
some  of  the  microstates  counted  by  the  statistics  are  devoid  of 
oxidizer  ! However,  as  the  statistics  are  steady  in  time,  the  ergodic 
surmise  applies  so  that  some  of  these  particle  absent  microstates 
persist  for  finite  time.  If  there  is  no  surface  melt,  continued  pyro- 
lysis of  the  binder  at  the  particle  absent  sites  can  expose  other  parti- 
cles which,  if  there  is  time,  can  ignite  and  burn.  It  is  readily 
apparent  that,  if  rQX  is  large  and  tjgn  very  small,  a multi-particle 
(deep)  penetration  is  possible.  On  the  other  hand,  suppose  a mobile 
melt  is  present.  The  possibility  now  exists  for  melt  to  fill  the  empty 
site  making  it  impossible  for  another  particle  to  be  uncovered. 
Response  to  the  question  "are  particles  uncovered  at  empty  sites?" 
depends  upon  more  than  must  the  presence  or  absence  of  a surface 
melt.  It  must  be  a dynamic  process  because  melt  motion  is  required. 
It  is  important  to  note  that  if  the  empty  sites  do  not  fill,  the  burning 
surface  is  no  longer  quasi-planar.  This  is  especially  true  if  other 
particles  are  uncovered  and  ignited.  Consequently,  a rough  surface 
must  be  allowed  for. 
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These  questions  are  not  simple;  however,  they  are  also  not 
unanswerable.  Indeed,  the  beauty  of  the  petite  ensemble  method  is  that 
answers  to  these  questions  appear  possible  within  its  framework  (a 
little  bending  may  be  necessary). 

Statistical  Framework  - In  Fundamentals  of  Statistical  Combus- 
tion Modeling  it  has  been  shown  that 
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where  m£J,  j is  the  mean  mass  flux  from  and  Sp^  j the  portion  of 
Sp  occupied  by  the  jth  distinguishable  particles.  Particles  on  the 
burning  surface  are  distinguishable  because  of  specie,  size,  shape, 
and  orientation.  If  the  particles  are  spherical,  their  shape  is  defined 
and  there  are  no  distinguishable  orientations.  Therefore,  only  specie 
and  size  represent  distinguishable  features.  However,  for  a general 
particle  three  parameters  must  be  specified  to  define  shape  and  size 
and  three  angles  must  be  specified  to  define  orientation.  Consequently, 
a total  of  seven  parameters  is  required  to  define  a general  distinguish- 
able particle.  Since 
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where  K is  the  total  number  of  parameters  and  Nj<  is  the  number  of  steps 
in  the  k1*1  parameter  required  to  adequately  describe  the  variation  with 
that  parameter,  it  is  seen  that  Qp  can  be  a very  large  number  in  propel- 
lant with  mixed,  polydi6perse  oxidizer  having  "general"  particles,  Ae 
noted  previously,  computation  time  is  directly  proportional  to  Qp. 
Therefore,  minimizing  Qp  while  retaining  the  essence  of  the  problem 
is  computationally  important.  Consequently,  it  shall  be  assumed  that 
the  particles  are  ;,.t  most  ellipsoids  of  revolution.  This  means  that 
"orientableness"  is  retained  but  only  one  angle  is  required  to  distin- 
guish that  orientation  and  that  asphericness  is  retained  but  only  two 
dimensions  are  required  to  characterize  size  and  shape.  Therefore, 
particle  characterization  parameters  are  at  most  specie(s),  size  (D), 
aspect  ratio  (a),  and  orientation  ( ).  The  size  of  the  ellipsoid  (D)  is 
taken  to  be  the  maximum  dimension  normal  to  the  axis  of  symmetry. 

The  single  sum  in  Eq.  90  is  then  replaced  with  the  quadruple  sum  (Q„  = 
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is  a distribution  function  such  that 
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If  Np^  s t a,  d , ^ is  the  numbe  r of  partic  les  on  Sp  with  £>-  D * 
a.sa  s i)  &•  <9*  and  specie  k per  unit  area  of  Sp, 


and  ^ ^ is  the  average  planar  surface  for  fuel  surface/ 

oxidizer  particle  pairs  with  s,  a,  D,  ^ parameters,  Eq . 92  becomes 
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Passing  to  the  limit  as  -»  oo  and  Aa.y  -»o  gives 
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For  spherical  particles  /m  A'a  is  independent  of  a and  $ so  that 
integration  on  these  variables  gives 
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Equation  95  is  a statistical  expression  relating  the  mean  burning  rate  of 
a propellant  with  mixed,  polydisperse,  aspheric  oxidizer  to  the  burning 
rate  of  monodisperse  psuedo-propellants . If  there  are  no  interactions 
among  fuel  surface /oxidize r particle  pairs,  the  monodisperse  pseudo- 
propellant rates  can  be  computed  one  at  a time  using  any  monodisperse 
propellant  code.  In  what  follows  it  is  assumed  that  there  are  no  inter- 


The  next  step  is  to  investigate  the  statistical  characteristics  of 
the  burning  surface  and  thereby  relate  A distribution 

function,  and  properties  of  the  monodisperse  psuedo-propellants  to 
propellant  formulation  variables.  From  formulation  variables  the  mass 
fraction  of  particles  in  the  propellant  with  specie  s,  tit  4 and 

DsD  S D+dD  is,  in  principle,  known  from  measurements.  * 

.Therefore,  the  volume  fraction  of  these  particles  is 
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and  the  number  of  these  particles  per  unit  volume  of  the  propellant  is 
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where  ^ is  the  volume  of  a particle  with  aspect  ratio  a and  size 

D. 


The  size  of  the  ellipsoid  is  taken  to  be  the  dimension  normal  to 
the  axis  of  revolution  as  this  represents  the  smallest  dimension  the 
particle  will  pass  through.  Therefore,  in  particle  centered,  cartesian 
coordinates  '*1  i CJ’,  ^ the  equation  describing  the  particles  surface 
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Consequently,  the  volume  of  a particle  with  a,  D characteristics  is 
a-Wi  0.0/j. 
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* This  information  must  be  known  for  if  we  don't  know  what  goes  in  we 
cannot  expect  to  compute  what  comes  out. 
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Assuming  henceforth  that  the  statistics  of  the  particle  arrange- 
ment in  the  condensed  phase  are  homogeneous  and  isotropic,  all  orienta- 
tions of  aspheric  particles  are  equally  probable.  Therefore,  the  number 
of  particles  per  unit  volume  with  s,  a,  D and  q i s is* 
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However,  particles  with  f = <9  and  <P  =Tf-  <$  are  indistinguishable 
because  of  the  assumed  particle  symmetry.  This  double  degeneracy  is 
accounted  for  by  doubling  the  probability  and  limiting  to0<<^  4 tr/2,  . 
Thus,  with  Eq.  100 
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With  number  per  unit  volume  known  the  next  question  is  "what 
are  the  numbers/unit  surface  on  the  burning  surface?"  Assuming  with 
BDP  that  the  burning  surface  is  a plane  dotted  with  concave  and/or 
convex  particles,  we  need  only  consider  the  particles  intersecting  a 
plane  Sp.  Consider  the  relationship  between  a single  particle  and  the 
plane  Sp.  Let  denote  the  angle  between  the  normal  to  Sp  and  the 
axis  of  revolution  of  the  ellipsoid.  Denote  the  location  of  the  particle 
by  the  coordinates  of  its  CG.  Clearly,  all  particles  with  cj.  = ! Xtt  d ^ £ 
will  intersect  Sp.  The  volume  containing  s,  a,  d,  <!)  particles  that  can 
therefore,  AV  - so  that 


intersect  Sp  is 
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The  next  question  is  how  many  a,  D,  dj>  particles  intersect  Sp 
at  a specific  distance  from  the  point  of  initial  intersection  with  Sp. 

That  is,  how  many  particles  intersect  Sp  with  Z‘ZsX^dX.  Tointer- 
sect  Sp  under  these  conditions  the  CG  of  the  particle  must  lie  in  a volume 
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Therefore,  the  number  of  a, 
is 
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This  shows  that  all  0^1  s X., 


are  equally  probable. 


The  distance  Xa^  ^ ^ is  the  normal  distance  between  two  planes 
parallel  to  Sp  and  tangent  to  a particle  with  a and  D.  At  the  first  quadrant 
tangency  point  (*)  A%/A^  = - ^ <5j  . Differentiating  Eq.  98  with 
respect  to  ^ (set  cr  1 o )*,  solving  for  <A^  / d , setting  that 


^Consider  a sphere  with  unit  radius  and  = being  angle  between  radius  vector 
and  north  polar  axis.  As  orientations  are  eq ually  probable,  probability  of 
orientation  with  angle  ^ is  (area  zone)/(area  sphere)  - ^4  / 4 if 

**Since  the  particle  is  a body  of  revolution,  there  is  no  need  to  work  in  three 
dimension. 
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so  that  from  Eqs.  108,  104,  and  112 
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Within  the  context  of  the  BDP  model  this  dimension  alone  would 
describe  characteristic  distances  from  the  initial  point  of  tangency  of 
the  particle  with  Sf  that  would,  in  turn,  define  the  doubly  degenerate 
mean  state.  However,  as  pointed  out  in  the  review  of  the  BDP  model 
this  approach  i6  correct  only  in  the  case  of  no  ignition  delay. 


The  mean  statistical  characteristics  of  the  fuel  surface  in  each 
monodisperse  "cut"  of  fuel  surface /oxidize r particle  pairs  cannot  be 
determined  exactly  without  detailed  knowledge  of  the  particle  packing. 
Since  this  information  is  unavailable  ' , an  approximation  is  introduced. 
Examination  of  the  characteristics  of  regular  geometric  packings(29) 
suggests  that  the  mean  volume  of  "binder"  associated  with  a specific 
particl_e  should  be  roughly  proportional  to  that  particles  surface  area 
or  cC  AS>6<  a ^ . To  allow  for  variability  it  is  assumed  that 
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where  n is  an  empirical  parameter  (to  be  defined  by  correlation  with 
data)  and  C is  the  proportionality  constant. 

The  surface  area  of  an  ellipsoid  of  revolution  is 
D/J 
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Differentiating  Eq.96  to  find  and  employing  v-  ’(j)-''!  as  a 

change  of  variable  enables  the  resulting  expression  to  be  integrated 
in  closed  form  yielding 
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The  volume  fraction  of  binder  is  - l-  (propellant  assumed 
here  to  be  solely  binder  and  oxidizer).  However,  it  is  also 
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With  Eqs.  102  and  114  this  defines  C as 
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With  C and  n known  the  mean  volume  of  binder  and  oxidizer 
associated  with  monodisperse  pseudo -propellants  with  s,  a,  d, 
haracteris tics  are  known.  Therefore,  the  volume  fraction  of  oxidizer 
in  the  pseudo-propellant  is 
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With  Eqs.  99  and  114 

this  is 
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Since  pseudo-propellants  are  oxidizer  and  binder 
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Translate  a plane  Sp  in  the  propellant  a distance  . The 
volume  swept  out  is  A^  . The  fraction  of  this  swept  out  volume 
that  is  oxidizer  is  A-j,  . However,  the  fraction  of  this  volume 

that  is  oxidizer  is  also  ^ 0l  , Consequently, 
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Therefore,  in  the  pseudo-propellant  (actually,  this  is  a definition) 
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"here  the  bar-overs  denotes  means  for  fuel  surface /oxidizer  particle 
pairs.  Therefore,  withEq.  113 
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Attention  in  this  program  has  been  largely  directed  at  additive 
free  propellants.  However,  most  practical  propellants  contain  additives. 
Moreover,  Cohen,  Derr,  Priced4  ) have  demonstrated  that  rate,  pres- 
sure characteristics  of  metallized  propellants  can  be  predicted  by 
assuming  the  additive  is  chemically  inert.  Therefore,  it  appears  that 
a useful  approximation  to  real  propellant  behavior  might  be  made  by 
assuming  that,  insofar  as  rate  defining  processes  go,  additives  are 
inert  and  are  apportioned  among  the  pseudo-propellants  as  a fixed 
fraction  of  binder.  Subdividing  binder  into  fuel,  catalyst,  and.  metal 
= <iV(  * dV^  * cW^  . Therefore,  /d\Jv  --  / ],  , 

i*  ! Tv  1 I • Consequently,  in 

any  pseudo-propellant 
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The  mass  fraction  of  oxidizer  in  the  pseudo  propellant  is 
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With  pseudo-propellant  density  known  the  mass  fractions  can  be 
computed  as 
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The  specific  heat  of  the  condensed  phase  is 
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Treating  additives  as  inert  the  flame  temperature  for  any 
pseudo -propellant  can  be  computed  from  a thermochemistry  code. 
However,  the  number  of  such  computations  required  to  define  flame 
temperature  for  all  possible  combinations  of  oxidizer,  metal,  fuel, 

< atalyst  and  mixture  ratio  is  excessive.  Therefore,  an  approximate 
strategy,  valid  for  small  additive  concentrations,  is  employed.  Basic- 
ally the  adiabatic  flame  temperature  of  just  fuel  and  oxidizer 

«.  " is  computed  with  the  thermochemistry  code  and 

tabulated.  This  can  be  done  once  and  for  all  for  a specific  fuel,  ox 
combination.  The  inerts  are  then  added  as  a diluent.  Consequently, 
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where  c'  is  the  specific  heal  at  constant  pressure  of  the  fuel  and  oxidizer 
and  Q,  the  latent  heat  of  fusion  of  the  j"1  inert  material. 

The  volume  fraction  of  oxidizer  with  s,  a,  d characteristics  is 
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Oxidizer  particles  in  a propellant  are  the  result  of  mixing  several  modes 
of  oxidizer  together.  Each  mode  is  characterized  separately  as  to  (a) 
the  distribution  of  sizes  (D)  within  that  mode  and  (b)  the  specific  surface 
in  that  mode.  Therefore, 
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is  known  by  measurement.  The  surface  area  in  a mode  is 
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where  the  surface  area  AS  ^ ^ of  a specific  particle  is  given 

by  Eq.  113.  The  mass  of  particles  in  the  k*-^  mode  with  s,  a,  d charac- 
teristics  is  related  to  the  number  with  that  characteristic  by 
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or  with  Eq  . 99 
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Thus,  Eq.  134  becomes 
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rherefore,  to  employ  surface  area  measurements  to  define  parameters 
in  the  modes  distribution  of  size  and  shape  the  distribution  of  shape 
must  be  known  because 
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This  information  is  just  not  aailable.  Consequently,  we  are  in  a posi- 
t.  n where  the  theoretical  model  is  more  general  than  the  available  input 
information.  Hence,  to  proceed  further  some  assumption  relative  to  the 
distribution  of  particle  shape  for  specific  size  is  required.  Assume 
this  distribution  is  log  normal  about  some  mean  as  j,.  with  standard 
deviation  ira  s [<.  Then 
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Consequently,  the  surface  area/unit  mass  of  oxidizer  mode  s,k  is 
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This  equation  leaves  one  free  parameter  for  each  mode.  This  parameter 
can  be  determined  by  correlation  to  a data  base. 

With  ae  k and  u g ^ known  in  principle, 
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where  Ks  is  the  number  of  oxidizer  modes  with  species  s.  Therefore, 
cot  • Eqs.  132,  140,  and  142  yields 
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With  Eqs.  107  and  143  this  becomes 
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Therefore,  with  Eq,  147,  Eq.  95a  becomes 
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For  spherical  particles  In  a = f*va  and  m"*  i func  (<f).  Therefore, 
integration  with  respect  to  <f  from  (^  = o to  ^ = TT/2  gives 
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rhis  is  the  resuit  previously  reported  for  spherical  particles  in  Ref.  28. 

Equation  148  relates  the  mean  burning  rate  of  composite  propel- 
lant with  n ixed,  polydisperse  spheroidal  oxidizer  to  the  burning  rate 
of  a sequence  of  monodisperse  psuedo-propellants.  The  properties  of 
these  monodisperse  psuedo-propellants  depend  upon  the  specified  pro- 
pellant ingredients.  The  remaining  theoretical  problem  is  to  compute 
the  burning  rate  of  a monodisperse  psuedo-propellant.  This  will  be 
handled  herein  with  a modified  BDP  model. 

MONODISPERSE  COMBUSTION  MODEL 

The  BDP  model  has  been  reviewed  previously;  several  errors 
were  noted.  These  are  listed  below; 

1.  mean  burning  rate  based  on  wrong  surface  area, 

2.  characteristic  dimension  for  oxidizers  mean  deflag  rating 
state  based  on  the  ignited  and  unignited  particles. 

characteristic  fuel  dimension  based  on  ordered  particle 
arrangement  when  particles  are  really  randomly  ordered. 

4.  Final  diffusion  flames  influence  on  the  A/PA  flame  was 
neglected. 

The  discussion  here  deals  solely  with  spherical  particles.  Portions 
dealing  with  ellipsoidal  particles  are  found  in  Appendix  II 

Surface  Area  Base  - Applying  conservation  of  mass  to  the  control 
volume  bounded  by  SQ  and  Sp  yields 
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Since  BDP  compute  rQ  the  correct  rate  is 

r = r (S  /S  ) 

p o o p cJs0) 

where  S /S  is  available  from  their  code, 
o p 

Characteristic  Dimension  for  Mean  Deflagrating  State  - A cen- 
tral problem  in  statistical  combustion  modeling  is  to  define  a particles 
mean  deflagrating  state.  It  has  been  shown  herein’1  that  for  a spherical 
particle  the  number  of  particles  intersecting  a planar  surface  is 


L A*  £) 


and  that  all  depths  of  intersections  are  equally  probable  or 

For  a monodisperse  propellant  d a A * li„.  . For  spherical  particles 
a = 1 and  ^ i D . Therefore,  integration  of  Eq.  143  from  ^ = o 

t’>  ^ = "n’/2  gives  the  result  shown. 
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The  BDP  model  assumes  that  binder  regression  rate  ( r f ) is  stationary 
while  oxidizer  regression  rate  is  a step  function.  Since,  the  ensemble 
of  exposed  particles  contain  samples  of  all  states  in  a particles  lifetime 
and  statistics,  rates,  and  initial  particle  geometry  known,  all  micro- 
states are  known  because  they  are  sequential.  Therefore,  mean  states 
can  be  computed  from  sums  over  all  states  on  the  surface.  This 
approach  to  definition  of  means  is  termed  "petite  ensemble  averaging" 
and  follows  the  spirit  of  integral  methods. 

The  deflagrating  surface  of  an  oxidizer  particle  is  assumed  to 
be  a spherical  segment.  Referring  to  Fig.  12  application  of  the 
Pathagorean  theorem  to  triangle  OabO  yields 
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Application  to  triangle  O'abO  gives  (with  Eq.  153) 

+•  x[hio)-x"\ 


The  distance  h is  the  difference  between  the  distances  traversed 
by  binder  and  oxidizer.  Therefore,  during  the  ignited  phase 
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If  h < o,  the  particle  is  consumed  before  x = D (o).  Therefore, 
oi  X i he)  + £l«)  * X-  e rf  be)  /rox 


Consequently,  for  all 
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mil  restates  will  contain  oxidizer  while  those  with  s X { 0(°j 

will  not.  A microstate  is  empty  lor  At  * ^Otoi  - ^ / r(  . During 

this  time  interval  (a)  continued  pyrolysis  may  expose  "new"  oxidizer 
particles  in  the  empty  microstate  or  (b)  the  empty  microstate  may 
"till  up"  with  surface  melt.  Although  either  (a)  or  (b)  are  really 
dvnarm  phenomena*,  clearlv  a surface  melt  is  required  for  (b)  and 
' t r v surfac  e guarantees  (a).  Since  melt  layer  thickness  dec  reases 
with  increasing  pressure,  the  probability  of  (a)  occurring  increases 
with  increasing  pressure.  It  should  also  increase  with  increasing 
particle  size.  These  "trends"  follow  the  behavior  observed  by  Cohen 
if.  nitramine  propellants . ( ^ * Situation  (a)  means  the  burning  surface 
is  no  longer  quasi-planar  (a  basic  assumption  in  this  study).  Therefore, 
"nly  (b ) will  be  considered  herein. 
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On  the  other  hand,  if 
♦ r;  D(o)/r0i  > , 

the  tinal  microstate  v.  ill  possess  oxidizer, 
residual  oxidizer  is 
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The  initial  size  of  this 
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Tiles e particles  will  be  either  (c)  freed  from  the  surface  (blown  away) 
or  (d)  retained  on  the  surface.  The  choice  of  (c)  or  (d)  also  hinges  on 
the  pressure  of  a surface  melt.  If  a melt  is  absent,  (c)  is  probable. 

If  a melt  is  present,  (d)  is  probable.  Since  residual  particles  occur 
with  large  particles  at  lower  pressures,  a melt  is  probably  present. 

Note  that  if  residual  particles  are  retained,  the  monodisperse  propellant 
becomes  effectively  bidisperse.  Both  of  these  situations  will  be 

considered.  -0000 ^ 

Considering  now  the  normally  accessible  microstates,  the  sur-  ^ 
face  area  of  a spherical  segment  is 
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Combining  Eqs.  155,  156,  and  161  yields 

T [®‘  iur  ro.0-ft,/r(')/r{  f Ot»)]x  *(r,Xr  /rjj  Out) 


The refo r<  , 


f 

\ AS 

— 

J 0VV> 

*'V 

T 

s how  long 

does  i 

the  site  and  the 

— Ax 


C".J) 


58 


With  Eqs.  151,  152,  and  162  this  gives 
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The  planar  area  of  an  oxidizer  surface  is 
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Therefore,  the  fraction  of  planar  area  occupied  by  ignited  particles  is 

^ S>,°VV"  / ^ J'F“  J ^ Jf,0‘  cisp4n  <JU) 

X'Y* 

With  Eqs.  151,  152  and  165  this  becomes 
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The  number  of  ignited  particles  is 
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Therefore,  the  mean  diameter  of  the  ignited  particles  is 

t,  - trtt  to  i -c*n  -Cfj] 
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Note  that  the  BDP  value  is  not  achieved  if  either  itl<¥„  or  If 

< D(°l  . Figures  13  and  14  v»*  S1V, /0(«jand  /d.% 

function  of  % .<*.  / M°i  for  several  rox/rf.  Appropriate  BDP  values 

are  indicated  for  comparison.  Note  that  significant  discrepancies  can 

occur. 


As  mentioned  previously,  the  situation  where  oxidizer  depleted 
microstates  occur  is  to  be  handled  at  present  by  as s urning  that  melt 
fills  the  site.*  The  situation  where  residual  oxidizer  exists  at  the  final 
microstate  can  be  treated  because  the  quasi-planar  surface  assumption 
is  not  violated. 


If  the  burning  surface  is  dry,  the  residual  particle  is  "free"  and 
will  blow  away.  Thus,  the  oxidizer  flowing  away  from  the  burning  sur- 
face consists  of  both  solid  and  gaseous  phases.  Application  of  mass  con 
servation  on  both  global  and  specie  bases  yields 
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where  mj^  , is  the  mass  flux  of  gaseous  oxidizer  and  mQX  r is  the  mass 
flow  rate  of  residual  particles.  Since  only  ignited  particles  evolve  gas, 
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Combining  Eqs.  171  and  172  and  employing  the  mean  value  theorem  for 
integrals  yields 


'°r 


OX, 


(.'■>3) 


The  mass  flow  rate  of  residual  particles  is  given  by  the  product 
of  their  rate  of  creation  and  a residual  particles  initial  mass.  Since 
the  BDP  model  assumes  all  deflagrating  oxidizer  surfaces  are  segments 
of  a sphere 


A rrf^t 

°*lr 


TPox 


OTM) 


la 


.CtlOOfeV. 


-4  1—  ilM 


Consider  now  a planar  at  t and  t + dt.  Since  the  distance  traversed  is 
rdt  the  propellant  volume  swept  out  is  dV  = r S dt.  The  number  of 
attached  ixidizer  sites  with  their  south  poles  in^his  volume  is  dN  = 
(dNox/dV)  dV  where  dkl /A V ■=■  <•  \OI  /\j  . Therefore,  the 

rate  of  creation  of  residual  particles  is 

= fo  \6).  r h’co)]  Ots) 

*This  is  a "crutch";  more  work  on  this  aspect  is  needed. 
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For  situations  where  h[D  (o)]/D  (o)  so  06  =06 

' 0 Ol 

On  the  other  hand,  if  the  burning  surface  is  wet,  the  residual 
particles  will  probably  stick  to  the  surface  and  be  consumed  thereon. 
This  presents  appreciable  theoretical  difficulties  (not  insurmountable) 
because  these  stuck  residual  particles  alter  the  population  of  particles 
that  exist  on  the  burning  surface.  This  means  that  distribution  functions 
and  pseudo-propellant  properties  are  no  longer  defined  solely  by  the 
propellant  recipe.  Environmental  variables,  since  they  influence  the 
existence  and  characteristics  of  the  residual  particles,  also  influence 
distribution  functions  and  pseudo -propellant  properties.  Appendix  I 
presents  the  treatment  (incomplete)  of  the  stuck  residual  particles  with 
a P.DP  like  model. 
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Fuel  Surface  Considerations  - The  fraction  of  Sp  that  is  fuel  is 
. Since  the  se  are  (ik/d Sp  - & ^ /Qrrcf ) oxidizer  sites  on  Sp  4 
the  mean  fuel  surface  per  site  is 


- TT 
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Arranging  this  fuel  surface  increment  in  an  annulus  about  the  mean 
ignited  particle  gives 


TT^=  AS,  t it  5 /q. 
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Employing  Eqs.  Ho  and  and  solving  for  b gives 
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The  effective  oxidizer  volume  fraction  0 has  been  introduced 

•'ox 

to  account  for  the  fact  that  the  effective  mixture  ratio  of  the  combus- 
tion reactions  are  altered  when  partially  consumed  oxidizer  particles 
escape  unreacted.  Since  ^ /p>i 
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Flame  Considerations  - The  gas  phase  temperature  field  for 
OS  X s Xpp  is  divided  into  two  parts.  The  final  flame  part 

s x s xrF  and  the  AP  part  os  xs  x*p  . For  the  approximation 
employed  by  Ref.  £ energy  release  is  concentrated  at  the  flame 
locations.  An  energy  balance  at  *pf.  gives 
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Since  TC*rF)  * integration  of  the  energy  equation  from 
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Hi  ve  s 


An  energy  balance  at  X 


% (**.)  = * ?Vv> 

Integration  of  the  energy  equation  from  6^  to  0 gives 
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Comparison  of  Eq.  (6)  of  Ref.  3 with  Eq.  (187)  shows  that  the  heat 
feedbacks  from  the  AP+FF  flame  sequence  is  identical.  Consequently, 
Eq.  (6i  of  Ref.  3 does  not  need  alteration.  The  only  differences  are 
that  i (*Rp)  is  not  the  adiabatic  flame  temperature  for  AH  and  T ( 0 ) 
is  available  for  comparison  with  Ts.  The  former  means  that 
in  Eq.  (18)  of  Ref.  3 should  be 


The  latter  means  that  a criteria  for  apportioning  energy  release 
between  condensed  phase  and  AP  flame  is  available. 

With  the  above  modifications  the  AP  flame  is  not  adiabatic. 
Therefore.  Eq.  20  of  Ref.  3 becomes 

where,  as  discussed  by  Ref.  3,  - BID  rdL/fynj 
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Results  and  Discussion 

Numerical  experiments  were  made  using  the  BDP  model  to  deter- 
mine its  ability  to  predict  mixture  ratio  effects.  Figure  15  presents 
results  for  both  the  original  BDP  model  and  a version  modified  to 
employ  the  fuel  dimension  based  on  random  particle  arrangement. 

These  results  clearly  show  that  the  BDP  model  exhibits  singular  be- 
havior at  high  and  low  mixture  ratios.  The  singular  hehavior  was 
traced  to  the  Burke-Srhuman  solution  for  diffusion  flame  height.  Two 
fixes  were  instituted.  First,  the  Burke-Schuman  solution  was  replaced 
by  a diffusion  flame  height  based  on  dimensional  arguments.  This 
removed  the  aforementioned  singularities  and  decreased  runtime. 
Second.  J.  A.  Condon  has  found  that  the  singular  behavior  can  be  over- 
come by  employing  more  terms  in  the  series  solution  present  in  the 
Burke-Schuman  solution.*  Near  stoichiometric  conditions  the  "single 
term  solution"  employed  by  BDP  is  adequate.  However,  at  extreme 
rich  or  lean  conditions  as  many  as  100  terms  may  be  required. 

Table  3 shows  the  effects  of  some  of  the  modifications  described 
herein  on  BDP  model  predictions  for  r(>)  of  the  monodisperse  PS/AP 
propellants  indicated.  In  this  table  BDP  (,r»)  denotes  the  original  BDP 
model,  BDP(rr)  denotes  the  original  model  with  rate  based  on 
rather  than  , BDP  (b)  denotes  the  BDP  model  with  the  b dimension 
based  on  random  rather  than  ordered  arrangement,  BDP  (B-S)  denotes 
the  BDP  model  with  the  multi-term  (501  Burke-Schuman  solution,  and 
BDP  (Geo.  ) denotes  the  BDP  model  with  4S  ; /cJs^. 

Figure  15  compares  results  for  BDP  CO  , BDP  (B-S),  and  BDP  (Geo.) 
graphically.  These  results  show  these  corrections  primarily  impact 
high  pressure  behavior. 

The  BDP  code  was  also  modified  to  include  residual  particles 
that  were  not  retained.  Numerical  results  at  high  pressure  indicated 
little  effect  (as  expected).  Numerical  results  at  low  pressure  with 
large  particles  "blew  up"  because  ■LLt>(o)l  / '■  This  result  was 
not  expected.  However,  it  merely  illustrates  what  is  known  from 
experiment.--Namely,  AP  particles  do  not  escape  from  the  deflagrating 
surface  unreacted.  Consequently,  theoretical  determination  of  resi- 
dual particle  effects  awaits  the  complete  development  and  inclusion  of 
retained  residual  particles. 

The  flame  modifications  noted  in  the  text  have  not  been  added  to 
the  code  yet. 

A code  for  polydisperse,  additive  free  propellants  with  spherical 
AP  oxidizer  was  assembled  and  made  operational.  In  this  code  it  is 
assumed  the  particle  size  distributions  in  each  oxidizer  mode  are  log 
normal.  That  is,  it  is  assumed  that 

♦Sammons  ^ employed  this  approach.  However,  Condon  has  found 
that  Sammons  employed  incorrect  values  for  the  appropriate  Bessel 
functions. 
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^Rates  in  in/sec. 


Appendix  III  presents  operating  instructions,  the  FORTRAN  IV  listing, 
and  a sample  problem  for  the  code. 
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T-BURNER  VENT  FLOW  STUDY 


Results  in  this  study  have  been  reported  in  the  open  literature 
(see  Ref.  30  and  31).  Therefore,  it  will  not  be  redescribed  in  detail 
herein.  However,  certain  comments  are  in  order. 

The  hydraulic  analogy  studies  demonstrate  that  flow  in  the  vent 
region  is  a sequence  of  flow  attachment  and  separation  from  alternate 
faces  of  the  vent.  This  process  creates  a Karman  vortex  street  down- 
stream of  the  vent  entrance.  Therefore,  the  available  acoustic  energy 
(kinetic)  at  the  vent  entrance  must  be  substantially  transformed  into  the 
kinetic  energy  of  the  vortex  system.  Consequently,  little,  if  any, 
amplification  appears  to  occur  at  the  vent. 

The  vent  gain  is  the  reverse  of  the  flow  turning  loss.  However, 
the  flow  turning  loss  is  an  irreversible  process.  Therefore,  Culick's 
analysis  "appears"  to  violate  the  second  law  of  thermodynamics. 
However,  recent  results  by  Culick  "appear"  to  refute  the  arguments  of 
Ref.  31.  Consequently,  little,  if  any,  amplification  appears  to  occur 
at  the  vent. 


FUTURE  PLANS 


This  work  has  formulated  a viable  general  approach  to  the  com- 
bustion of  composite  propellants.  However,  specific  application  of  the 
methodology  (operational  computer  code)  has  been  to  polydisperse,  addi- 
tive free  AP  propellants  with  spherical  oxidizer.  Moreover,  that  code 
does  account  for  neither  retained  residual  oxidizer  articles  nor  multiple 
particle  penetrations.  Consequently,  the  general  direction  of  future 
efforts  is  reasonably  clear. 

o Implement  retained  residual  particle  methodology 

o Implement  flame  modifications 

o Develop  methodology  to  treat  multi-particle  penetrations 
o Extend  modeling  to  nitramine  "oxidizers" 
o Implement  metal  and  catalyst  additive  methodology 
o Develop  methodology  to  treat  particle -particle  interactions 
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NOMENCLATURE 


LA  TIN  SYMBOLS 

a aspect  ratio  of  ellipsoidal  particle  (length/diameter) 

b D + S 


c 

C 

D 

e 

E 

F 

h 

k 

K 

s 

L 

m 

m 

m' 

n 

N 

P 

q" 

Q 

r 

R 

S 

t 


specific  heat 

constant  defined  by  Eq . 118 
diametei 

specific  internal  energy 
internal  energy 
distribution  function 

specific  enthalpy  or  height  of  spherical  cap  of 
oxidizer  particle 

Thermal  conductivity  ^ 

number  of  oxidizer  modes  associated  with  s oxidize 
specie 

normal  distance  from  plane  S to  plane  S 

P ° 

mass 

mass  flow  rate 
mass  flux 

empirical  exponent  defined  by  Eq. 
number  of  particles,  number 
pressure 
heat  flux 

number  of  microstates  or  energy  to  gasify 

burning  rate 

radius 

surface  area 
time 
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r 


X 

y 

GREEK  SYMBOLS 


temperature 

velocity 

volume 

mass  fraction 

spatial  coordinate  as  noted 

spatial  coordinate  as  noted 

spatial  coordinate  as  noted 


OC 

h 

\ 

J 

S 

A 

A. 


"1 

^ % 


mass  fraction,  mass /mass  propellant 
intercept  defined  by  Eq.  1 06b 
volume  fraction,  volume/volume  propellant 
mean  width  of  fuel  around  oxidizer  particle 
liquid  layer  thickness 
denotes  an  increment 

distance  from  solid/liquid  interface  to  center  of  curva- 
ture of  a retained  residuals  burning  surface 

denotes  an  area  element  or  ratio  of  outer  diameter 
of  excluded  zone  to  particle  diameter 

parameter  in  Eq.  89 


SPECIAL  SYMBOLS 


orthogonal  coordinates  as  noted 
density 


( ) 

* 

( ) 

SUBSCRIPTS 


denotes  a mean  or  value  for  a particle 
denotes  a monodisperse,  psuedo-propellant 

denotes  particle  aspect  ratio  or  particles  with 
denotes  binder 
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c 

Cait 

* 

:r 

A 

I 

N 

0 

on 

1r 

1 

ir 

s 

<¥ 

CD 

t 


denotes  condensed  phase  or  characteristic: 
denotes  catalyst 

denotes  oxidizer  particles  with  s Cl  - 
denotes  fuel 

denotes  ignited  particles 
denotes  j*-^1  oxidizer  particles 
denotes  k1*1  oxidizer  model 
denotes  metal 
denotes  northern 
denotes  burning  surface 
denotes  oxidizer 
denotes  planar  surface 
denotes  residual  particles 
denotes  heterogeneous  reaction 
denotes  southern 

denotes  particles  with  *■  ^ 

denotes  initial  propellant  temperature 
denotes  point  of  tangency 
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APPENDIX  I 


RESIDUAL  PARTICLES 


When  residual  particles  are  retained,  something  holds  them  to 
the  burning  surface.  The  most  likely  candidate  is  a surface  melt.  That 
is,  the  oxidizer  particles  are  stuck  in  a surface  melt.  If  this  is  correct, 
the  evolutionary  trend  in  the  combustion  of  the  stuck  residual  will  be 
something  like  that  sketched  in  the  margin.  Note  that  if  an  interfacial 
phalanx  flame  "frees"  the  particle  it  will  simply  float  away.  Indeed,  it 
would  not  have  become  stuck  in  the  first  place.  This  suggests  that 
interfacial  flames  may  not  be  important. 
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Consider  now  a specific  evolutionary  state  in  some  detail, 
equation  of  the  original  residual  particle  is* 
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layer  must  lie  on  this  surface  so  that 
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Therefore,  a and  hence  the  mean  intersection  diameter  r = 2a  are 

known  and  " * 
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Note  that  in  this  situation  all  residual  particles  are  ignited. 


Following  BDP  the  deflagrating  surface  is  assumed  to  be  always  / 
a segment  of  a sphere.  Therefore,  (see  marginal  sketch) 


1-4 


where  is  the  radius  of  burning  surface  of  the  residual  particle  and 
(°»  ) a-e  the  coordinates  of  the  center  of  that  spherical  segment. 

This  surface  must  pass  through  (O,  h [D  (o )]  - y)  and  (a,  £L  ).  This 
gives  two  equations  for  and  Rf.  Solving  gives 
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* Only  two  dimensions  need  be  considered  because  the  particle  is  a 
body  of  revolution. 
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The  area  of  a spherical  cap  is 


AV  . 2TTR  -K 


Substituting  Eqs.  1-5  and  1-6  this  becomes 


A^r  --  *TT  \ + WM  - /l 


Assuming,  as  in  the  BDP  model,  that  instantaneous  oxidizer  regression 
rate  is  constant  in  time  for  a specific  particle,  all  microstates  are 
equally  probable.  Therefore, 

ffe,  • ^ % 


Consequently, 
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Employing  Eqs.  1-8  and  1-9  and  integrating 
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From  Eq . 174  the  creation  rate  of  residuals  is 
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Since  mox>  r (o)  = IT  po>  [ G,  the  mass  flux  of  solid 

residual  oxidizer  being  created  at  the  burning  surface  is 
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The  rate  at  which  residual  oxidizer  is  being  gasified  is 


( | f - r ^ ^ \\,r  - tt  r n ^ 
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At  steady-state  these  rates  are  equivalent  so  that 
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Equations  1-3  and  1-16  define  the  geometry  of  the  residual  particles  for 
a residual  BDP  model. 


The  existence  of  retained  residual  particles  means  there  are  more 
particles  on  the  burning  surface  than  one  computes  from  the  propellant 
ingredients.  In  the  petite  ensemble  approach  each  distinguishable 
particle  is  recognized  in  the  statistical  framework  and  a monodisperse 
psuedo-propellant  is  generated  for  that  distinguishable  particle. 

Therefore,  retained  residual  particles  mean  that  the  statistics  of  the 
burning  surface  are  altered.  Although  this  aspect  of  the  problem  has 
not  been  explored  at  this  time,  no  great  "theoretical"  difficulty  is 
forseen  (just  time).  However,  as  the  distribution  function  is  now  depen- 
dent u£on  the  combustion  model,  the  present  model  becomes  the  interior 
01  an  iteration  loop  involving  the  distribution  function.  Consequently, 
’exact"  treatment  of  retained  residuals  will  probably  triple  or  quadruple 
run  time. 
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APPENDIX  II 


ELLIPTICAL  PARTICLES 


A central  difficulty  with  aspherical  particles  is  finding  an 
expression  for  the  deflagrating  surface  of  the  particle  that  satisfies 
the  following  criteria:  deflagrating  surface  possess  same  intersec- 
tion with  Sp  as  ellipsoidal  oxidizer  particle,  deflagrating  surface  lies 
within  ellipsoidal  particle,  and  deflagrating  surface  possess  a tangent 
point  within  ellipsoidal  particle  to  Sp.  This  "exercise"  in  Analytical 
Geometry  appears  reasonably  6imple.  However,  it  has  proven  to  be 
difficult. 

The  equation  for  an  ellipsoid  of  revolution  is 
(r [+t)/<£  + *-  1 I- 1 

In  an  ^ - constant  plane  assume  that  deflagrating  surface  has  the 
form 
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A and  B are 
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where  % denotes  the  value  of  g at 
the  value  of  % at  the  intersection  of 
Solving  for  A and  B gives 
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The  curve  in  the  O'  = O plane  must  approach  an  ellipse  as 
' -*■  ° . Thus,  the  simplest  possible  curve  must  approach  a 

segment  of  an  ellipse.  Assuming  that  this  ellipse  is  translated  and 
stretched  but  axis  of  symmetry  is  always  parallel  to  the  particles 

+ c 4 - 

Then  - aj  and  A.  = so  that 
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This  curve  must  pass  through  the  intersection  of  with  the 
particle  at  CT»o  ('f\^J%x  > ^ ) • Thus,  (ft,  = Ax  =a£) 

A,  ('V|1  - v ) * Ax  ( £r  - a,  - I 2-la, 


Solving  for  Aj  and  A£  yields 
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Differentiating  Eq.  II. 6 with  respect  to  7|_  and  setting  that  \| 

result  to  - W;  ^ gives  Ji.  £ 
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Therefore,  since  the  point  where  d^/d1^  - - dev*,  must  lie  on  ! Ei- 

this  line  the  coordinates  of  the  tangency  point  are 
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Since  s must  also  be  on  the  ellipse 
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Solving  for  s gives 
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As  ^ lies  on  the  particles  surface 


so  that 
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The  T - C curve  is  given  by  Eq.  11-6  so  that 
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Since  and  (-*\  is  known,  the  equation  for  the 

deflagrating  surface  is  known.  * 
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CARD  1:  Data  Card 


NUMBER  REQUIRED:  One  per  run 

FUNCTION:  Specify  number  of  oxidizer  types,  pressure  range,  print 
parameter,  accuracy  parameter  in  Bessel  series  expan- 
sion and  limit  on  number  of  terms  in  Bessel  series, 

FORMAT:  (15,  2F10.5,  215,  F10.5) 

Columns  1-5:  NOX,  number  of  oxidizer  species 

Columns  6-15;  PSTART,  pressure  to  start  incremental  calcu- 
lations, atm 

Columns  16-25:  PSTOP,  pre  ssure  to  stop  incremental  calcu- 
lations, atm 

Columns  26-30;  IPRINT,  = 0 suppresses  OUTPUT  of  rate  v.  s. 
D.  , IPRINT  = 1 prints  v vs.  Do 

Columns  31-35:  LIMBES,  Limit  on  number  of  terms  on  Bessel 
series  solution 

Columns  36-45:  ERRBES,  minimum  error  allowed  in  Bessel 
series  solution 

CARD  2:  Data  Card 

NUMBER  REQUIRED:  One  per  run 

FUNCTION:  Specify  integration  parameters 

FORMAT:  (215,  3F10.01 

Columns  1-5:  NCOUNT,  number  of  intervals  in  the  numerical 
integration  of  total  propellant  mass  flux.  Typically  30 

Columns  6-10:  NXCOUN,  number  of  intervals  in  the  numerical 
integration  for  the  proportionality  constant  C in  the 
equation  for  the  volume  of  fuel  associated  with  a particle 
in  a polydisperse  packing.  Typically  300 


84 


Columns  11-20:  XN,  diameter  exponent  in  the  equation  for  the 

volume  of  fuel  associated  with  a particle  in  a polydisperse 
packing.  2.  0) 

Columns  21-30:  DSTART,  Particle  diameter  to  start  integration 
Columns  31-40:  DSTOP,  Particle  diameter  to  stop  integration 
CARD  3:  Data  Card 

NUMBER  REQUIRED:  One  card  per  run 
FUNCTION:  Specify  fuel  type 
FORMAT:  (15) 

Columns  1-5:  IFUEL,  IFUEL  = 1 if  HTPB,  IFUEL  = 2 if  PB AN 
CARD  4:  Data  Card 

NUMBER  REQUIRED:  One  card  per  run 
FUNCTION:  Specify  propellant  parameters 
FORMAT:  (6E12.  6) 

Columns  1-12:  TZERO,  initial  propellant  temperature,  deg  K 

Columns  13-24:  XALFA,  Oxidizer  mass  fraction 

Columns  25-36:  QFUEL,  heat  of  pyrolysis  of  the  fuel  binder, 
cal  /g 

Columns  37-48:  RHOF,  density  of  fuel  binder,  g/cm^ 

Columns  49-60:  AF,  arrhenius  frequency  factor  of  fuel  binder, 
g/cm^-sec 

Columns  61-72:  EF,  activation  energy  of  the  fuel  binder,  cal/mole 
CARD  5:  Data  Card 
NUMBER  REQUIRED:  One  per  run 
FUNCTION:  Specify  propellant  parameters 
FORMAT:  (6E12.6) 


Columns  1-12:  BETA,  mass  fraction  of  metal 


,A 
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Columns  13-24:  RHOM,  density  of  metal,  g/cm^ 

Columns  25-36:  QM,  heat  release  of  metal  combustion,  cal/g 

Columns  37-48:  PCSUBP,  specific  heat  of  propellant,  cal/g  °K 

Columns  49-60:  PLAMB,  thermal  conductivity  of  propellant, 
cal/cm-sec-°K 

Columns  61-72:  XK,  proportionality  constant  for  responce 
function  peak 

CARD  6:  Data  Card 

NUMBER  REQUIRED:  One  per  run 

FUNCTION:  Specify  oxidizer  type 

FORMAT:  (15) 

Columns  1-5:  IOXID,  IOXID  = 1 if  AP,  2 if  KP,  and  3 if  HMX 
CARD  7:  Data  Card 

NUMBER  REQUIRED:  One  card  per  oxidizer  type 
FUNCTION:  Specify  propellant  parameters 
FORMAT:  (6E12.6) 

Columns  1-12:  Not  Used 

Columns  13-24:  GMWD,  Molecular  weight  of  final  flame  products 

Columns  25-36:  XNOID,  Proportionality  constant  for  final 
diffusion  flame 

Columns  37-48:  XNUPD,  proportionality  constant  for  primary 
diffusion  flame 

Columns  49-60:  PMWD,  molecular  weight  of  primary  flame  products 
Columns  61-72:  Not  used 
CARD  8:  Data  Card 

NUMBER  REQUIRED:  One  card  per  oxidizer  type 
FUNCTION:  Specify  propellant  parameters 
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FORMAT:  (6E12.6) 

Columns  1-12:  QLD,  latent  heat  of  vaporization  of  the  oxidizer, 
cal/g 

Columns  13-24:  RHOXD,  density  of  the  oxidizer,  g/cm^ 

Columns  25-36:  AOXD,  arrhenius  frequency  factor  of  the  oxi- 
dizer g/cm^-sec 

Columns  37-48:  EOXD,  activation  energy  of  the  oxidizer  cal/mole 

Columns  49-60:  TAPD,  temperature  of  the  AP  flame,  deg  K 

Columns  61-72:  Not  used 
CARD  9:  Data  Card 

NUMBER  REQUIRED:  One  card  per  oxidizer  type 
FUNCTION:  Specify  propellant  parameters 
FORMAT:  (6E12.6) 

Columns  1-12:  CIGND,  oxidizer  ignition  delay  paramet  r, 
sec(atm)m  cm  where  m*POWIGN  and  n=  PCWD 

Columns  13-24:  POWIGD,  pressure  exponent  in  oxidizer  parti- 
cle ignition  delay  term 

Columns  25-36:  POWDD,  diameter  exponent  in  oxidizer  particle 
ignition  delay  term 

Columns  37-48:  CONFD,  CONF  = 0 if  parabolic  flame  assumed, 
CONF  = 1 if  conical  flame  assumed 

Columns  49-72:  Not  used 
CARD  10:  Data  Card 

NUMBER  REQUIRED:  One  card  per  oxidizer  type 
FUNCTION:  Specify  propellant  parameters 
FORMAT:  (6E12.6) 

Columns  1-12:  KPFD,  rate  constant  of  primary  flame,  g/cm^- 
sec  -atm) 
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Columns  13-24:  KAP1D,  rate  constant  of  AP  flame  at  low  pressure, 
g/(cm^-sec-atm) 

Columns  23-36:  KAP2D,  rate  constant  of  AP  flame  at  high 
pressure,  g/cm^-sec-atm) 

Columns  37-48:  XNlD,  reaction  order  of  primary  flame 

Columns  49-60:  XN2D,  reaction  order  of  AP  flame  at  low 
pressure 

Columns  61-72:  XN3D,  reaction  order  of  AP  flame  at  high 
pressure 

CARD  11:  Data  Card 

NUMBER  REQUIRED:  One  care  per  oxidizer  type 
FUNCTION:  Specify  propellant  parameters 
FORMAT:  (6E12.6) 

Columns  1-12:  CSUBPD,  average  heat  capacity  of  solids  and 
gases,  cal/g-°K 

Columns  13-24:  XLAMBD,  average  thermal  conductivity  of  the 
combustion  gases,  cal /cm- sec-°K 

Columns  25-36:  GAMMAD,  diffusion  parameter,  cm^/sec 

Columns  37-48:  AFHD,  flame  height  factor 

Columns  49-60:  EPSD,  exponent  for  diffusion  pressure 
dependence 

Columns  61-72:  YD  proportionality  constant  for  short  diffusion 

flame  (not  used) 

CARD  12:  Data  Card 

NUMBER  REQUIRED:  One  card  per  oxidizer  type 
FUNCTION:  Specify  number  of  particle  size  distribution  modes 
FORMAT:  (15) 

Columns  1-5:  MODES,  number  of  particle  size  distribution 
modes 
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CARDS  13  — > 12  + MODES:  Data  Card 

NUMBER  REQUIRED:  One  card  per  mode  per  oxidizer  type 
B UNCTION:  Specify  oxidizer  size  distribution  paramete  rs 
FORMAT:  (3FI0.0) 

Columns  1-10:  SIGMA,  standard  deviation  of  oxidizer  size 
distribution  for  a particular  mode 

Columns  11-20:  DBAR,  mean  oxidizer  crystal  size  for  a 
particular  mode,  microns 

Columns  21-30:  ALFAI,  mass  fraction  of  oxidizer  in  a particular 
mode  relative  to  propellant  mass 

CARD  13  + MODES:  Control  Card 


NUMBER  REQUIRED:  One  card  per  case 
FUNCTION:  Program  terminator 
FORMAT:  (15) 

Columns  1-5:  NSTOP,  If  (NSTOP  ,LT,  1)  START  NEXT  CASE, 
if  (NSTOP  . GE,  1)  STOP  EXECUTION. 
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UUUU5H 
UUUU3D 
UUUU3D 
OUUU13 
UUUU^t) 
UUUUb^ 
UUUUb* 
UUUUD3 
UUUUb^ 
UUUU63 

UUUUb ( 

JUUU  f 1 
UUUU  ts 

UUUU 

UUUU  ft) 
UUU V f t 
U U U 1 U 1 

UUUioA 
UUU13* 
U U U 1 3 3 


- 


Program  listing 


L X I LKNAL  UKbANl 


»<L«L  KAP1  « KAP2  • KPF 

« MO A * M | < 

i KAHILI « KAPirU  i 

i t\bt u 

COMMON 

Ml  , 

A ^ i 

Ah  . 

A F FI  « 

A L F A S 1 

1 

dLTAF  * 

BSLK  , 

CIC>H» 

CONI  * 

CSUbP  * 

2 

LAP  , 

UtLUl , 

b^LKU i 

tP  • 

tux , 

3 

vjAMMA  * 

blvl»w  • 

HDfJ  • 

HLP  * 

F.  A P 1 * 

A 

f\PP  « 

PMW  * 

i 1 Q U u * 

P 0 * 1 b N * 

PSI AK  | 

3 

UAp  i 

OFF* 

CjF  ULL  . 

L L « 

OPF  * 

b 

^Ap  t 

Kp  * 

KHl/F  • 

K Ft  0 S P « 

KHOX  * 

COMMON 

sux  * 

1 AP  * 

TfiV  i 

TP  * 

1 ZtFU i 

1 

ALAMO  * 

Xfsil  t 

xhk  • 

XCMUP » 

ANUS  1 ♦ 

2 

AS  I PF  t 

Ab  1 PO  « 

aSTAKU. 

XS  1 AP 

LOfiliUf’J/UUUBLL/ 

P-OA  , 

iviT  , 

IS 

COMMON. 

/ ALL/ 

P 

UOMMON/UP 1 1 N/ 

N « 

Ut-Lb  • 

ULM4N * 

1TL1M. 

COMMON/  X i1*  I / 

K 

COriMON/OAl*/ 

NUA 

common 

/ K WN/ 

bt  1 A , 

KHUF  t 

UM 

common 

/ 1'^P  IK/ 

PLSUbP « 

i PLAfitJ. 

A F * 

1PK1IM  | 

COMMON/bLPl/  1FU(3)«  bNOu(3)*  XNU10l3)t  XNUPU  l 3 ) • 

1 PM#L’l3l*  ULL‘(3)«  HHUXU(5)t  AOXU(3)*  LJXU<3l.  1APU(3)* 

2 ClbNU(3)*  PUwibUl3)»  P0*UU(3).  LUNFu(3)«  IaVu(3>* 

0 KPPu(3)<  r\  API  U ( 3 ) « rvAP^uli).  XMUlilt  XN2U13J*  XN3U(b)» 

4 C S U LJ  P U l 3 > » ALAMbUt3J*  b A MM  AD  ( 3 ) * APHU13)*  tPaL  ( 3 ) ♦ YU13) 

common/kii/  ps*p< iuut3j * xmuoo)»  oiiuo).  pki^ot*  kki2ui* 

1 ♦ IMLUUNI 

UlMLNSION  A(2SI*  ISV(lUU)*  XlllO  0»*  PUllUUJ 
LAlL  1NPU  I It  JJ*  U *U*F  SKP  • nlUUImI  » OUU  * iUA  ) 

CALL  LONLAHU) 

JAK  = U 
Kb AK  = U*U 
*P  = 1 

px  = m • /*p 

wK  A T t t b * bU  U2  ) PXt  xox 
kvKifLIbOUUU 
K = 1 
bO  10  9 
25  CONTlNUt 

wKllLloOUU^)  PX«  10X 
Ip < KP.tb. U i bU  TO  bU 

hKl  I L T b » 3 U U 1 ) 

bU  lon T j rjut 
JAIN  s U 

22  LALL  lIMPu  M JJ  • 2 • 0 * P SkP  • NcOUN  F ♦ UUU  • XUA  I 
lM*P,tb*U)  Is  = ISJlJU) 

V CALL  bDPlAux«JAK) 

TSV(JJ)  = IS 
IP IKP.LN.U / bU  TU  bl 
XSIPUU  = ACjTPU*iUUUU  • 

XSIPPU  = AoTPf-  *1UUUU, 

XSIAPU  = aotah*1UUL0. 

XSfAK  = AO  I AKU*10UUU  # 

MHiltlbiDUUQi  U2LKu « K,  rs*  XuUSl • ALF  AST  * KHUSP«  TP*  bLIAF 
lUt  XSIPPU*  XSIAPU*  XSlA* 

bl  lontxnul 

JAk  = 1 

CALL  XSl  UNI  JJ  * AMT  ♦ ANUS  I • ill  ) 


AO X , 
U3P* 

LPS  * 
KAP2  • 

P S I UP  « 
K * 

KUN 
XALP  A • 
XMU 1 . 


IP  I 


N S A M P 


• XSfPO 


90 


t 


UUU1HU 

1F1NCUUNI  ' JJ)  23.22.22 

UUU  1*0 

Hi 

CALL  INTloI  1 OX. AMI  .O.PbKp.UDU.NCUUNT  .AVAL) 

uuun^ 

2K1 TL l fa • 4UU ) AVAL 

UUUlbU 

XK  = t 1 . U/«HUXU < IOX ) ) *A VAL 

UUU  IbO 

KPiAK  = Ah  ♦ Kb AH 

UUUXbS 

iPlNOX  - Ul*l  90.90.29 

UUUl^U 

29 

1(JX  = IOX  + 1 

UUU1  ! £ 

JJ  = 1 

UUU1/5 

00  10  29 

UUU1  ( b 

90 

JJ  = 1 

UUUi/h 

K = KBAK 

UUUi /D 

KBAK  = 0.0 

U U U 1 / b 

CALL  OUTPJ 1 ( 1 . IUX ) 

uuu^uu 

IF  (P-PSTUP)  31,92.92 

UUU^Ui 

91 

CALL  CUNCAL(l) 

UUU«iU3 

IOX  = 1 

UUU^Ub 

KP  = 1PK1NI 

uuu*u  / 

PX  = 14  , / *P 

UUUiiX* 

00  TO  23 

UUU^lC 

92 

LALL  OUTPUI (2.I0XI 

UUU^IH 

Kt«U(3.109)  NbTOP 

UUU^<f 

blUP 

UUUid^H 

109 

POKMA l(13» 

UUUi't'H 

40U 

POKMA  ( ( 1UA  . 24HVALUL  OP  KATL  INTLbKAL  =.P«.3//) 

UUU^h 

1U01 

FOKMAT  ( 2H0 . 0 ) 

UUU22H 

1U02 

F OKMA 1 ( 313 | 

3U0U 

FOKMAT  (F12.2  .F  11 .4  .F9 . 0 , 3F'10 .4  ,F9. 0 ,P10 .4 . 4H0.2  ( 

uuu zz* 

3001 

FOKMAT  (3X.3HU2LK0.3X.bM  KATL  .4X.4M  13  .4X.bH  A f IUS  T . 4 X . bH  ALP  A b 1 . 4 X 

l.bPi  KH(jbP'4X.4H  TP  ,4X.bH  BL  T AF  . 3 A , 3PIX  b 1 PU  . 3 X , 3HX  S T PF  . 3 A . 3PI A b 1 AP  . 3 
IX. 3H  X3TU • / ) 

UUU**H 

3002 

FOKMAT  ( 12Pt  PKLbSUKL  1 b . F 7 . 1 . 1 0 X . 33H  1 HL  0X1UI2LR  bLINb  CONblULKLU 

lib. 13// ) 

JUU*** 

LNJ 

INPUT 


PKUbK AM  CLNbTPt  1 r J U L. U L)  1 N 0 1/U  BUFFLKb 
UNUbLU  COMP  I LLK  SP  ACL 


414b 
2 TOO 


bUBkOU  1 1NL  INPUT  l JOiMi  J . F'bKP  . NCOUN  I .UUU«  10*  ) 
UUUU11  KLAL  KAP1 .KAP2 .KPF ,MUX , M| . KAPlD . KAP2U , KPPU 


0UUU11 

COMMON 

A1  t 

A2  • 

AF  . 

AF  H t 

ALF  AbT . 

AO  A « 

1 

Ot.T  Ah  « 

BbOK  , 

CION, 

LUf\U  t 

CbUbP. 

LbP  « 

2 

LHP  , 

ULLUi  , 

u2LK0. 

tF  « 

LOX  . 

t_Pb  * 

9 

bAMi^A  i 

GMW  . 

HUN , 

HUP  « 

kAF  1 « 

KAP<f  , 

4 

*PF  t 

PMk  . 

POWU. 

PUWlbN  t 

PbTAKT , 

PblOP. 

3 

« 

OF  p . 

JPULL . 

QL  * 

OPP  . 

K 9 

b 

* A P • 

KP  , 

KHOF  . 

khu^p  « 

khox  . 

KUN 

UUUUXX 

COMMON 

aux  t 

TAP, 

TAV  . 

TF  t 

T2LK0. 

xalf  A « 

1 

ALAHb • 

XN1  , 

XN2  . 

XNUK  t 

XNUbl , 

ANUX  « 

2 

AbTPh  t 

XblPo, 

XSTAKU , 

Xbl  AP 

UUUUXX 

LOI^MUN/UUUOLL/ 

MOX  . 

r-.T . 

Tb 

UUUUXX 

C.Or'ihON/ALL/ 

P 

UUUUXX 

L QIVI 11  UN 

/A^UT/ 

NXCUUn 

UUUUX1 

UOMMON/UAi*/ 

NOX 

UUUU11 

LOi-'MON 

/KWH/ 

BL  1 A. 

KHOM ' 

Un 

UUUUX1 

COMMON 

/ 1 nk  7 k / 

PCbUbP , 

PLAMB. 

XK  * 

1 PK 1 N 1 

UUUUXX 

common 

/bLbF  AL/ 

LIMBLb. 

LRKdLb 
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J J U U 1 1 


CU^MON/bU^l/  IPO(3)*  bl%b(5)  ! XbUlUlb).  X«\jUPU4b), 

i P-^wUlb)*  lJLU  ( 5 ) » KHUAJ(^)*  A (j  X U ( ^ * • LOXU(b).  I APu  ( 3 ) * 

<;  UoMliM  PUwlb(j4b>*  PO*UU(3)»  OUuPUtbJ*  UVUlili 
4 KPT-U43)*  AAPlU(3)«  , XfvlLMbJ*  XluUli)  . XNbU(3)t 

4 CbUbPU43*«  XLAhbU  ( 3 ) t 0ANr'Al)(3)  • A p M U ( 3 ) • tfJoi,(3)  1 Y U 4 3 ) 

UUUU11  LUWMON/F  LAIN/  1UXN  ( 3 ) • iPbc_L 

UUUUll  COnMu(\i/UU  • / ALFAl(b«3>*  olbMl  ( b * 3 ) * UcAr’i  ID»0)  , f.OULb 

UUUUll  LO^^O^/'L/  t-blAKI,  LlAUt  t.lAL.  jA^Ui  UM^L.  O , ^6T*  L1AK,  L I Af  » 

1 IhLAb 

UUUUll  UinLNblUN  AUtbUU).  b ( iOU  ) « NMuUtblb)*  FbKP4iOUt3)*  PuPltoUU)* 

1 X h br\P  ( bUU  i 3 ) 

u u u u 1 1 i M n - l)  it  l « lb 

uuuuli  199  HjKttAl  (ibf^Flu.bf^ib.HQ.bi 

UUUUli  1 Kt  MU  ( 5 , i.9?  J ,JUX»  KblAKl  , KbT  OF  « IPKlNT,  Li^dt^t  tKhbtb 

UUUU33  KLAU(b»luUl  IVCOUM  • NXLQu1**  X.g*  bil  A^Tt  ubIOP*  fSTAKT 

UUUUbH  XNLUUU  = •’*LOUIUl 

UJUUbD  ArjXLUU  =■  Wxcuuw 

UUUUbb  UUO  = tUbtUH/UbTAK  | ) **l  l./ANC-UUi4> 

UUUU/4  XUuO  = 4L3  t UP/US7  AKT  ) **  ( 1 , /XImXCCU  ) 

UUUiU  i KLAU(t)*^.UU/  If-  ULL 

0 IPUtL  = 1 / MlPb  /»  2 /PbAU  / 

Ul'UlU/  Kt  AU(b.UU<i  ) l^tKO*  A MLh  M « Uf-ULl*  KHuh  « Ah*  tF 

UUUliiJ  KLAU{t),«U<i)  Li  L 1 A « KhUHt  jM,  POSUbP*  PLMFltJ  * AK 

UUUltJi  OALL  OUTPU  I ( 4 , IUX  ) 

uuuibu  uo  au  iux-i*uux 

UUUibD  KtAUtb.^UU)  1UX1U 

L 2(jU  FQKrtAT  (ah; 

U U U 1 / i ^UU  POKI1A  T ( lb  I 

C 1MI0XIU.LU.1AP)  00  [0  HU 

0 ih  ( 10X1U.LU.  iKP  ),  00  1 <j  4 £ 

0 IP  4 1UX1U.C.U.  innx  ) bO  lu  4b 

UUU1//  IF4IOXlU.LU.il  bo  T 0 4b 

UUU^Ul  IPUOXIU.LU.^)  bO  10  4 b 

UUU^Uc  IP  4 1UAJ,U»LU,3)  bo  ib  4b 

UbU*U4  4(1  iUXN(iUX)  = i 

UUU*U/  bO  10  bU 

UUU^U/  4,5  1 0 A iM  ( iUA  ) - ^ 

UUUi^lC  00  TO  DU 

UUUJM<£  4b  10AIM410X)  = 3 

UUUUli  bO  L 0*  ‘ 1 1 IMUL 

UUU^lb  KLAU(bittU^)  ll-u(lOX)*  b>.*U(IOX)*  AuuiudO*!  » Ar*UP0ll0X) 

1 . Pr-IWUCIUA) 

UUU^HU  Kf.AJ  ( b . ) ULU(IUX).  KHuXU(IuA)*  AOAUilt'X).  Lb  XD  4 I 0 X > » T APO  4 1 0 X ) 

UUU^bb  KLAU(b,ttU^>  L10NU410X)*  PUw  1 bu  ( 1 uX  ) * POWUU t 1 0 A I , CUNFUllUX) 

U4JU311  KtAj(5»ttU^)  l\PFU(10X)*  KaPIlK  iOX)  • FmP^U ( 1 U A ) « XM1U41UXI* 

1 XN4U  < 10A  I * A u ^ ( j ( XUX  ) 

UUU34*  KL  AU  ( b * ttU<£  ) ObUbPu(lUX)*  A i_  A -ibb  4 1 b A ) . o AMI*  Atj  ( 1 U X ) * ApnjllOX)* 

1 LPOUUUA/*  T 0 4 1 0 A ) 

UUU3//  3 TAUU(iOA)  = 4TAMU410X)  ♦ IPJ(  l Ox)  )/<i»UL)U 

UUUHlb  KLAU(b*19^)  noULb 

U U U 4 <sl  4 IJ^iUULb  ( 1 JA  < = MOUtb 

UUU4 dt  UO  1U  II  = lfMUULb 

UUU4b4  KLAU4b.  lUi)  bittMAl ( 1 I * IUA ) • Ub A K 1 4 1 1 • 1 0 X ) * ALF A J ( 1 I t IUX J 

U'jU4t>/  10  CONTINUL 

UUU4bl  i-IOULS  = NHUDLb(IOX) 

U 0 U 4 b 4 call  UUTPU I ( 3* lux ) 

UUUH/*  HU  LON  r 1 NUL 

UUU4/D  JJ  = 1 


92 


UUU4 

bun  = u.u 

UUUH/b 

Hbun  = u»u 

U U U / f 

bun*  = u#u 

UUU4  It 

c 

II 

C 

C 

uuuduu 

UO  3 IUX1  =l.NOX 

UUU3U<J 

MOULS  = NMUULSt  IOX1  ) 

UUUSU4 

UU  4 IJ  = i.MUULS 

UUUb^l 

XSUM  = ( U »U/KHOXU( IUX1 ) ) - ( 1 .0/KHOF ) )*ALF A1 t 10. 10X1 ) 

UUUD<d^ 

bun  = SUM  ♦ XSUM 

UUUtJ^b 

4 

CONTINUL 

UUUS24 

SUM2  = SUM2  + SUM 

UUUtJ^b 

bun  = O.u 

UUUDb^ 

S 

CONTINUL 

uuuabD 

KHUT  = l.U/l ll.Q/KHOFI  + SUM2) 

UUUD4U 

WKlTfc.tb.4UU)  KHOT 

UUU546 

400 

FOKMAT ( lhi . 1UX .2UMPK0PtLLAf\IT  QENS1 1 T =.F7.4/> 

UUU5t)C 

SUM  = U.U 

SUM2  = U.U 

uuussa 

UO  7 10X2  = l.NOX 

UUUDDD 

MOULb  = NMUDLSt  10X2) 

ouut>b  t 

UO  b 1JJ  = UM0ULS 

uuub  t $ 

XSUM  = t ALfc A1 ( 1 JJ. 10X2 ) /KHOXDl 10X2) ) 

uuub  ^4 

SUM  = SUM  + XbUM 

UUUb  /D 

b 

COM  T INUt 

UUUb  l b 

SUM2  = SUM*  + SUM 

UUUbUU 

SUM  = O.U 

UUU6U4 

7 

CONTINUL 

UUUbU  t 

XNU  = KHUI*SUM2 

UUUbll 

WKlTfc.tb.4Ui)  XNU 

UUUblb 

401 

FOKMAT  1 1 U x . 21HV0LUML  F KAC  OF  OXIU  =iF7.4/l 

UUUb^f^ 

FSUMKP  = u.O 

UUUb** 

SUM 4 = U.u 

UUUblb 

UI  = USTAKI*XUUO 

UUUb*D 

UO  14  IIl=i,NXCOUN 

UUUb*  ( 

UO  12  1UX0'=  l.NOX 

UUUbJU 

MOUfc.S  = NMUDLSt  10X4) 

UUUbb* 

UO  11  JUJ  - 1 • MOUt S 

UUUbbH 

CALL  UISTF  (Ui,  UUJ. 11 1 . 10X3 . YVLOK.UbAKl .blUMAl ) 

UUUb4D 

XF  SUM  = ALF  Ai  ( JUUU0X4) /YVLCH 

U U Ub4b 

FSUMKK  = FSUMKF1  + XF  SUM 

UUUbbH 

11 

CONTINUL 

UUUbbl 

XFSKFM111U0X3)  = FSUMKP 

UUUbb* 

FSUMKF1  = U.O 

UUUbbb 

XSUM4  = tKHOT/KHOXUt  IUX3)  ) *XFSKPU1I  • IUX4) 

JUUbbb 

SUM4  = SUM4  + XSUM4 

UUUb  t $ 

12 

CONTINUL 

UUUblb 

FUPl 11 I ) = SUMS 

UUUb  f b 

SUMS  = U.U 

UUUBUU 

XU ( 1 I 1 ) = Ui 

uuu^ui 

ui  = ui*xouo 

uuu  nj 4 

14 

CON  t INUt 

UUU  \ Ob 

CALL  CLALtFUP,XU,X(j.XNU,xUUO.t) 

uuu  tit. 

WK1 TL ( b . 402 ) 0 

UUU  ( *b 

4U2 

FOKilAT  t 1UX.23MC  FKUM  VF  = C*U2tK0**N  1S.FB.S/) 

uuu  rdf 

UI  = USIAKUUUU 

UUU  f bU 

UO  24  I IL  =1.NOOUNI 

UUU  M* 

UO  22  10X4=1, uUX 

UUUMi 

FbUMKP  = U.O 

uuu/io 

MOULS  = NROOtbl  10X4  ) 

UUU  f $t> 

UO  21  022=i,m0ULS 

UUU  /4U 

CALL  UXb IF (01 .022  <122. 10X9 .TVLCR.UbARi. S1GMA1 ) 

U U U 7 d 1 

XFbUM  = MLFftl (022> 10XAI/TVLCR 

UUU  / 

FSUMKR  = FbUMKR  ♦ XFbUM 

UUU  / bU 

21 

CONT 1NUL 

UUU  fbC 

FbRRl 122. 10X4 ) = FbUWAR 

uuu  r b f 

22 

CONTINUL 

UUU  f f 1 

U( 1221  = Ui 

UUU  f ft 

U1  - UI*UUU 

U U U 7 t *♦ 

23 

CONTINUL 

UUU 

XOX  = 1 

uuu  f t / 

lb 

CONTINUL 

UU1UU3 

U2LK0  = U(OU) 

uuiuu/ 

XNUbt  = l.U/ll.O  + b.  U*C*L)2E:ko**  < XN-3.  U ) /3. 14ib92bb4  ) 

UUlUlD 

ALF  AST  = 1.0/IX.0  + b . U*L*RHOF  *D2tR0*»  ( XN-3 . 0 ) / 1 3 . lRlb92bbH*KH0XD  l 

110X ) ) ) 

UUlUrff  / 

KHUSR  = KdUXUUOX)*XNObT/ALFAST 

UulU.il 

3b 

CONI 1NUL 

UUlU^l 

ML  1 URN 

UUluJ^ 

1U0 

FORMAT  (2XO>5H0.5) 

UUlu4«^ 

1U1 

FORMAT (3F1U.3) 

Uuiu^d 

OU2 

FORMAT  (bLiii.b) 

Hb* 

UUiUi^ 

C 

LNU 

*«***»*»  CCAL 

******** 

U22i 


bUbRKOGRA M LLNGTFI  11UA3 

UNUbLU  COMRILLR  SRACL  aUO 


bUbHOUJ 1NL  CUALlFUR.AU.XN'XNu.XMULI  «U  ) 


UUUUXU 

COnMON  / X m u T / NXCOUN 

UUUU1U 

UXMLNblON  FUR(bOO>«  XUltfOC) 

UUUUlU 

A = O.U 

UUUUXU 

b 1 = 3.141092bb4*t 1.0  - xNU ) / b • 0 

UUUUI* 

02  = 1 

UUUUX  ( 

2 

XL1  = FUR(O/)*XU(02)**lXN  - 9,0) 

uuuu^-i 

XL2  = FURtO*  ♦ 1>»XD102  + 1)**(XN  - 

UUUU.il 

XOUO  = XUTU2  + 1I  - XO(u2) 

UUUU.it) 

ARART  = l XL 1 + XL2)*XUUU/2.0 

UUUU4U 

02  = 0^  ♦ 1 

UUUU^l 

A = A + ARART 

uuuu^a 

IF  (02  - NACUUN)  2,  A,  A 

UUUU<+3 

4 

C = bi/A 

uuuuh  f 

RL 1 URN 

UUUU4  ( 

LNU 

c 

**•»*«**  UISTF 

bUBHKUGNAf 

LLNGTH  102 

- 4.U) 


UNUbLU  CUMRlLt-K  SRAUL 


a3uo 


SUbKOUT  lilt  OiSTF  (Ui  • JJJ«  ill  ,XOX.TVLCR.UbAKl.SiGMAX) 


UUUUXl 

UJMLNS10N  UBAK X ( b , 3 ) , bluMAIIb 

UUUU11 

X - 

ALOG 1 Ul ) 

UUUU12 

XM 

= ALUO(UBAKKUUU.IUX)  ) 

UUUU2  / 

Sib 

= ALUolSiUMAl 1 UOO ,10X1 ) 

U 0 U U 3 A 

L\y 

= ,b*> (X-XM)/blG)»»2 

94 


I — ™ 


0UUU4U  if-  t ZiP.b  I • IbU.  ) *1*  - XbU« 

OUUUbo  Yv/LLK  = Sib**?  . SObb*  to*f  / Su  *L  xP  ( 2 XH  ) 

UUUUbU  *L lUKN 

UUUUbU  tNU 

C *••••*  C Ui,L  Al 


SUtfPPUGK  Art  Lt'^blH  A.'* 

UNU6LU  CU^KILLH  SHALL  HbUO 


bUOMOul  lot.  CUKLAL 

( J 1 

4fe* 

L 

rMli>  iUbHUUtJUL  lfcLKLI’lt.lJli* 

THL  FKLbbUKt 

46* 

U u U U U £ 

umtNsiuN  »*  al  ( i u ) 

4b* 

uu  u u ut 

KLAL  KAP1*  -SAP 2%  KHh  • I'lUA* 

Mr 

U U U U U £ 

LOMMON  ax. 

A if  * 

*«f  • 

Ah  M • 

ALh  AST* 

AUX  . 

1 ULltth , 

bswr< , 

C.  I &IX « 

LUNA  » 

CSUPPf 

L3P  • 

^ L4p , 

ULLJ1 , 

I 

Lh  . 

LOX  , 

LPS# 

$ o«rnA< 

bP'i*  » 

MUIn 

HUP  * 

♦v  API  * 

KAP2  * 

4 *Hh  , 

PM*  • 

HOWlm 

Potion* 

PSFak  T * 

PS1 UP* 

b «» A H « 

Oh  b t 

jf-Ut-L. 

ul  * 

l»Hf  t 

K « 

b ft  AH  f 

Kh  t 

khuf  * 

KMOSP  » 

khux  * 

KUN 

UUUUU<! 

COMMON  WXi 

1 AP  t 

U>\l , 

lb  * 

1 ZLKU, 

X ALh  A * 

1 ALAMO* 

Xf4i  * 

xrji . 

Xi^UP  • 

XNUS 1 « 

XNUl « 

2 A b TPH  * 

XilPu, 

xSTAKu. 

XS1  AH 

UUUUUC 

LOP''Murj/UUubLL/ 

MO  A ♦ 

hT  i 

1 s 

U U U U U C 

LOMMON/ ALL/ 

P 

U U U u u t 

1 

I F < J > 2%  2>  lb 

46* 

UUUUUH 

2 

PAL(1)S1*U 

UUUUUb 

P AL l 2 I =1 . t t 

uuuuu  / 

KflL(4l=iti/ 

UUUUiU 

h AL ( 4 ) =S • 02 

uuuuu 

b C.L  t b > =iU • U 

UUUUii 

JJ  = U 

Hb^ 

UUUU  1* 

1 = 1 

46* 

uuuu  i a 

lb 

If  (JJ  - 11  lb,  20 

* 2U 

•+6* 

uuuutu 

lb 

XMuLF  = PSIAKT 

4b* 

1 J U U U 2 i 

JJ  - 1 

4b* 

U UU  Ud2 

1 = 1 

4b* 

UUUU23 

2D 

lb  (i-6)  Ou, 21*21 

UUUUtb 

21 

XMUL  1 = A •*> U L 1 ♦ 10  • 

4b* 

UUUUBU 

l = 2 

46* 

UUUU-31 

bO 

P = XMULI*HAL<1) 

46* 

UUUUii 

1 = 1 + 1 

4b* 

U U U U 0 b 

«L 1 UKN 

4b* 

JUUUiS 

LNU 

46* 

L 

******** 

tjUP 

******* 

* 

SubPKUGHAl^  LLNblH  Qb 

UNUSLU  CUi'iPlLLK  SPALL  4<fU0 


U U A U 
UUHl 


UUS  f 


UUbO 

UUb* 

UU/U 

uu  f 1 
uO  ft 

uu  f s 

uu  o 
uu  fb 

UU  f t 

uu  fb 
uu  /* 

UUttU 


SubKQbllUL  BUP<lO**JAK» 

UUUUU4  HLML  K API *KAP* *KPh  *r,'UAt'-!rfKAPlU*PAPiUtKPf-U 


COMMON 

Al  , 

A 2% 

rtF  . 

Ah  H* 

ALFftSI  . 

AUX* 

1 

OLTAh • 

bSUK  , 

C 1 G to  * 

L ON  1 • 

CSUbF. 

LAP  • 

2 

LHP  % 

ULLui  • 

UZt*0« 

Lh  • 

LOX  , 

LPS* 

S 

V»AMMA  • 

Cm*  ♦ 

HUN* 

HUP  • 

KttFi  « 

MAP* » 

4 

rvpf-  • 

PM*  • 

P 0 w u * 

POWibN* 

HS1 ANT. 

PSIOP* 
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r 


b ^*P* 

OFF  . 

jFUfc-L* 

UL  * 

UPh  « 

b kmr, 

NF  . 

khof  • 

KH03P, 

KHOX  t 

UUUUUH 

COMMON  5UX» 

T AH  i 

TAV* 

Th  • 

T^LKU  t 

1 XLAHb* 

XIM1  * 

xr*2  * 

XNUP  t 

XNUb  1 « 

2 xbTP)-  « 

xstpq, 

aSTAku. 

XbT  AP 

UUUUUH 

LOHMUN/UObdLL/ 

MOX  * 

MT  . 

Tb 

UUUUUH 

COMMON/ ALL/ 

p 

UUUUUH 

COMMON  /APUT/ 

NXLUUl* 

UUUUUH 

LOMMON  /Hfti'k/ 

bLl  « < 

hFtOM  • 

UM 

UUUUUH 

LOnMON/bUPl/  T F U ( 3 ) 

. oriWLMi) 

. XMjIu ( i ) 

« XiJUr'U  l 3 ) 

« 

1 PM*Ul3)f  «l_U(3)t  KH0*U(3>*  AqXUO)*  LOXUCOIt  IAPL)(3)t 
Z C1GNIM3M  PU*lGLM3)t  F0wUU(3)*  L0NPLK3)*  IAVU13)« 

3 KHFLM3)*  *AP1U(3)»  KAP2U13),  XI\tlUldl«  XN2UI3)*  XN30(3)* 

H CbUbPUt3)»  XLAWtfUl3>»  GaMMAI) ( 3 ) » Af-MUl3)»  LPoU(3>«  YU<3) 


UUUUUH 
UUUUUD 
UUUUUb 
UUUUl  A 
UUUU^U 
UUUUUb 


IMJA*)  1*1,2 
Tb  = VUU. 

A1  = • b* l 1 • 0 ♦ 1,0/bWKI 13.0) ) 
A2  = U.b*»l.U  - 1.0/bUHT (3«0  ) ) 
K = ,Ub*£>UKTlP) 
box  = xrjuoi 


u u u u 2 t 
uuuuou 
UUUU3* 


OM*  = bMKU l I UX ) 
XlMUl  = XUUIO(IOX) 
XNUP  = XNUPuliOX) 


UUUU33 
UUUU3b 
U U 0 U 3b 


PM  * = PM*U  \ IOX  ) 

OL  = ULUlA^X) 

NhUX  = HhUXUllOX) 


UUUUHU 
UUUUH1 
UUUu‘+3 
UUUUHH 
UUUUHb 
UUUU**  ( 


AO  X = AOXUIIOXI 
tox  = tOXUUOX) 

TAP  = TAPUUUX) 

LION  = C !<»!<□  I IOX  ) 
POWlbN  = PUWibUUOX) 
PO*L>  = POWL'DI  IOX  ) 


uuuubi 

uuuud* 


CONF  = CUNhOllOX) 
KPr  = KPPU'IOX) 


UUUUbH 


KAMI  = KAKlO(iOX) 


UUUUbb 
UUUUD  ✓ 
UUUUbU 
UUUUbH 
U U U U 63 
UUUUbb 
UUUUbb 
UUUU  / U 

uuuu  r l 

UU  UU  /3 

UUUU  t 3 
UUUU  / H 
UUUU  f / 

uuuxuu 

OUUlUl 


KAP2  = KAP2DULX) 

XN1  = XNIOIIOX) 

XN2  = XN2UII0X) 

XN3  = XN3U 1 1UX ) 

CSUbP  = LDUBPUl IOX) 

XLAMB  = ALAMBUlIOX) 

GAMMA  = bAHMAU(IOX) 

AFH  = fib hU l I OX  ) 

LPb  = tPDU»  IUX) 

2 CONTINUL 

L LALL  hLAML  I ( ALFAS1  * IOX • TF ) 

TF  r TFUdUX) 

TAV  = C f P ♦ 1 AP ) /2 . U 
LALL  STLMP 
HE  I URN 
LNl) 

L ******  sixth 


bUbPKUGR AM  LLWGFH  133 

UNUbLU  COMPiLtK  SPALL  3bU0 


SUBROUTINE.  STLMP 


K t 
RUN 
XALH  A i 
XNU1  » 
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UUUUUi 

Nt*L  KAH1*  KAH2* 

kpp  . nox . 

MT 

U 0 U U O 1 

LUHMON 

A1  « 

A2  t 

AF  . 

AH  H • 

ALF Ab 1 . 

AUX  • 

1 

DL  T Ah  • 

bbon . 

CION. 

CONI. 

Cboop* 

LAH. 

d 

C*P  , 

OLLOl , 

0/t-KO  * 

tH  . 

tox . 

tHb. 

5 

UAMHrt  « 

0M1«.  . 

HON* 

HUH* 

K API  . 

*AH*  . 

* 

. 

PMA  . 

P0F.O* 

HUWlbN. 

PbT AKT . 

Hb 1 UH. 

s 

WMH  t 

OFF  . 

wF  OtL • 

OL  * 

OPF  . 

rt  • 

b 

NMH  , 

KF  • 

KFlOF  . 

KHUbH ♦ 

KHOX. 

KUN 

UUUUUI 

LOnrtON 

bUX  • 

T AP  • 

T AV  « 

Th  • 

T 2t  KU < 

x alh a • 

1 

*UAHb • 

Xul  • 

Xfj2  < 

XNur . 

XNUb 1 • 

XNU1  . 

AbT  • 

XbIPu. 

X ST  AKu . 

XJ>  1 Ap 

UUUUUI 

LO^'MON/UUUbLt./ 

MOX  • 

MT  i 

Tb 

UUUUUI 
UUUUUi 
u UUUUl 

uuuuu/ 

uouulu 

UUUUl 1 
UUUUl^ 

uuuulo 
uuuu  1* 


UUUUJJ 

UUUU-3^ 

UUL'Uio 

uuuu*u 
u u u u * c 

uuuu** 

UUUUHH 

OUUU*D 

uuuu*  / 
uuuu*  f 

UUUUOl 
UUUUDO 
UU'JUDt) 
Uuuubu 
JUUUb^ 
UUUUb* 
U U U U b f 
UUUUl* 

uuuuf i 

UUU 1 u* 
UUUlUO 

uuu l u* 

UUU1 1* 
UUUlib 
UUU1 1 i 
UUU 1*U 
UUUi^C 

U U U 1 1 D 
UUU  l*b 
UUU 1 *b 
uuuuu 
UUUlib 


L 


T2tK0)  - AlF  A*  I I AP  - T2tK0)  + ll.U  - ALFA)/ 


LOrtilON/ALL/  P 

LOllflON/t  A/  KP  • KP»  tTAF,  tTAP 
UIl-itNSION  WHILST  I 20  > . «U I S ) 

ALMA  = MLI-Hbl 

xrjo  = XNObi 

KHUP  = KHW3H 
KP  = U 

KP  = U 

UFP  = LbOOP*|(TF  - T2LK0)  - AlFA*UAP  - 1 

1 C AUBP*UF ULU t /ALF A 

UPP  = LbODP»lltP  - T2tK(J)  ♦ ALF A»UL/tbUbP 
1 LbOOP) 

UAP  = CbOOP»i|AP  - T2tK0)  ♦ UL 
XfJL  = XN2 
KAIL  = FAPi 
i'll  = HFiOP»« 

*ni  = u . u 

btOXONXNO  UF  LOMPtTXNO  FLAME.  CALLULA  I ION 
AKoLb 1 ( 1 ) = T b 
AKOLb)  t 2 > = K A T L 

AKOLSHAI  = Xl'JL 

I [I  Lk  = u 

LONVtKOt.Ul.l-  LaLLOLaTXON  on  Tb  FOLLOwb 

II  — IS 

t i = DUN  l A*|,Lb1  ) 

I ? = I A - tl 

12  = AF A A 1 l T2 < T2LKU ) 

12  = APXNA I T2 . TAP » 

AKoLb 1 II)  = 12 

L 2 = UOU I ««OLb I I 

I A - 1 1 1 * u - t2* 1 1 ) / t tX  - t2) 

| A = Af.BUI  T A , I 2LKU  ) 
lb  = A I-  11.11  1 A.  TAPI 
AKOLb I I X I = I A 
tA  = UOUlAKOLbl) 

IPIAobtl.  “ (tA+lb)/lb)  - ,L0])  blKtit’iul 

XF l Aribtt2 I .L I . AUbltl ) ) uu  fu  70 
L2  = tA 
12  = I A 

XT  ILK  = 1 M tK  ♦ 1 
IF(lTltK.oi.S(j)  00  TO  OU 

00  lO  ID 
tl  = tA 
IJ  = IS 

1 1 ItK  = 1 I ItK  + 1 
iFlllltK.ol.su)  00  TO  OU 


Abb  U097 


- ALFA)*OFUtL/ 


Abb  U X U 2 
Abb  U1UA 


Abb  U1US 
Abb  U1UO 


Aob  U 1 X A 


. U 0 1 > bU*oU*ul 
To  71) 


00  TO  OU 


00  TO  OU 
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UUU 1 OS 
UUUl ss 

bO 

bO  lb  HU 
M H i T L 1 1>  « ) 

UUU l a 1 
UUU1H1 
UUUl A 1 

bU 

► UNMAl  ( SX  t • 1 1 1 tHA 1 A OU  Liril  1 

► UMi'iAl  l >LlO,6,Hu.<r) 

LON  1 lNUt 

tXCLLULU* I 

UUU 1 A 1 
UUUiHS 

w = MT/KHUH 
KL  1 UKIM 

AbS 

UUUl'FA 

L 

LNU 

****** 

DON 

AbS 

****** 

SbbMKUbNAri  LtNGIM 
UNUSLU  LOMHILLN  SHALL 


FUNCTION  UUN  (AKbLSl) 

LALLULAlt  IS  F UK  LUPlHLU^b  FLAKES 
KLAL  KAMI*  *Ap2,  KHF  , nux*  MT 


Hb*  U2*H 


uuu  uu^ 

COMMON 

AX  , 

Ai>, 

At-  . 

AF  H • 

ALF  AST  * 

AUX« 

1 

bLTAF  , 

bbu«  , 

CIGI*. 

CUM  . 

CbUUH, 

L^H, 

£ 

L AH  • 

L'LLuI  . 

jZL«0* 

LF  * 

LUX « 

L HS  t 

A 

b aMM a • 

bMl*  * 

mUN  • 

HUH  * 

K At  1 • 

KAHi?  * 

A 

• 

He'll*  • 

HOt*U  * 

HUMlUM* 

HS1 AKT . 

HS 1 UH. 

D 

WAP  t 

UFt  • 

y F ULL • 

OL  * 

UHF  • 

K • 

o 

KAH  • 

KF  , 

khgf  • 

KHO^H , 

FthUX. 

KUl^ 

uuuuu* 

COMMON 

iux  • 

T AH  • 

TAV* 

T h • 

F XLHU • 

XALF  A • 

1 

XLAMb  * 

XNi  • 

AhU  • 

XljUH. 

XNUb t • 

Aigul « 

XOTHF  » 

XbIHU* 

XSTAKU* 

Xil«H 

uuuuu* 

LOHMUN/UUUti  LL/ 

hox  • 

FIT. 

lb 

uuu  uu^ 

COMMON/ ALL/ 

H 

uuuuud 

COMMON 

/K#N/ 

at  1 A. 

KHGF'  * 

UM 

uuuuu* 

UOnMON/t/  c-SIhHI, 

t 1 au * t.T«L 

< oselUt 

JA^Lt  0 * Xbl 

UUUUO* 

uuuuud 

uuuuuc 

UUUUU3 

UUUUUA 

UUUUUb 
UUUUU  / 
UUUUl 1 
UUUU13 
JUUVdL 
UUllU^i 
UUUU*r> 
UUUU SI 
UUUUHU 
UUUUHU 
UUUUHD 

uuuus* 

UUUUSb 

UUUUbS 

UUOUbD 

uuuu/u 
UUUU  t D 
uuu 1UU 
UUU1US 
UUUlUo 
UUUlll 
uuu  i lb 


LOMWON/t^/  KF  « KH*  LTaI-»  lJAH 
DIMENSION  AKoLSKifU) 

T S = AKbLa 1(1/ 

KAIL  = AKOLS 1 ( d ) 

XNL  = AKbLbT  V 6 ) 

HHUH  r KMU^P 
ALF  A = aLFaSI 
XMU  = XNU^I 

MO*  = AU^»LXK ( -LUX/ ( 1 .**tt 7* 1 S) ) 

LALL  SOXLAl  ( I S*  XN'JtSuA  ) 

MT  = ttUA*oux/ALFA 

IFIKF.tb.U/  L 1 AH  s LSIAKT 

XSIHO  = L I MCALl XNU« 1 tL I AH ) *b/^Sl*AHh 

KH  = 1 

X S I HI-  = M I / ( KHF*H**XN1  / 

XSIAH  = r.UX/ ( KATC*H**Xl\iC  ) 

U Uh.tUiUI  LIAF  s LSIAKT 

XSIANU  = B*t 1 ALALI XNU*U,lTAF )*AFFi/*S1 

KF  = 1 

titlAI  =(X^IAH  - XS  IFF  / *AHH/XSIHL 
IF  ( bcT AF  * 1 • U ) bLTAF  = 1,0 
IF (dLTAF  .ul.U.U ) bL  T AF  = U.O 

zah  - csunH*nux*xs r ah/xlmho 

/Al  = CSUbH*rl  | * ( XS  I AH  ♦ XS I AND ) /XLAMd 
IF  (bLTAF -1.0  I 97t*L«yu 

/HF  = LSUOH*nl /XLAKib*(  ASTHF  ♦ XSIHU/ 

bO  TO 


H69  U241 
Hb*  U £*£ 
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UOUliD 

it 

ZHF  - LdUbH*nl /XLAMb*  l ♦ XblHF) 

UUU  1*.'M 

Vti 

IF t^^l.Ll.-lUU.I 

^Al  = - i u u , 

UUU1</  / 

1FIZAH ,LI«”1UU.I 

A AH  = -iUU. 

UUU1SC 

If-  1 /P^  .L  1 --lUO.  f 

ZHF  = -J.OU. 

UUU IDS 

11-  Unl.bl  tlUU.  1 

Z A 1 = 100. 

UuU  1*4  U 

it  (4  At'.  1>I  .ICO.  > 

^ A P = 1 Jo  • 

UUU  1*0 

if-  i ^Pf-  .0  1 • 10U  . ) 

ZHF  = 1U0. 

UUUl*4h 

ATS  = IZL*U  - ALf  Ahul/LMJHH  - (1.  - ALI-A  - bL  1 A ) *v.t-  L'LL/C  dUHH 

1 - bL 1 A*nn/CdUbH 

♦ 11.  - bLI AF > * l ALF A*UAH/LdUbl *E AH l -Z AH » ♦ UK  t / 

?.  CiJbH*LXH  l -ZA1  ) *ALI  A 1 ♦ dLT«F  *01  F ZLdUbH*LAH  ( -*;hf  j 

UUUciU 

XTS  = API  A a i IMS 

. IZLnUI 

UUU^ID 

DUN  = AKbLSTIl) 

- XTb 

Hb^ 

U<t?b 

UUU^IH 

bO  T U 1UU 

UUU^i 3 

't'i 

LION  s AKbLdTll) 

UUU^ib 

1UU 

HL 1 UKN 

Hod 

\J<iV  ( 

UUU**U 

LNU 

•♦to* 

Oc !*0 

L 

****** 

SUXCaL  **•*» 

* 

SUttPHUGHAM  LLNbfH  <>/<> 

UNUSLU  CUn^lLt-K  SP  ACL  o**UU 


bOoKUUl  X ML  SUXCAL  U bl  tXiMU*  bOXi  ) 
UUUUUD  «t«L  «AHH  KAP2*  KPP  « MUA*  MI 


UUUUUD 

COMMON 

*lt 

A*;  « 

«F. 

Af-Mt 

ALF  Ad  r • 

AUX  . 

1 

ULTAf-  , 

bSUH  # 

C 1 bl.  « 

tot'll  . 

LbUbH, 

LiH, 

2 

f*p, 

ULLUi  t 

LlZt  HU  t 

LF  . 

LUX. 

l.Hd  . 

5 

bA%HA  « 

bhy  • 

Hlll.t 

HijH  t 

K AH1  . 

AAHZ. 

H 

HPf-  • 

PrtW  « 

t OwU  t 

HOwibNt 

Hbl A«  I . 

HblUH. 

D 

y ap  « 

WP  P • 

jMjLL* 

xL  • 

I.HF  . 

t<  • 

o 

HAP  • 

HP  « 

(,H0F  » 

KFlOdH  , 

NMOXt 

KUH 

UJUUU3 

COMMON 

DUX  t 

1 AP» 

TAV  t 

IF  . 

TZLHU. 

xalf a. 

i 

ALAMO  « 

Al\lX  • 

\Nl  • 

Xi  jUH  t 

XMJb  1 • 

XIMUl  . 

2 

Ab 1 Ph • 

AS  1 PU« 

» ST  A AO « 

XS  1 AH 

UUUUUD 

COP’ihON/UOUbLL/ 

hOA  • 

fll  . 

lb 

UUUUUD  COMMON/ ALL/  P 

UU'JCUS  Kh  = AP*LAP<-tP/<l.VO/*lbl)  )/KHOK 

UUUU<fU  HAP  = MJA/HMUX 

UUUU*<i  T1»H  = c lOfj*  \ uZLHU*U  . UUUi  ) **  ( POfcU+1  . I /p**PU*Ibb 

UUUUDO  XIJU  = Hh*  I ibn/ iu^lkuau  «o  JU1  ) 

U'JUCSD  VUU  = XUU  ♦ Kh/KAP 


UUOUJ/  UUC  = AMINA  l 1 • U • VUU  ) 

UUUUAO  KAPQKJ-  = n«p/Kp 

UUUU*0  bUAi  = f).U*ANU*  ((  l A • U-KAf  UHF  ) **2-1  • u I • l UUO**D-XPO*»3 ) /$.U  ♦ 

1 l 2 . U^XCU*PAPUKf  • ( 1 . U-Krt»*C'Hh  ) ♦ X . U I * l UU0*OUU-XUU*Xl>0  ) /2  • 0 ♦ 

2 l KAPO*»-  *A«JU  > *♦**  l buU-XuU  ) ) 

UUUU/U  If-  ( OL'X  i . L • • 0 . U ) SUX1  - . j 

UUUU/*  hLIUKN 

UUUU/D  LNU 


bUBPNUCjH  AM  LH^bVH  ID/ 

UNUSLU  CU^PILLh  bPACL  4UUU 
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*•**•*•*•**  FLaMlT 

bUbKOUT  INC  FLAMt T I MLF  rtb  r f XOX  • TF  ) 
uuuuua  mnn0K/FLMn/iuxi\i(3t « ikulL 

uuuuub  uintubluu  IFA12«.S«2X  I « ti  I PoAf(  2i  I ' F’tXAUAPt2XI 

UUUUUS  IF  = uuu 


UUUUUo  X = XUXIJtXUA) 

UU0U1U  liAIAI  UTPO«P(  J|  . 0= 1 • 21 ) / *bV . 0 , lU2b.U,  1UBO.U.  114B.U.  12Ub.U. 

1 X2b«.U,  liUO.u,  libb.u,  X^y^.U*  X44U.U,  X47  f . 0 » lbXX.U.  lbH2.0« 

£ 13'S.Ui  IbbU.Ui  xe2b.U,  tbUu.O.  FUb^.O.  A2XU.0,  2bbU.U.  1HUU.0/. 

A l PHAIJAPl  JJI  ,Jj=l,21l/lUb9.U.  lUb9.U.  X1U1.U.  lliO.U,  lXb/.U. 

A XIOA.U.  X2UH.U,  X2S2.U,  X2bb . 0 « 120U.U.  XbUH.O.  IbiU.U.  lib/.O, 

b XbbO.U.  Abby.u,  2U2o.U,  2SXX.G.  2939.  0,  299U.0.  2**0i.U.  132U.U/ 

uuuuiu  uo  iu  j=x*2x 

UUUU2U  IFAIi.X.O*  = hthhahiji 

UUUU2X  TFA(2«X.0I  = KbAI.AP(J) 

UUUU22  XU  L0NT1NUL 

UUUUJU  AKA  = (ALKAbl  ♦ ,Ub)*2U.U 

UOUU32  AK  = XF  K ♦ 1 

UUUU-lb  JJ  = IKK-iltO 

UUUUJI  A JO  = JJ/XuO.U 

UUUU4X  IF  = I F A l irgt.L  > 1 ,KA)  - 2 0 . U * ( * J J- ALF  Ah  ( ) * II  F A l 1 F UEL  • I • KK  ) - 

X I F A < IhUtt • x t AA-1  ) I 
UUUUb2  At  I O A I J 


UUUU32  tNU 

u asTok  ******** 


bObPAUGK AI*I  LL^UlH  bbe> 

U:4(JbL0  CUF.F  ILLK  bFALL  4100 


buoAouixGt.  *bruH(jo»xni  .anust.MT) 

UUUUUb  At AL  T 

UUUUUb  UIMtGbXUU  AMIIXUOI 

UUUUUb  AMI  (JO)  = Mt/XHUSl 

UUUUIU  JJ  s 00  ♦ X 

UUUU1X  AtlUAIN. 

uuuuxx  tr.u 

L ♦••****«  IfjTtb 


SUbAKUbA  All  LtWbtH 
UINlUbLU  COMPXLtA  SPACt 


2U 

Abuo 


uuuuxx 
uuuuxx 
uuuuxx 
u u u u 2 i 

UUUU32 
UUUUbb 
UUUUbb 
U U 0 U 3 I 
UUUUAo 
U U U U J u 

UUUU  ft> 
U U U X U b 
uuuxxu 


SUbMOUI  XUL  imtUUUXiXni  . U»FS>hf’»UUU«NtUUNT»XVAL) 
uinfcHbAUi.  o(b)*C(b.b)iV(i).XMT(lUU) .FbAP( XUUO t 3 ) ,U< 1UU ) 

A = U.  U 
UO  X LL=X«3 

blLL)  = AH  I lUL0UNl+LL-b)*FbKP(WC0UN1-FLL-i»X0X)/n(NU0UI*T'»LL-b) 

L ( LL. 1 ) = X . u 
L ( LL  • 2 > - U(N0UUMI+LL-3) 

L(LLO)  = U<N(.UUNI+LL-3>**2 

X LOMl loot 

U*lL  UUAU(D.L.V) 

a I J A K I 3 = U(  NCUUNT-X  IMLUUU  T - X ) /2  . ♦ V ( A I *U  ( NLUUN1 -X  ) **2 

X /A.)  - U(|«CUUNT)*(  V(X)+(M2)*U(NU0UNI  1/2.+V (3)*0(NL0UNI  >**2/3.  ) 

UO  2 LL  = 1 • 3 

H ( LL ) = xni (LL)*FbKP(LLtXOX|/u(LL) 

C ( LL  * 1 ) = X.U 
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UOUlll 
UUU113 
UUU1 13 
ULUi*J 
UUUi^D 

uuum* 
0 U U 1 A A 

uuui<o 

UUU ib3 

UUUlb / 
U 'J  'J  1 u. 
UUU1  o 
UUU1 fo 
UUU1 f / 

UUUffUb 

U JU*U  / 
UU  U*** 

'JUU^^b 
OL'U**  / 
uuu*3u 

UUU^Oo 

UuU^H3 

UUU*4b 


ULL.^I  = U(LL  I 
ULL*i)  = 

CALL  UUMUlOict V I 

ClH  r 1 fJl.'L 

AHAhll  = U i 2 I • * V t 1 » *V 1 *»  *Jl<i ) /* . *V t :> > *u\ * » ***/ b.  ) - 

1 ui  1 I • t VC  1 )♦  VC*  I • UC  1 )/*.*V  t dl*L*<  A > ***/0.  ) 

LLL  = 1 

cu««  r iniit 

AH  ak  T = (U(  LI  L^2  I • t V u I ♦ v I ^ | *u(LLL^*  !/*•♦¥  I 3 I • UCLLL*’*  I ***/3. 
1 - J ILIL  > • 1 W l 1 ) ♦V  ( c I *U  l LlL  ) /*  . ♦ v ( b I #U  l LLL  > )> /*  • 

A — A ♦ AH AH  | 

if-  l LLL  *bL  »'<LOUN  T -*  ) JO  fo  ( 

LLL  = LLL  ♦ 1 

UU  4 1 =1  • * 
ci(  i ) = M i i ♦ 1 ) 

JO  4 JSltO 

Cl  i • J > = L l I ♦! * J ) 

COM  I If.UL 

t3<3)  = XM  iUL*^)nbhnLLL^2»IUX)/UILLL^<!) 

CCd«l)si«U 
L ( A • * ) = U C LLL** ) 

L C 3 « d ) = lm LlL** ) *»* 

CALL  UUALML3,L«V) 

<*u  ro  ^ 

AVAL  = AKAH  T 1/*  • 0 ♦ AHAp<T^/*«U  ♦ A 

«L  I UNfj 

LNU 


blibHKUGMAIi  LL'^blH 
UNUbLU  CUHH I LtH  bHALt 


O** 

bOOO 
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UuUUUD 
UUUU  1 U 


uuuudi 

UUU UDD 

uuuu:>t 


»•«»•«»***»  I^UmU  *****»*»**•**♦• 

bUbkOUliNL  UUaU ( b «l « V ) 

Ui  NLFVbiU'v  u(  Jl  . C<d«D)f  w<o) 

i ( i c i i • a ) -l ( * • d j » * c c i i • *: » -l ( a « 2 ) )-a»i»5i-u3o)  i 1 1 

V < «r  I = I 1 1»  1 1 « .-> ) -c:  l 4 « 4 M • ( h ( 1 ) -b  ( 2 > > ~ 1 1 1 1 • £ > -C  l 2 ' 4 ) ) • < b t 1 ) -h  < 4 M J / 

1 ( l*IUiO)-L(3,3)  l-ICUii'l-Uiu)  )*IC(1*3)-L(**b)  I I 

V » X ) = hill  - V ( 2 ) * L ( 1 ' <;  ) - V ( 3 ) * L t 1 * J J 

KL  I UHIJ 

cNU 


bUB^KUOHMPi  Ltl.olh 

uuUbLU  cum xllk  bn  all 


hUULTlON  L I ALAL l XNUt 1KUiL I A1 ) 


UUUUUD 

COMMON 

AX  , 

aF. 

AF*  h • 

ALF  Abl  » 

AUXi 

1 

ulja^  • 

bb*K. 

LJLiI'M 

L ui  jX  • 

CbUtiK, 

LDH  • 

* 

U4H * 

ullui  , 

^LLKU. 

fc>  • 

u 0 X • 

t.Hb  • 

3 

UHMriA  f 

bN*  • 

HUl»« 

HUH  • 

KALI  • 

KAH<!  , 

4 

rvHL  • 

• 

F'  0 U.  0 • 

H 0 •*»  4 U N • 

(>S>  1 A»<  1 * 

HblOH* 

D 

UMH  • 

UF-  h t 

uFUtL« 

OL  • 

►'»  » 

K * 

n 

KAM  t 

Kl  t 

KHOF  • 

HhODH  « 

KHtyX  i 

HUN 

UUUUUD 

LOnMON 

DUX  • 

1 AH* 

1 A V . 

1 »-  • 

1 C 1 1<  U . 

X AL  Y A « 

1 

ALAND • 

XimX  % 

Xll2  « 

XMIH* 

XUUb 1 • 

XNUl  t 

A d r HF-  • 

Xb  f HUf 

A S T A K U • 

Xb  ( AH 

UUUUUD 

LOr'iMurj/UUuuLL/ 

hux » 

.<1  • 

fb 

UUUUUD 

UUUUUD 

UUUUUD 
UUUU  ID 
UUUU  lt> 
UUUU 1 / 
UUUU*'* 
UUUUUD 
UUUU/fb 
UUUUDU 

Uuuu 
UUUUD* 
U U l*  I.1  D D 
UUUUb / 
UUUUHD 
UUUUbl 
UUUUb*. 
UUUU  / D 
UUUU  /A 
UUUU  / D 
UUUU// 

uuuiuu 
U u U 1 u b 
UUUU/ 
uuunu 
UUUl ID 

UUU11H 

UUU11D 


ttlMrtUN/ ALL/  H 

LOMMUN/L/  Lb  I AH | « LIAU*  lTAL*  JA^U*  uA^L*  U*  ^bl«  L I AH  * t F Ah  * 
1 lHAti 

iF-  i XImU  • L I • U , i .OH  • XhU.b  I . u • W » bU  10  1 

Lb  > AH  T = LI A1 

CALL  L I ALUN(  Ai4li«  IKU'LI  A1  ) 

LTACAL  = L I A 1 
KL I UNN 

i F l Xi'jU  • L I • u • i ) bU  I a iU 

IF-  l JA^U.LO.  0 ) LMAbT  = LlAl 

AI\|UU  = U • *** 

L ALL  LFALbw{XnUU«lKU»L|Au) 

Lb  I AH  I = LIAU 
AfiUL  = U.*0 

LALL  L T ALUN( ANULt IHUtL I AL ) 

A C = ALU-»U)rtU/Ll«L)/4LUblll.-XUUUl/a(-XIJUU  I 
Hi.  = L 1 AU*  i 1 •/  t 1 .-XJMUUI  ) #*A/ 

L T A l = Hi.*  \ 1 • •»XI>IU)**A£ 

L T ALAL  = L I A1 
JIM  J = i 
KL  i UK|* 

XNUU  = • ti 

IF ( JAll »LMi U I LbTAKI  = tl«l 
CALL  L 1 ALUH  ( XUUU  • U'L)  t L I MU  1 
JA*L  = 1 

t T«i  = XUU*[_  | nU/XNUU 
tTuLAL  = A.IA1 
KL I UKN 

trju 


bUbKKULKAI-i  L{  HG  I M 
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ur»U5tu  conpiLtK  dp  Act 


muo 


UUUUUD 

UUUUUD 


UUUUUD 


UUUUUD 

uuuuud 

UUUUUD 

UUUUUD 
UUUUUD 
UUUUUD 
UUUUUb 
UUUU  U ( 
UUUUll 

uuuuiz 

UUUU1H 
UUUU1D 
UUUUlb 
UUUUZU 
UUUUll 
uuuuzd 
UUUUZH 
UUUUUD 
UUUUUb 
UUUUZ  t 
J L'  U UZ  t 
UUUUDZ 
UUUUDD 
UUUUDD 
UUUUDb 
UUUUD  f 
UUUUH* 

UUUU4D 
UUUU4& 
UUUUH  / 
UUUUD  t 
U J r t b D 
UUUU  / 1 
UUUU  f b 

UUU1UU 
U V U 1 U D 
UUUlli 
UUUi  ID 
UUUlZb 
UUlUb 
UUUI  * / 
U U U 1 D U 


SUbKOUTlNL  El  ACOMXNU*  iPj'tTAl  ) 
KtAL  Ml 


COMMON 

A 1 . 

A 2. 

AF  . 

AFH< 

ALF  Abl 

• AOXt 

1 

bLTAF « 

bbUK, 

C 1 ON  • 

CONI  i 

CDUbP « 

can. 

2 

C4P, 

ullui t 

u2CwU  t 

LF  , 

EOX, 

t-Pb, 

i 

bAMMA . 

OMW  « 

HUN  ■ 

HOF  « 

KAP1  t 

KAPii. 

4 

fNPF  • 

PMb  . 

POwo « 

F’UWibN* 

PM AK1 

. PblOP, 

a 

y Ap . 

OFF  t 

,*F  OtC • 

OL  • 

« 

K . 

b 

WAP. 

KF  . 

kHOF  • 

KHObF , 

KhOX  • 

KON 

COMMON 

DUX  t 

T AF  « 

1 AV. 

TF  . 

1 ZtKUt 

X ALF  A . 

1 

ALAMd« 

xni . 

XI, 2 , 

XNUF  . 

XNUS  1 « 

XNOl . 

2 

ADTPI-  • 

XalFo. 

J^ST  (»KUi 

i X S 1 AH 

COMMUN/UObbCL/ 

MOX  • 

MT  , 

lb 

COMMON/ ALL/ 

p 

COMMON/t/  Lb  1 AK  1 , 

LIAO,  clAL 

i JAiUi 

OALL • 0.  2bl 

. ETAP. 

LTAF  . 

1 IFLAb 

COMMON 

/bLiF  «c/ 

Lirioti* 

LhKbLb 

UIMtNblON  A*JD14) 

IF  1 1 HU  ) DO,lU 
D bM*l  = bMb 
XNUX  = XNUA 
IA  = fAH 
Tbb  = TF 
bO  TO  11 

10  bhwl  = PM* 

XNUX  = ANUK 
l A = lb 
1Gb  = TAP 

11  1 TLIMZ  = DU 
JX  = Z 

dZ  = -1 
K = 0 

FULL  = U.U 

TED T = XNUA/iXNUX  ♦ 1.1 
IF ( XNU.b 1 • I EDT ) GO  TO  DU 
bO  TO  bl 
DU  JX  = 1 
JZ  = 1 

bl  TbP  = ( I bb  + I A )/Z#U 

KHobP  = P*bMl^l/(ttZ.UbD*TbP) 

UbP  = r T /Ntt(JbP 
b = bLAL(ANU) 

Zbl  = bftHflM* lbP**l .7D/ tP*UbP*B> 

L = bUK I l anU I 

Cl  = bUNliA.u  ♦ bb.  7z /b8ZDe*ZSI*Zbi  ) 

CZ  = Z . U*^*  ( X.jUXfl  . ) / ( XfvuX-C*C*  ( XNUX  + l . ) ) 

AKb  = b .001  fUD**  ✓♦C 

CALL  btSdLLt  U.U  • AKb*  AND) 

L - CZ*A.vbi2  l/i  D.Obl  fUD*7*(  ( -.40Z ^D*b9D7 )**dX ) I 
IF  <Z«LT . 1 • U | bU  TO  DU 

LSlMHl  = ^••ZDi*Zbi/(Cl  - 1 * U ) * Al Ub  U I 
DU  CUNllNUt 

t TAi  = U I ah  I 
N = 1 
1 T ILK  = U 
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UUUI  )1 

S’* 

J = U 

UU  U 1 

b[’ 

U - J ♦ 1 

UUUlit 

IF  (LIAl.Ll  .U.  ) LTA1  = ul aULD*.0*»u* 

UUUiH^ 

UHbi  = IlMbU  1 I J ) 

UUUI^H 

(.1  = ONKItl.U  ♦ 4.0*  t a*01*AKO1 **0i»40i  ) ) 

uuuirx; 

DKb^  = AKUi*L 

uuuibb 

LALL  bLOOLL  ( Altbl  . AKO*  . ANa  ) 

UUUlbl 

LO  = AnOT*  1 / 1 ak(,i»  ANO(  1 ) »*uX  ) 

UUUi bo 

(.4  = I 1 ,-U  I *L  1 Al/  l * . *40i*XOi  ) 

UUUI  1 $ 

PaKTI  = L0*l.XP(L4> 

UUUdUl 

PA«[*  = I'tnl  l*L4/t  1 Al 

uuu^u^ 

MJLL1  = PULL1  ♦ PAKT1 

UUU*UH 

FULL*  = FULL*  + P AK  1 * 

UUU^Ub 

IF  < AbOl  Pm*  • 1 1 .OT  .LKHyLO.  ANU.  J.L1  . LlMbtO  1 Oil 

U U U id  1 b 

tlAOLU  = LlAi 

UUUidlb 

L f A 1 = L 1 A 1 - ( F UL  L 1 - (XNUX/L  - L»( ANUAFl . 
1 FULL* 

UUUddb 

FULL1  = U.4 

UUUddb 

FULL*  = U»U 

UUUdd f 

1 I ILK  = i M L«  + 1 

UUUdil 

IF  UlltK.UI  ,1  ILXFi*)  Ou  To  75 

UUUd^H 

oO  TO  / 4 

UUUdiH 

n 

AKllt(b.**4)  LIA1,  LIAULU 

UUU*H* 

■»*  4 

FOKMAT  ( SX«  ' I II  tKA  1 ION  Llflil  UCLLUtUi  ' . 5Xi  ' 
1 ‘LTAULU  =’ iLlfc.fli/l 

UUOd4b 

00  TO  7U 

UUU*H  f 

/ 4 

IF  l Abb  ( 1 .-LTAI/LT  AOLO  1 .L  T .1  .t-b)  00  10  70 

UUU*bH 

OO  TO  Ob 

Uuu*dd 

ft 

LS1AKT  = LlAi 

UUUdib 

70 

KL 1 UKN 

UUU*b  f 

LNU 

OUbPKUGPAr'i  LtNfclH 
UNUOLU  C0HP1LLK  OPALt 


FACTION  bUAL(XNu> 


UUUUUd 

KLAL  mux 

• nr 

UUUUUd 

common 

A i • 

A*. 

AF  . 

AF  FI  • 

ALF  AO  I t 

AUX  t 

1 

ot-TAF  , 

bOUK  , 

CION* 

CONi  • 

COUbP. 

L3P, 

d 

CHHf 

ObLui  . 

UbtKU. 

LF  . 

LOX, 

LKb* 

3 

UAMMA  • 

Ohm  f 

mUN  » 

HOP  « 

K API . 

KAP*  * 

4 

*KF  • 

PMFi  • 

POWU  • 

POwlON. 

POT AKT . 

Kb  1 UK, 

b 

WMPt 

OFF  « 

OFULL  « 

UL  * 

UPF  i 

K t 

b 

kap, 

KF  * 

KHOF  • 

HHOOP. 

KHOX  . 

KUN 

UUUUUd 

COMMON 

aux» 

TAPi 

1 AV» 

TF  . 

ULKUi 

xalh  a , 

1 

ALAMO, 

XM1  . 

xN*  • 

XNUPt 

XNUOT . 

ANU1  , 

2 

AbTKF  • 

XOlPU, 

XSTAKU, 

XOl  AP 

UUUUU* 

LOi'IMON/UOUOLL/' 

MUX  . 

llT  t 

10 

UUUUU* 

UUUOUH 
UUUUUO 
UUUUU  7 
UUU011 

ouuu^a 

uuuu*o 

UUUU* / 
UIHJU  AO 


LOI'iMUIM/ ALL/  P 

khop  = khuop 
KF1  = P.  1/bMOK 
K API  = MOX/KMOX 

I ION  = LiON*  (UZtKO*U  .UUUi  ) **  (POWU+1.  I /P**POWlON 
XUU  = KH»IIbN/»UXLHO*0,uU01) 
i/UO  = XUU  ♦ KF1/KAP1 
UUU  = AfUNl(l.U»VUO) 

OION  = 2./O.»u*E.KU»U*tK0*lA.0*(UUU*UUU-XUU*XUU) 


- *.0*(UUO**O  - 


104 


J 


UUUUHb 

X xuu**3M/iuuu  - x u i 
IFlUlbNI  **3»3 

UUUU3U 

2 

UlbN  = 0.0 

UuUU^i 

bU  TU  4 

UUUU3C 

3 

UlbN  = SbK 1 ( UlbN ) 

UUUUDb 

4 

LUNT 1NUL 

UUUU5b 

USh  = (X.U  - XNU ) *U< 

000063 

ttCAL  = 3**k  » t Ub»-  ♦ U 

uuuu/u 

Kt 1 UKN 

uuuu/u 

LNU 

5U0KKUbKAM 

LtNGIM 

UNUStU  CUHH1LLK  SHALL 

UUUUU* 

FUNCTION  OKOU1  INI  ) 
if i^l.bT.OU ) bO  1U  l 

uuoous 

bU  TU  i X 1 * « 3 • 4 • b « b f 

ouuo ft 

1 

X 7«*tt.29  . 

X*4b«4/»4d#49tbU)  • N 1 
bKUOT  = 3.O31/039' 

ooou  / 4 

KE 1 UKN 

UUUU  / 4 

2 

bKUUT  = ^»UX5t)bbb/ 

UUUU  f b 

KL 1 UKN 

UUUU  / b 

5 

bKUUT  = 10.1/346614 

u juiuu 

KL 1 UKN 

UUU1UU 

4 

bKUUT  = 13.3*369194 

ouoiu* 

Kt 1 UKN 

uouiu* 

6 

bHUOT  * lb.4/ub5U03 

UUU1UH 

KE 1 UKN 

UUU1UH 

6 

bKUUT  = 19. 61363631 

UUUXUo 

KL 1 UKN 

UUUlUb 

7 

bKUUT  = ** • /6006436 

UUU11U 

Kt 1 UKN 

UUU11U 

S 

bKUUI  = *3  • 90367206 

UUUXX* 

KE  1 UKN 

UUUXX* 

V 

bKUOT  = **•  04662633 

OUUXX4 

KE 1 UKN 

000X14 

10 

bKUUT  = 3*. 16967991 

000116 

Kt 1 UKN 

UUUXXb 

11 

bKUUT  = 33.33*40735 

OUUl*U 

Kt 1 UKN 

uuux^u 

12 

bKUOT  = 30 « 4 /4  / fob* 3 

UUUi  <LtL 

Kt 1 UKN 

U UMl'tt 

15 

bKUOT  = 41.61/094*1 

0001*4 

KE 1 UKN 

OOUl*4 

14 

bKUOT  = 44.7593X900 

UUUXXb 

Kt l UKN 

UUUXXb 

lb 

bKUUT  = 4/. 90X46009 

000130 

Kt IUKN 

uooi3u 

16 

bKUUT  = 31 . U4333bX0 

uuux 3* 

KE 1 UKN 

U0013* 

17 

bKUOT  = 34.10333364 

U0U134 

Kt  f UKN 

000134 

16 

bKUUT  = 3 ( • 32 / 32344 

UUUXXb 

KE IUKN 

UUUXXb 

1* 

bKUUT  = bO. 46943703 

UUU140 

KE 1 UKN 

j.45 

HlUO 
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Ouu l*u 

20 

UHUUI  = 

bA.fellAOb/U 

ouoi*»c 

«t  1 UKN 

UUUIH^ 

21 

bKUOT  = 

ob.  70A22b 13 

UUU iHH 

kl 1 ukn 

UUU 1HH 

22 

bKUUT  = 

b7.070U710H 

DUUiHb 

Kt  » UKN 

UUUlHb 

di 

UHUUI  - 

/O.  0AOB902A 

UUU1DU 

Kt  » UKN 

•JUUl'jU 

2H 

bKUUT  = 

< O.  X /Ub990O 

Jl'U  13* 

Kt 1 UKN 

JUUX3* 

23 

UKUUT  = 

/7. A2UHB71B 

JUU1JH 

«t 1 UKN 

J U U X 3** 

db 

bKUUT  - 

t)<  .4b220971 

JUU  X3b 

Kt  * UKN 

JUUlbb 

27 

dkuut  = 

OO.bOHUl9HH 

JUUXdU 

Kt 1 UKN 

JUVlbV 

do 

bKUUT  = 

BO.  7<4b/b7lH 

JU  0 ib^ 

KL 1 UKN 

JOUlb^ 

29 

bKUUT  = 

71 . 00/00 H20 

JUUlbH 

KL  | UKN 

JUt/lb* 

3u 

bKUU?  = 

70. 0272A1B1 

JUUlbb 

Kfc.  | yKN 

JUUibb 

31 

bKUOT  = 

70.1 /u70U / 3 

JUUX  f U 

KL 1 UKN 

JUUX fv 

3d 

UKUOI  - 

lUi.AX2bbl02 

Juux  r* 

KL  1 UKN 

JUUX 

33 

bKUOT  = 

XU1*  .HOHAbO  /9 

juux  f+ 

KL 1 UKN 

JUUX  '*+ 

AH 

bKUU?  = 

lU/,37bUbA2b 

JUUX  ^ b 

Kt 1 UKN 

JUUX  /b 

33 

bKUUl  = 

110  , 13/  7t>H  te 

JUU^UU 

KL 1 UKN 

JUU2UU 

3b 

dWUUl  = 

llA,0/9HHU0b 

JUU^Uc 

KL I UKN 

JUU2U* 

3r 

bHUUl  = 

11 / . U2 1 12 19U 

JUU^UH 

Kt 1 UKN 

JUUkiU* 

30 

bKUUT  = 

ICO . 102796AA 

JUU^Ub 

KL 1 UKN 

JUU^Ub 

A9 

bKUUT  - 

120.  AUHH7UH<* 

JUU21U 

KL ( UKN 

JUU*lU 

HO 

bKUUT  = 

12b.HHblAB/U 

JUU21* 

KL 1 UKN 

JUUkil* 

HI 

bKUUT  = 

127.30/00020 

JUU^XH 

KL 1 UKN 

JUU^XH 

H2 

bKUUl  = 

102,  /29HfaHA9 

JUU^Xb 

KL 1 UKN 

JUUkJXb 

Hi 

bKUUT  = 

lAO,B/1122Ab 

JUU*<iU 

KL 1 UKN 

JUUk'j'U 

HH 

bKUOT  = 

1A7.U12/7/A9 

JUU*<!* 

KL 1 UKN 

JUtf*** 

Hb 

bKUUT  = 

AH2.10HH29bb 

JUU**H 

KL 1 UKN 

Jl'U24* 

Hb 

bKUUT  = 

lHD.27bU77AH 

JUUkf^b 

KL 1 UKN 

JUU^kib 

H7 

bKUUT  = 

1HB . H A 772bb2 

JUU*3U 

KL 1 UKN 

JUUkfAU 

HB 

bKUUT  = 

101 , 0 /7A 7XbA 

JUU*3* 

KL 1 UKN 

JUU^^H 

4* 

t3KUU  1 = 13H,  f*lUl*02 

KL  l UKN 

JUU^OH 

JUU^^b 

luu^d  / 

30 

bU 

DHUUT  = 13  f • Ob2bt)^HU 
Kt i UKN 

xni  s m 

JUU^HU 

)UU**0 

XK  = 2 .0*3, 1 ,o  ♦ H#0*ANl) 

ONU0T  = Xm*i.  l4lD^2fa3H.  ( 1 . 0 ♦ 1.U/CH.0*XN1I 

1 12.U/XK**3 

JUU^DO 

KL 1 UKN 

LNU 

SUttKKU&HAM  LLNGIH  4bb 

JNUbtU  COH^ILLK  bKALt  53U0 

bUBKUUHNL  BLbbtL  1 x . T . ANb  1 

JUUUUD 

UIMLNblON  ANb(2l 

r r 

ON  NtlUKN 

L 

ANb (11=  0 2LK0 

(X) 

L 

1 

A N b ( 2 ) = o ONt 

( T 1 

tuuuuo 

w 

IF (X-3.  l2»2.4 

IUUUU  f 

2 

2=A*X/*. 

lUUUll 

IF (XI2U.2U.3U 

»uuuic 

20 

ANb ( ll=x. 

»UUU  iH 

00  10  b 

IUUU1H 

30 

CONUNOt 

3.U/XK  + 


IUUU14 

IUUU2  7 
IUUU3U 


ttl\lb(  1 ) = ( < M ( . UU021 *2-. 003*444  ) *Z  + . U 4444 7*  ) *2-.31b3bbb  » *2  + 
*i.2b3b2uai  *2-2.24***o7>*2+i. 
on  iu  b 

4 4=3. /X 


IUUU31 

IUUU43 

t U U U Db 
IUUU  fU 


F2LKO=l((><.UUU1447b*2-.0UU72e0b)*2+.UU137237l*/-,UUUU*312l*2 
*-.  0033274  l*z-.  7 7L-b  1*2  + . 7*7b643b 
TH2t KO=x+M(((. 00013330*2-. 0002*3331*2-. UUU 34 X23)*2+.U02b2373)*4 
♦-.U0UU3*34)*2-,041b63*7>*2-.7o33*Olb 
ANbl i »=F2Lrto*C0b<  Hi2LKO)/bOKT( X » 
b k-  i r-3. io.o.iu 


iuuu  7 3 
IUUU  73 

iuuu  ft 

'UUU  7 7 
IUU1UU 
‘UU1UU 

'UU113 

'UU113 


o 2=r*r/*. 

IF  l T l*0.>tU.5U 
4 u AWb(2)=U. 

30  10  ttU 
30  CON  I I NUt 

ANb\2 > = ( t »» U ,0000110**2-. 0 0031761 ) *2+ . 0044331* ) *2- . 0 3*342B* ) *4 
»+. 210*33 70 1*2 -.3624*^03 | *2*.3)*Y 
60  KtlUKN 
lu  2 = 3. /T 


'UUllb 

1UU13U 

00*43 
UU133 
U U 1 33 


F ONL=  (((('-. UUU20U33*2+.UU113b33>  *2-.  0024*311  >*Z-F. 000171031*2 
*+.Ulb3*bb7  )*2+1.3bC-b)  *2  + • 7976643b 
THUNt=T+l <(<(-. OOU2*lbb*<;+.OUU790  24)*2*-. 00074  34  61*2*. UUb3  76  7*1*2 
*♦ . UU0U3b3U 1 *2*. 124**bl2 1 *2-2 . 35bl*43 
ANb ( 2 • =F  UNL*CUb ( TNONt I /SilK  f ( 7) 

Kt I UKN 
LNU 

******  0UTPU1  ****** 


.UBKKOGKAH  LFNb  TH  273 

iNUbtU  C CJ71K  I Lt-K  bFACt  3bU0 
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subKuui  lut.  ou  TPUT  t rt  1 1 lo*  i 

UUUU4  KtAL  KAPHKAP2iKPF  iMOXiMl  .KAPlUiKAP2UtKPFU 


UUUUH 

COMMON 

«1, 

A<i  i 

aF  i 

AP  M « 

ALF  AbT  i 

AUX  , 

1 

OL  TAJ-  » 

060K  , 

CifeNi 

CONI  % 

CbUbPi 

CiP« 

LHP  . 

ULLUl  , 

U/t-KUl 

LP  • 

EOX. 

LPbt 

OAMMA  « 

bMw  « 

HUN  i 

HUP* 

K API  i 

HAP*  , 

4 

*PP‘  • 

PM#«  • 

POwU  i 

PUwiOlJ  » 

PbTAKI  , 

Pb 1 UP  « 

ID 

«Ap  « 

OPP  t 

oF  ULL i 

OL  • 

OPF  i 

H t 

b 

KMP  « 

HP  « 

hFlOF  • 

HHDDP  • 

KFiUX  i 

HUN 

UUUUH 

common 

iUX  « 

Tap* 

t AV  • 

TP  . 

T2LK0 i 

X ALP  a « 

X 

ALAMO * 

XNi  • 

XN2  • 

XijUH  , 

XNUb  1 i 

XNU 1 « 

2 

AST  PJ-  , 

XilPU, 

XST  AKU i 

Xb  1 AP 

UUUUH 

COr'ittUh/UUUbLL/ 

MUX  « 

MTi 

Ti> 

UUUU«4 

COMMON/ all/ 

P 

UUUU4 

COMMON/  X 1'*  1 / 

K 

UUUUH 

CCP'MUU  /KW  Uf 

bLIA« 

KHUM  1 

UM 

UUUU1* 

COMMON/bU^i-/  1 h U ( ^>  ) • 

• XLU1UI3) 

1 XIJUPU  t 3 ) i 

1 PMWU ( 3 ) « 

*»LU  (31. 

KHOAU ( 5 ) * 

AuXU(3)i  E0XU13), 

1 APU ( 3 ) • 

^ LlbrjU(4» 

• PUwlOUt 

ill  POhUUtbli  CUNFUtbti  1 AVu  t 3 ) i 

3 *PP  U l 3 i « 

K AP1U ( 3 ) 

i K«P2o'3) 

* X f 4 X U ( 3 J « 

XU2UC3I 

t XN3U l $ ) » 

a CbUBPut3)i  XLAMbUl3li  OAMMAU  ( 3 ) • AFPiutil*  LPbu(3)i  YUI3) 

UUUU4  COMMON/Ou  I t ALFAIlbi3li  bl  bl“i  A 1 < 5 i 3 > i UbAKItbi3l.  MOUEb 

UUUU4  U 1 ntNblOl'i  bl<(3UI,  PY(3U),  bK  A ( 30  ) * PIA13U) 

UIMjut  1,0  10  < lUi*Uibi2)  iM 

UUU13  2 LONTINUt 

U U U 1 3 mKI  lEtbillU) 

UUU17  11U  F OKMA  I I Ini  • 13X  i 13b  Fut-L  u A T A IS,//) 

0UU21  «K  1 TL  1 b 1 1 11  ) I2LKU 

UUU2b  ill  FOKMAI  (2X.C1H  PKGP  INiliAL  1 Lf.P  lbiFb.li/) 

UUU3U  WKi tt t b i ll* ) X ALF  A 

UUU33  112  FUKMa  l < HA  • 2bH  Mabb  FKaO  OF  OXiUlAtK  lb  .FS.2i/J 

U U 0 3 / wHi  I H b » H3  ) OFUEL 

auutt  113  F OKMA  t t 2x  « 29FI  FILAl  OF  PYKULTS1S  OF  FULL  ISiFb.li/l 
uUU Ab  Ah  1 Tt l 6 . 114  ) KHQF 

UUUD3  114  FOKMAI  (2XH9H  UENblTY  OF  FUEL  Ib<F3.<'</l 

UUUbb  IaKI  lL<btl!3)  aF 

ItJUbi  lib  FOKMAI  I2X.3UFI  AKKhLNlub  FKEQ  KAO  OF  FUEL  lbiE9.3i/l 
UULK4  WKiTt(b.llb)  LF 

UUU/1  lib  FOKMAI (2X.29M  AO  1 1 V A f 1 UN  ENEKbY  OF  FULL  1 b i F 7 , U , // / / ) 

UUU / 3 KL I UKN 

U U 0 / 3 b CONTINUE 

U U U 7 3 WKITLtbillM  1UX 

OUlui  117  FQKFA  ' l2X'lbH  0X101 ZEK  NUMbtK, I3»3H  UAtA.//) 

JU1U3  WhltLtb.llO)  OLU(IUX) 

uum  na  f okma  r 1 2x  < <:ah  latent  mlai  of  uxiu  ibiFb.ii/i 

UU113  WKllElb.119)  KFIOXUtlUX) 

UU121  119  FOKMAI 12XH9H  UENblTY  OF  UXIU  ISiF3.2i/1 

UU123  WKlfElb'lCO)  AOXOtlOX) 

U0131  121)  FOKMAI  (2X|30F1  AHhhLMo.  KLO  FAC  OF  uxiu  Ib.t*.3«/I 

UU133  AK 1 t L < b • 1F1  I LOXOtlOX) 

UU141  121  FOKMAI  <2Xi29FI  ACIIVAIIUN  ENEKbY  OF  UXIU  lb«F7.U./l 

UU143  wHlTElbtlZZ)  TAPUdOx) 

UUlbl  122  FOKMAI (2X.17M  AP  FLAME  TEMP  IbiFb.Ui/) 

UU133  WK1 IE < bl 1Z3 ) ClONUtlOX) 

UUlbl  123  FOKMaT (2AI29M  0X1U  IGNITION  ULLAY  PAKam  IbiFb.U,/) 

UUlb3  »Kl  tL  ( b.  1F4  ) F'OWt)U(IUX) 

U U 1 / 1 124  F OKMAT ( 2X • 41H  OIAMLTek  EXP  IN  0X1U  lONlT  ULLAT  TEKM  Ib.Fb.3./> 
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001/3  wmrttb«l*3)  FOKlGUTluX) 

U0201  12S  K)KH»U^»0»M  FKLSS  LXF  IN  ox  IGN  ULLMT  TLKM  XS.Fb.3./) 

UUiJUi  WK1  1 1 T b<  12b  ) kFFO  T 10X ) 

U0211  12b  FOKMAT 12X.31H  FKIMAKT  flAML  KATL  CONSTANT  XS.Fs.l,/) 

0 0213  .KlTtlbili:/)  KAFlUlXOXI 

0U221  127  FOKMAT  (2Ai2bM  AH  F C AML  KaIL  CONSTANT  lS.Fb.ii/) 

00223  WK1TLT6.120)  XNIU(IOX) 

00231  12B  F OKMA I ( 2X • Obr»  KLACIIUN  OkOLH  OF  FKIMAKT  FLAML  iS.FS.l./l 

UUiiO  wH  1 1 1 ( b .12b  ) XN20TX0X) 

0 0 2X1  12b  F OKMA  I ( 2X  . 30H  KLACIION  OKULK  OF  AF  FLAML  Xa.FA.l./) 

002X3  wN 1 IL l b . iOU ) CSUBFUTIOA) 

0 0 2 31  130  FOKMAT  (2X.21F1  AVL  ML  A 1 L AH AC  1 T T lb.Fb.3./) 

00233  mK ITLtb.lAl)  XLAMbU(IOX) 

U U 2 b 1 131  FOKMaT  (2X.3.5H  THLKMAL  COnUUCUVIIT  OF  OASLS  1S.F8.3./) 

U 0 2 b 3 hK 1 TL  T b « 132  ) GAKMAU(iUA) 

0 U 2 / X X 32  F OKMA T < 2X • 23M  UIFFUSIOn  HAKAMLTLK  13tLb.2./> 

002  /3  *K1 ILT6.133)  WOOLS 


UU3UU 

133 

FOKMAT (2A»xbH  NUMbLK  OF  MOOES  IS.XX./) 

UU5UC 

00  ISA  X ■ =X i WOOLS 

UU3U3 

h* 1 T L l b . 133 ) SXGMAX ( 1 1 • I jX  ) 

UU313 

166 

F0KMATl2X.2bH  STANO  OcV  uF  TFiIS  MOUL  13.Fb.2l 

JUi  ID 

*K 1 TL 1 b . 13b ) UBAK  1 l I T . XOX  ) 

UU323 

16b 

FOKMaT (2X.37M  WT  MLAN  OI«M  OF  0 X10  XN  THIS  WOOL  IS.Fb.U) 

UU32  ( 

WKX  ILIbiXOM  ALF  AX  ( I T « XOX  ) 

UU33  ( 

167 

FOKMAT  (2x«3bFI  MASS  F K AC  oF  THIS  OAXO  I ia  THXS  MODE  1S.F7.3./I 

UU341 

134 

CON  I X NUL 

UU5H0 

Kt 1 OKN 

UUilH 

10 

bKtIU  = K 

UU34b 

FTTH)  = H 

UUADU 

HT  A ( K 1 = H | ( A.  1 *14 . 7 

uu3D2 

*jK  A ( K ) = tJK  ( K ) /2  ,bA 

UU^DH 

KM«X  = K 

UU334 

K = K + 1 

U U ADb 

KL 1 OKN 

UU3t>b 

20 

CONTXNUL 

U IJ  D D b 

«K1 TL I b . 102 ) 

UU362 

KKlTLTb.XOX) 

0U3b  f 

00  bb  H = 1|  KMX  X 

UU3  ( £ 

*KlTL)b.lU3)  HI  (K)  1 F T A ( K ) . BK ( K > , OK AIK) 

UG4U  i 

b4* 

CONTINUE 

UUHl  1 

KL 1 UKN 

UU412 

102 

FOKMAT  <ttA.4MFKLS.3X.4MFKES.7X.2FlbK.bA.2HBK) 

UUHU 

103 

TOKMxl  OX*F7.2.2X.F7.X,2(2X.f9.A)) 

uum2 

104 

FOKMA  1 ( ttX  * AHA  1 MS.  3X  . 4MFSX  A « SX  , bFtCM/sLL  . 5X  > bH I N/ SEC / 1 

uu**l<i 

END 

.UBHKUGHAM  LLNGIH  10/7 

.NOSCO  COWHXLLK  SHACL  1300 


xfob  U333 
Abb  033X 


Abb  Ujb/ 
Abb  03bt> 
Abb  03/1 


Abb  OX/t) 


I 
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FULL  UAI«  lb  Sample  Problems 

PKUP  lNillAL  T 1 1-1 H lb  SUU.U 
nflbb  F KAL  Ot-  UXiUUt  « lb  .70 

«L«t  ot-  PlKOLlSib  Uf-  tULL  ib  200. U 
otNaiir  ot-  t- ut.L  lb  1.2' 

AKKHtNlOb  FFIlu  t-AL  ot-  FOf.L  1S2.'UUL  + U.i 
ACilVAIluN  t-NLKOT  OF  FULL  lb  1S0UU 

UA1012LFC  NUFIpLK  1 L)AI« 

LAU.N1  ML  A I OF  UA  10  lb  U.O 
ULNblTY  OF  0X10  lb  1,-XD 

ArinMLNIUb  PnLU  FAC  OF  oaiu  1 S4 . 0 U Ut ♦ 0 s 
AL11VAI10N  LNtKOT  OF  OAIU  lb  22UUU 
AM  FLAFlL  I LnP  lb  140U 
0A1U  I ONI  1 1 ON  UtLAT  PAH  AFI  lb  l'XU 

01AI*H.TtK  LAP  IN  0 A 1 [)  1 ON  1 | UtLAT  ILKFi  lb  ,60U 

PfXLSb  LAP  IN  OX  1 014  ULLAY  ILKFI  ib  .7pU 

PpiFlApY  FLAFlL  mail  CONblANi  lb  20.0 

AP  FLAFit  KAIl  CUNblANI  is  2. ISO 

KtALlIUN  OKOlk  OF  PHIFiaat  FLAMt  lb  1.6 

KLALT10N  OKOLK  UF  AP  FLAFlL  IS  1.6 

A V L Ht  A I LAPAL1IT  lb  .AUU 

IHLKFIAL  LONOUCIIVIIT  OF  GAbtS  lb  .UOubU 

O1FFUS10N  PAKAFILltH  lb  /.bUfUb 

NUFItiLK  OF  FIOOLb  lb  1 

blANU  UtV  OF  T Ml  S FIOUL  IS  1 • 0 U 

Wl  F>F  AN  OlAil  OF  OX  1 0 1 N THlb  FlOUt  ib  20U 

Fi  Abb  P K A L UF  THIS  0X10  IN  IHIb  FiOUL  lb  .700 


►'KUKELLAN  I OLNdill  = 1 « bB  U 1 
V WLO^L  F*AL  UF  WAlu  = .b031 

L FKUM  VK  s t*UALHU**N  lb  .344S3 
JHtbbUHt  lb  29. H The  UX1U12LH  bt  1 NO  LONblUtKLu  lb  1 


U2c.HU 

KAIL 

1b 

XNU5T 

ALF  AbT 

KrtUbH 

TF 

ttt  TaF 

XbTHU 

X3TKF 

XbT  mP 

XblU 

199.02 

. 0993 

79  b 

• 00  31 

,7000 

1 .bOUl 

2393 

1.0000 

33.07 

23.9b 

191.90 

lb. 3b 

199.0b 

• 099b 

79  b 

• Oubl 

. 7U00 

l.oeoi 

2393 

1.0000 

33  , bO 

23.93 

191.9b 

lb. 30 

199, 0 3 

,U99b 

79b 

.0031 

. 7000 

1 . t»oui 

2393 

1.0000 

33,00 

23.93 

191.93 

10.39 

199 , Ob 

. l)99b 

79b 

• OUbl 

, 7000 

1 . bOU 1 

2393 

1 . uoou 

33,09 

23.93 

191 ,9b 

10.39 

199,00 

« 0 993 

/9b 

• oObl 

, 7000 

1,9001 

2390 

1 .0000 

33.09 

23.93 

191.92 

10.39 

199,911 

. U99b 

7 9b 

.oubl 

. 7000 

l.bbol 

2393 

1 .UOOU 

33,90 

23.93 

191.91 

10.39 

199.91 

,U99b 

79b 

.oubl 

, 7000 

1 .btJOl 

2393 

1 .0000 

33, 9U 

23.93 

191.40 

10.39 

199,9b 

. U99b 

79b 

.oubl 

. 70UU 

l,b»01 

2393 

1 .UOOU 

33.91 

23.93 

191,39 

10.39 

199. 99 

. U99b 

79b 

• oobl 

. 7000 

1 . t»eui 

2393 

1.0000 

33.91 

23.93 

191 .30 

10.39 

199.90 

,U99b 

79  b 

• OObl 

. 7000 

l.bOUl 

2593 

1.0000 

33.92 

23.99 

191.37 

10.90 

199.90 

. U992 

79  b 

• OQbl 

, 7000 

1 . beui 

2593 

1 .uuuu 

33.92 

23,99 

191 . 33 

10.90 

199,9* 

.0992 

79b 

.0031 

, 7uOU 

1 .OOul 

2393 

1.0000 

33,9b 

23.99 

191 , 39 

10.90 

2U0.  Ol 

.0992 

79b 

.oobl 

. 7000 

1 .bttOl 

2393 

1 • ouuo 

33.9b 

23,99 

191.33 

10.90 

200. U2 

.0992 

793 

.0031 

. 7000 

1 .bOOl 

2393 

l.UOUU 

33.99 

23.99 

191.32 

10.90 

2UC . UH 

.0992 

793 

.oobl 

. 7000 

l.bOOl 

2593 

l .uoou 

33.99 

23.99 

191.31 

10.90 

200  • Ut 

.0992 

793 

.0031 

. 7000 

l.bOUl 

2393 

1 .0000 

33.93 

23.99 

191,30 

10.90 

*00  . 0 t 

,099c 

793 

.0031 

. 7000 

1 .0601 

2393 

1 . uoou 

33,93 

23.99 

191,29 

10.91 

*00. U* 

.0992 

793 

.oobl 

. 7000 

1 .bttOl 

2393 

1 ,0000 

33.96 

23.99 

191,27 

10.91 

*UU. 1U 

.0992 

795 

.0031 

. 7000 

1 .6001 

2593 

1.0000 

33.9b 

23.93 

191.2b 

10.91 

*JU. 1* 

.0992 

795 

.0031 

. 7000 

1 . bOUl 

2393 

1. uoou 

33.9  7 

23.93 

191.23 

10.91 

*U0,14 

.0992 

793 

.0031 

. 7000 

1 ,bt>0l 

2393 

1.0000 

33.97 

23.93 

191.29 

10.91 

*U0.  lb 

.099c 

793 

.oubl 

.7000 

l.bOUl 

2393 

l.uoou 

33.90 

23.93 

191.23 

10.91 

*U0. 1 • 

.0992 

793 

.0031 

, 7000 

1 .bOOl 

2393 

1 . ouuo 

33.90 

23,93 

191,22 

10.91 

*00  . 10 

.0992 

793 

.oQbl 

, 7000 

l.bOUl 

2393 

1 .0000 

33,99 

23.93 

191.21 

10.92 

*UU • *U 

VA^Ut 

.0992 
UF  KAIL 

793 

1 IV  1 L U K A L = 

.0J31 
• 1933*7 

. 700U 

1 .bOOl 

2393 

1 .uuuu 

33.99 

23.93 

191.19 

10.92 

lUKL  lb 

52.0 

I ML 

OXIulZt-K 

ttt 1 NO  CUNOlULKLU  lb 

1 

(J2LK0 

KATE 

10 

xivub  1 

ALF AbT 

KhOSH 

TF 

ttt  lAF 

XOIHU 

XOTHF 

XMmH 

Xb  1 u 

199.02 

.11  72 

739 

.oubl 

,7uuu 

1 .ObOl 

2343 

.0239 

02.34 

10 . 12 

Ob . 40 

10.35 

199.03 

. 11  72 

7 39 

.0031 

. 7UO0 

1 .obUi 

2343 

• 0230 

02.34 

1U  . 12 

(J8.HU 

10.36 

199,03 

.1172 

739 

.0031 

. 7UUU 

l.bOUl 

2343 

.6237 

02.33 

10.12 

5b. HU 

10.36 

199.00 

.1172 

739 

.0031 

, 700U 

l.ObOl 

2343 

• 6236 

02. 3b 

10.12 

Ob  . 39 

10,36 

199.03 

. 11  72 

739 

.0031 

. 7 U u 0 

1 .ObOl 

2343 

.0236 

02. 3b 

10.12 

Ob  , 39 

18. 3t 

199, 9U 

.1172 

739 

.0031 

. 7000 

1 .ObOl 

2343 

.0233 

02.37 

10.12 

00,39 

10.36 

199.91 

.1172 

739 

.0031 

. 7000 

1.0001 

2543 

• 0234 

02.30 

10.12 

00.30 

10,37 

199.93 

.11  72 

739 

.0031 

. 7000 

l.bOUl 

2343 

.0233 

02. 3b 

10.12 

00.30 

10,3? 

19^,99 

.11  72 

739 

.0031 

, 7000 

1 .bbOl 

2543 

.0232 

02.39 

1U.12 

00,30 

10.37 

199.9b 

• 11  72 

739 

.00*1 

. 7000 

1 .0  01)1 

2343 
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